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19 .  ABSTRACT 


EXECUTIVE  SUMMARY 


This  Phase  I  Small  Business  Innovative  Research  (SBIR)  contract  produced  a  Design 
Manual,  a  guidebook  intended  for  use  by  designers  of  two-phase  fluid  flow  and  heat  transfer 
systems  for  space  applications.  Design  guides  are  available  for  one  gravity  applications  of 
two-phase  systems,  but  none  of  these  are  currently  applicable  to  microgravity  thermal 
management  systems.  This  contract  produced  two  documents:  a  Final  Report 
(AL-TR-89-026)  and  a  Design  Manual  (AL-TR  89-027).  This  Design  Manual  will  assist 
in  the  confident  design  of  two— phase  system  components.  It  provides  for  quick  scoping 
calculations  and  more  detailed  component  calculations.  Although  not  covered  under  this 
contract,  these  design  methods  can  be  integrated  with  existing  compuiwi  >-odes  for  transient, 
system— level  performance  predictions  of  thermal  management  systems.  The  Design  Manual 
will  save  considerable  time  and  effort  in  the  development  of  modelling  tools  by  the  Air  Force. 
As  part  of  this  contract,  recommended  areas  for  follow-on  mechanistic  model  design  analysis 
were  identified.  These  efforts  would  expand  on  the  topics  already  presented  in  the  Design 
Manual  and  permit  the  analysis  of  an  entire  thermal  management  system. 

Future  Air  Force  space  missions  will  require  power  levels  and  heat  transport  rates  at  one 
to  two  orders  of  magnitude  larger  than  present  systems  (10  to  100  kW).  Since  two-phase 
systems  are  a  candidate  approach  to  meeting  heat  transport  capabilities  for  those  higher  power 
levels,  future  thermal  management  systems  will  involve  vapor— liquid  flows.  As  the  missions 
are  specified  and  two— phase  systems  come  to  the  drawing  board,  a  framework  of  design 
methods  common  to  the  designers  (primarily  the  Air  Force,  NASA  and  the  aerospace  industry) 
will  be  valuable.  In  this  Design  Manual,  the  plethora  of  multiphase  flow  literature,  which  can 
easily  overwhelm  an  analyst  who  is  a  novice  to  two— phase  flows,  is  distilled  into  practical 
design  methods  which  can  be  conveniently  used.  This  Design  Manual  is  intended  to  be  a 
working  document  illustrating  selected  topics  of  basic  two-phase  flow  regime,  pressure  drop, 
and  heat  transfer  methods  across  the  range  from  earth  gravity  to  microgravity  conditions. 

The  focus  and  organization  of  this  document  is  at  the  component  level.  Typical 
two-phase  transport  systems  will  involve  components  such  as  transport  lines  (pipes), 
evaporators,  condensers,  pumps,  separators,  and  radiators.  Separate  sections  of  this  Design 
Manual  will  present  design  methods  for  pipes,  evaporators  and  condensers  at  a  level  consistent 
with  a  Phase  I  SBIR  contract.  Alternate  approaches  can  be  used  in  the  design  of  each  of  these 
types  of  components,  and  methods  which  are  of  general  application  to  different  approaches  are 
presented  in  the  subsections. 

The  sections  of  this  Design  Manual  present  the  methods  at  two  levels: 

•  Preliminary  sizing  calculations  can  be  performed  for  the  components 
using  simple  calculations  and  dimensionless  maps  presented  in  the  text, 
and 

•  Detailed  equations  are  presented  so  that  the  methods  can  be  incorporated 
into  design  software  for  detailed  component  analysis  or  in  system 
analysis  codes. 
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Each  section  includes  information  which: 


A.  Describes  the  background  and  discusses  the  limitations  of  the  use  of  the 
methods  presented, 

B.  Presents  the  design  information  (including  maps  or  graphs)  for  simple, 
scoping  design  calculations, 

C.  Describes  in  detail  the  analysis  methods  (which  would  be  used  to 
develop  detailed,  computerized  analysis  routines),  and 

D.  Validates  the  design  methods  against  available  microgravity  data. 

The  scoping  design  information  for  each  topic  is  also  useful  because  the  dimensionless 
parameters  illustrate  the  governing  physics  and  the  relationships  between  the  key  design 
parameters. 

This  document  considers  steady— state  modelling  of  the  two— phase  flow.  Transient 
two-phase  flows  can  be  addressed  by  implementing  these  methods  into  transient  analysis 
codes.  The  methods  incorporated  in  this  Design  Manual  are  all  based  upon  mechanistic 
analysis  methods  to  the  extent  possible.  There  are  several  advantages  to  the  mechanistic 
analysis  approach  used  here  versus  approaches  based  upon  empirical  correlations.  Use  of 
mechanistic  methods: 

•  Ensures  that  the  design  methods  scale  appropriately  with  size,  pressure 
(gas  density),  and  fluid  properties. 

•  Uncertainties  in  the  design  methods  can  be  quantified  by  specifying 
ranges  for  key  modelling  parameters. 

•  The  methods  can  be  readily  improved  and  refined  as  new  experiments 
are  performed  to  verify  the  models,  without  disturbing  the  basic 
framework  of  the  models.* 

The  ability  to  quantify  uncertainties  in  modelling  the  component  behavior  is  particularly  useful 
in  the  system  design  context.  The  key  uncertainties  in  the  system  operation  can  be  identified. 
Therefore,  good  designs  can  be  made  and  the  range  of  uncertainty,  potential  risk,  or  margin  in 
the  design  can  be  identified  in  order  to  minimize  the  risk. 

By  analogy  with  analysis  codes  used  for  reactor  safety  calculations  by  the  nuclear 
industry,  the  component  design  methods  presented  here  are  suitable  for  incorporation  into 
system  codes  which  can  be  used  for  steady— state  or  transient  analysis  of  space  thermal 
management  systems.  For  example,  the  analysis  methods  are  compatible  with  the  ATHENA 
computer  code  presently  being  used  by  the  Air  Force  for  preliminary  system  design  and  with 
the  SIND  A  '85/FLUINT  code  being  used  by  NASA  to  design  thermal  management  systems  for 
the  Space  Station. 


*Two— phase,  microgravity  research  is  being  performed  by  a  number  of  investigators. 
The  methods  presented  in  this  Design  Manual  have  immediate  application  in  the  design  of  the 
experiments  and  the  interpretation  of  the  data  produced  by  these  experiments. 
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Tn  this  Design  Manual,  the  four  topics  included  are: 

(1)  Flow  regimes  in  adiabatic  pipes,  unifying  methods  to  determine  flow 
regimes  at  earth  gravity  and  microgravity, 

(2)  Pressure  drops  in  adiabatic  pipes,  using  two-phase  flow  models  for  four 
flow  regimes  of  interest  (bubbly,  slug,  annular,  and  stratified)  where  the 
methods  are  simple,  easy— to— use,  and  sufficiently  accurate  for  design 
purposes  (compared  with  alternate,  detailed  analysis  methods), 

(3)  Condensation  heat  transfer  in  a  condenser  controlled  by  the  shear 
between  the  vapor  and  liquid  phases  (in  a  pipe  of  uniform,  circular 
cross-section),  and 

(4)  Boiling  heat  transfer  in  an  evaporator  dominated  by  forced  convection 
heat  transfer  (in  a  pipe  of  uniform,  circular  cross-section). 


The  documentation  for  each  topic  includes  comparisons  of  the  methods  with  available 
microgravity  dati .  Thus,  the  manual  also  serves  as  a  summary  and  review  of  available 
microgravity  research  results. 
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1  TRANSPORT  LINES:  Cylindrical  Pipes  with  Wetted  Walls 


Most  thermal  transport  systems  for  spacecraft  will  include  two— phase  flow  in  the 
transport  line  between  the  evaporator  and  condenser  components  or  in  the  radiator 
components.  Some  systems  will  permit  two— phase  flow  at  the  evaporator  exit.  Others  will 
have  single— phase  flow  at  the  evaporator  exit,  but  will  have  two-phase  flow  in  the  transport 
line,  because  pressure  drops  and  heat  losses  from  the  transport  line  will  lead  to  condensate 
formation  between  the  evaporator  and  the  condenser. 

The  purpose  of  this  section  of  the  Design  Manual  is  to  provide  methods  for  determining 
the  two— phase  flow  regime  and  for  calculating  the  pressure  drop  and  void  fraction  in  straight 
segments  of  piping  without  phase  change.  Knowledge  of  the  flow  regime  is  crucial  to  the 
calculation  of  the  pressure  drop  and  void  fraction.  Void  fraction  is  especially  important  when 
calculating  transient  system  behavior.  The  basic  data  which  must  be  supplied  by  the  designer 
are  the  fluid  properties,  the  pipe  size,  the  acceleration  (magnitude  and  direction),  and  the 
volumetric  flow  rate  of  the  gas  and  the  liquid  phases. 

The  methods  in  this  section  apply  to  adiabatic  flow  in  pipes  of  uniform,  circular 
cross-section.  Alternate  geometries  such  as  square,  rectangular,  or  triangular  channels,  which 
may  be  used  to  enhance  heat  transfer,  must  be  treated  separately. 

This  section  assumes  that  the  liquid  phase  will  wet  the  material  of  the  wall.  If  the  liquid 
does  not  wet  the  wall  surface  (such  as  mercury  in  a  glass  tube  for  example),  then  the  flow 
regimes  and  pressure  drop  models  may  be  different  than  described  here.  This  situation  may  be 
treated  later. 


1.1  Flow  Regimes 

This  section  deals  with  the  determination  of  the  flow  regime.  A  key  feature  of  the 
design  methods  recommended  and  used  here  is  that  the  methods  are  applicable  to  a  wide  range 
of  acceleration  levels,  including  earth  gravity  and  microgravity,  and  all  "inclinations".  The 
upper  limit  on  the  acceleration  level  would  be  the  point  where  fluid  properties  vary 
significantly  across  the  flow  channel  due  to  hydrostatic  pressures.  Once  the  flow  regime  is 
selected  from  four  possible  choices,  the  pressure  drop  specific  to  that  regime  can  then  be 
calculated  (see  Section  1.2). 

Section  1.1.1  provides  a  brief  introduction  to  the  concept  of  flow  regimes.  It  also 
identifies  the  information  which  the  designer  must  have  in  hand  in  order  to  determine  the  flow 
regime  and  the  pressure  drop. 

Section  1.1.2  presents  dimensionless  design  maps  for  the  flow  regime  determination. 

This  section  is  set  up  so  that  the  user  does  not  need  to  know  the  details  of  the  theory.  The 
designer  approaches  the  flow  regime  calculations  with  knowledge  of  the  phase  flow  rates, 
physical  geometry,  and  fluid  conditions.  Simple  flow  regime  maps  arc  then  used  where  the 
designer  can  calculate  a  few  dimensionless  parameters,  plot  a  coordinate  on  a  map,  and 
determine  the  flow  regime. 


1 


Section  1.1.3  presents  the  basic  equations  behind  each  of  three  flow  regime  boundaries 
which  determine  the  four  flow  regimes.  The  equations  are  presented  in  sufficient  detail  that 
the  designer  could  recreate  the  dimensionless  design  maps  presented  here  or  generate  design 
maps  in  dimensional  coordinates  (such  as  phase  velocities)  if  desired.  Examples  of 
dimensional  design  maps  are  provided. 

Finally,  Section  1.1.4  compares  the  flow  regime  maps  for  the  recommended  design 
methods  against  selected  data  from  microgravity  experiments. 


1.1.1  Introduction 

Figure  1.1  illustrates  the  four  primary  flow  patterns  of  interest  for  two— phase  flow  in 
pipes.  The  patterns  include: 

•  Stratified, 

•  Slug, 

•  Annular,  and 

•  Bubbly. 

These  patterns  are  referred  to  as  "flow  regimes." 

In  the  stratified  regime  the  gas  and  liquid  phases  are  segregated  by  body  forces  due  to 
gravity.  Liquid  flows  near  the  bottom  of  the  pipe  and  the  gas  flows  near  the  top.  This  regime 
does  not  occur  in  a  "vertical"  pipe  on  the  ground.  (On  a  spacecraft,  the  acceleration  vector 
would  have  to  have  a  component  perpendicular  to  the  pipe  for  this  flow  regime  to  occur.)  In 
the  slug  flow  regime,  the  vapor  forms  long  bubbles  so  that  the  vapor  phase  is  not  contiguous  in 
the  axial  (flow)  direction  but  is  separated  by  liquid  slugs.  The  vapor  bubbles  may  be  several 
pipe  diameters  in  length  and  may  be  either  asymmetrical  or  symmetrical  with  respect  to  the 
central  axis  of  the  pipe.  In  the  annular  flow  regime,  the  vapor  flows  as  a  continuous  gas  core, 
and  the  liquid  flows  as  a  thin  film  on  the  wall  of  the  pipe.  In  the  bubbly  flow  regime,  the 
vapor  flows  as  discrete  bubbles  within  the  liquid  phase. 

The  key  to  determining  the  flow  regime  is  to  first  identify  the  dominant  forces  acting  on 
the  fluid.  Balances  between  the  dominant  forces  are  then  used  to  describe  the  boundaries  (the 
transitions)  between  the  four  flow  patterns  described  above.  Three  flow  regime  transitions  are 
needed  to  bound  the  four  regimes: 

•  Slug  to  stratified, 

•  Slug  to  bubbly,  and 

•  Slug  to  annular. 

These  three  transitions  depend  upon  different  balances  of  forces,  as  described  in  the 
subsequent  discussions. 
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HORIZONTAL  GAS  LIQUID  FLOWS 


BU8BLY  SLUG  ANNULAR 

VERTICAL  GAS  LIQUID  FLOWS 


Figure  1.1.  ILLUSTRATING  FLOW  PATTERNS  (REGIMES)  OF  INTEREST 
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For  this  Design  Manual,  the  model  used  for  each  of  these  flow  regime  transitions  is 
mechanistic.  In  contrast  to  empirical  correlations,  the  mechanistic  approach  considers  the 
governing  physics  for  each  transition,  therefore  the  mechanistic  models  should  scale  well  with 
all  acceleration  magnitudes  (g— levels)  and  vectors.  There  are  three  basic  forces  and  two  basic 
force  balances  to  consider  in  the  determination  of  the  flow  regimes: 

•  Balances  of  buoyancy  and  inertia  forces,  and 

•  Balances  of  surface  tension  and  inertia  forces. 

References  which  consider  the  balances  of  buoyancy  and  inertia  forces  include  Taitel  and 
Dukler  (1976),  Bamea  (1986),  and  Dukler  (1988).  References  which  consider  the  balances  ol 
surface  tension  and  inertia  forces  include  Lee  (1987)  and  Kachnik  et  al.  (1987). 

There  is  a  tendency  in  current  literature  to  draw  a  distinction  between  flow  regime 
models  for  earth  gravity  and  for  microgravity.  Based  upon  the  experimental  evidence 
presented  thus  far,  the  models  developed  for  earth-gravity  conditions  reduce  satisfactorily  to 
the  appropriate  limits  for  microgravity  conditions  without  resorting  to  separate  models. 
Examples  are  given  for  selected  flow  regime  transitions  by  Reddy  Karri  (1988)  and  Revar.kar 
(1989).  In  this  Design  Manual,  it  will  be  shown  that  all  of  the  significant  regime  transitions  at 
earth  gravity  reduce  to  microgravity  conditions,  hence  a  single  set  of  flow  regime  transitions 
can  be  used  at  all  acceleration  levels  and  orientations.  Table  1.1  summarizes  the  models 
recommended  for  the  Design  Manual. 

The  recommended  models  have  been  extensively  compared  against  experimental  data  at 
earth  gravity  in  various  flow  orientations  ranging  from  vertical  upward  to  vertical  downward. 

It  will  be  shown  that  the  models  for  each  orientation  reduce  to  the  same  model  under 
microgravity  conditions  and  compare  well  with  the  available  microgravity  data. 


Table  l.l.  RECOMMENDED  MODELS  FOR  FLOW  REGIME  TRANSITIONS 

TRANSITION 
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RANGE  OF 
UNCERTAINTY  FOR 
KEY  PARAMETER 
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!  i 
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TURBULENCE  '14  0. 725  to  3. 7 

COEFFICIENT  <k) 

CRITICAL  VOID  0.45  0.40  to  0.5: 

FRACTION  (a..) 

BUBBLE  D I SnilBimON  i  1.2  j  1.0  to  1.2  : 

'  1  _ 1 

SLUG-TO- ANNULAR 

BARNEA  (1986) 

''NEUTRAL 

STABILITY’ 

INTERFACIAL  SHEAR 

(Vf.,) 

[1  +  ?5(l-a)l 

1  CO 

[I  +  75(1-0)1 

or 

"PHYSICAL 

BLOCKAGE" 


1.1.2  Dimensionless  Design  Maps 

This  section  illustrates  how  the  flow  regime  can  be  quickly  determined  from  known 
conditions.  The  information  required  in  order  to  determine  the  flow  regime  includes: 


Fluid  Properties 

p  i  Liquid  phase  density 
pg  Gas  phase  density 
Pi  Liquid  phase  viscosity 
|ig  Gas  phase  viscosity 
a  Surface  tension 


Geometry 
D  Pipe  diameter 


Acceleration 

a  Magnitude  of  acceleration 

0  Direction  (angle  between  acceleration  vector  and  the  flow 
direction) 

Flow  Rates 

j  i  Liquid  phase  volumetric  flux  (Q  /A) 
jg  Gas  phase  volumetric  flux  (Qg/A) 


Each  of  the  three  flow  regime  transitions  can  be  represented  on  dimensionless  design 
maps.  The  dimensionless  coordinates  for  these  regime  maps  can  be  determined  from  the 
information  listed  above. 

The  order  of  examining  the  transitions  is  important  to  the  determination  of  the  flow 
regime,  since  certain  regimes  supersede  others.  The  order  should  be: 

•  Slug-to-Stratified  Flow  Regime  Transition, 

•  Slug-to-Bubbly  Flow  Regime  Transition,  and 

•  Slug— to— Annular  Flow  Regime  Transition. 

Figure  1.2  illustrates  the  logic  for  the  flow  regime  selection.  The  dimensionless  design  maps 
for  flow  regimes  are  presented  and  discussed  in  that  order  below. 

Slug— to— Stratified  Flow  Regime  Transition.  To  determine  first  whether  or  not  the  flow 
regime  is  stratified,  the  proposed  design  point  is  mapped  on  a  dimensionless  plot  of  gas  Froude 
number  Fg  versus  the  Martinelli  parameter  X  as  shown  in  Figures  1.3a  and  1.3b.  The 
definitions  of  these  dimensionless  parameters  used  in  the  maps  are: 
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INPUTS 

*  FLUID  PROPERTIES 
(density, viscosity, 
surface  tension) 

*  GEOMETRY 
(pipe  dianeter) 

»  GRAVITY 

(magnitude,  vector) 

*  ROW  RATES 
(each  phase) 


IF  NOT 
SATISFIED 


T 


(STOP) 
REGIME  IS 
SLUG 


Figure  1.2.  LOGIC  FOR  FLOW  REGIME  DETERMINATION 


GAS  FROUDE  NUMBER, 


b)  Smooth  vapor/liquid  interface 


Fieurc  1  3  DIMENSIONLESS  DESIGN  MAP  FOR  SLUG-TO-STRATIFIED  FLOW  REGIME 
TRANSITION  (TURBULENT  FLOW) 


(1-1) 


where 

£ 

Y 

a 

Pi 

Pg 

Ji 

Jd 

0 


Pg°'5jS 


Fg  = 


[aD(p  i  -  pg)]0-5 


rr4c,i 

[“IT 

fPiQi.il 
l  Hi  J 

-n  1 

j 

(4CJ 

_[t> 

Pg  J 

-m 

hH 

a(Pi  -  pg)cos0 

.[#] 

Pgfile 

PgJ 

-m 

(1-2) 


(1-3) 


is  the  dimensionless  Froude  number  for  the  gas  phase, 

is  the  dimensionless  Martinelli  parameter, 

is  the  dimensionless  pipe  inclination  parameter 

is  the  magnitude  of  the  acceleration  (m/s2), 

is  the  gas  phase  density  (kg/m3), 

is  the  liquid  phase  density  (kg/m3), 

is  the  superficial  velocity  of  the  liquid  phase  (m/s), 

is  the  superficial  velocity  of  the  g?s  phase  (m/s), 

is  the  pipe  diameter  (m),  and 

is  the  angle  betweer  the  flow  direction  and  the 

acceleration  vector  (°)  (e.g.  0°  is  colinear,  90°  is 

perpendicular,  and  180°  is  opposing,  corresponding 

to  downflow,  horizontal,  and  upflow  on  the  ground). 


Table  1.2  indicates  how  to  determine  the  values  of  the  friction  factor  coefficients  C\,  Cg,  n  and 
m  to  use  in  Equations  1—2  and  1-3. 


The  procedure  in  using  Figures  1.3a  and  1.3b  is  to  first  calculate  the  three  dimensionless 
quantities  in  Equations  1—1,  1—2,  and  1—3.  These  quantities  depend  upon  the  superficial  phase 
velocities  (or  volumetric  flow  rates;  ji  and  jg  of  each  phase  and  the  fluid  properties.  The 
desired  operating  point  is  plotted  on  the  design  map  by  the  coordinates  (X,  F„).  Then  the  line 
of  constant  Y  for  the  proposed  design  conditions  is  located.  If  the  design  point  falls  below  the 
line  for  the  value  of  Y  (Equation  1—3),  then  the  flow  regime  is  stratified.  If  the  design  point 
falls  above  the  line  for  the  value  of  Y,  then  the  slug— to— bubbly  flow  regime  transition  is 
assessed  as  described  below.  Note  also  that  if  the  acceleration  level  is  very  small,  for  example 
at  microgravity  conditions  where  a  =  KHg,  then  Y  =  0. 

The  difference  between  Figures  1.3a  and  1.3b  represents  the  range  of  uncertainty  in 
modelling  the  regime  transition.  (The  uncertainty  specifically  lies  in  the  interfacial  friction 
between  the  gas  and  liquid  phases,  as  described  in  Section  1.1.3.)  The  proposed  operating 
point  should  be  located  on  each  of  these  figures.  If  both  figures  indicate  that  the  flow  is 
stratified,  then  the  flow  is  probably  stratified.  If  Figure  1.3a  indicates  that  the  flow  is  stratified 
and  Figure  1.3b  does  not,  then  the  flow  regime  may  be  uncertain. 
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Tabic  1.2.  SUMMARY  OF  PARAMETER  VALUES  FOR 
FRICTION  FACTOR  CALCULATION 


PHASE 

REGIME 

CRITERION 

c 

g 

m 

VAPOR 

LAMINAR 

Reg,  - 

l  Hg 

<  1500 

16 

1.0 

TURBULENT 

Regs  > 

1500 

0.046 

0.2 

PHASE 

REGIME 

CRITERION 

C1 

□ 

LIQUID 

LAMINAR 

Re  is  = 

ptDj  r 

Hi 

<  1500 

16 

1.0 

TURBULENT 

Re  i,  >  1500 

0.046 

0.2 

Slug— to-Bubblv  Flow  Regime  Transition.  The  transition  from  the  slug  to  the  bubbly 
flow  regime  should  be  examined  as  the  second  step.  The  slug-to— bubbly  flow  regime 
transition  consists  of  two  parts: 

•  A  transition  based  upon  fluid  turbulence,  and 

•  A  transition  based  upon  a  critical  void  fraction. 


and  both  of  these  criteria  need  to  be  checked  to  determine  the  flow  regime. 

Figure  1.4  is  the  dimensionless  design  map  for  the  transition  based  upon  turbulence.  To 
determine  the  flow  regime  a  dimensionless  Weber  number 


Wej  = 


rP  iDj2- 

pP(jg  +  ji)2 

a 

a 

(1-1) 


and  a  dimensionless  Bond  number 


Bo 


aD2(p  l  -  Pg) 


are  calculated  for  the  proposed  design  conditions,  where 


Pi  is  the  liquid  phase  density  (kg/m3) 
pg  is  the  gas  phase  density  (kg/m3) 

D  is  the  pipe  diameter  (m), 
o  is  the  surface  tension  (N/m), 

jg  is  the  superficial  velocity  of  the  gas  phase  (m/s),  and 

j  i  is  the  superficial  velocity  of  the  liquid  phase  (m/s). 


(1-5) 
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Figure  1.4.  DIMENSIONLESS  DESIGN  MAP  FOR  SLUG-TO-BUBBLY  FLOW  REGIME 
TRANSITION  (TURBULENCE  CRITERION) 


The  point  with  coordinates  (Bo,  Wej)  is  then  located  on  Figure  1.4.  If  the  point  is  above  the 
transition  line  in  the  figure,  then  the  regime  may  be  bubbly,  but  the  second  (void  fraction) 
transition  criterion  must  be  checked  first.  The  three  lines  in  Figure  1.4  illustrate  the  range  of 
uncertainty  for  the  slug— to— bubbly  flow  regime  transition.  If  the  point  is  below  all  three 
transition  lines,  then  the  slug— to— annular  flow  regime  transition  must  be  checked,  as  described 
below. 

Figure  1.5  illustrates  the  second  slug-to— bubbly  flow  regime  transition  criterion,  based 
upon  a  critical  void  fraction.  For  this  "ansition.  the  dimensionless  parameters  are: 

M  =  (ji/jg)  l-o  i 


(ac(pi  -  pglp-25 


P~  =  tpt/pg)  0-8) 

where  M  is  the  dimensionless  superficial  velocity  ratio, 

Kg  is  the  dimensionless  gas  Froude  number, 
p*  is  the  dimensionless  density  ratio, 

and  the  other  parameters  are  defined  above. 
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Figure  1.5.  DIMENSIONLESS  DESIGN  MAP  FOR  SLUG-TO-BEBBLY  ROW  REGIME 
TRANSITION  (VOID  FRACTION  CRITERION) 


The  proposed  design  point  is  plotted  in  Figure  1.5  at  the  coordinates  (1.4/K,  vP*.  Ml.  If 
the  design  point  is  below  the  transidon  line,  then  the  flow  regime  is  bubbly.  If  the  design 
point  is  above  this  line,  then  the  slug— to— annular  flow  regime  transition  needs  to  be  evaluated. 
In  Figure  1.5,  the  shaded  region  represents  the  range  of  uncertainty  in  the  flow  regime 
determination  by  this  criterion. 

Slug— to—Annular  Flow  Regime  Transition.  The  slug-to— annular  flow  regime  transition 
involves  two  of  the  same  dimensionless  groups  as  the  stratified-to— slug  transition;  namely,  the 
parameters  X  and  Y.  The  proposed  design  point  is  plotted  on  the  graph  of  Figure  l.b.  for  the 
coordinates  (X.Y).  If  the  location  of  the  design  point  is  beiow  the  curve  .n  the  figure,  tttett  :::e 
flow  regime  is  annular.  If  the  location  of  the  design  point  is  above  the  curve  in  the  figure, 
then  the  flow  regime  is  slug.  The  two  lines  in  the  figure  illustrate  the  range  of  uncertainty  tor 
this  transition. 

It  should  also  be  noted  that  for  microgravity  conditions  where  a  =  !(H'g  me  value  ot  Y  > 
approximateiv  0,  so  that  onlv  the  X  coordinate  :s  unDonani.  That  is.  the  uoint  cun  'e  ciottcc 
astX.O). 


1.1.3  Equations  for  Flow  Regime  Transitions 

This  section  describes  the  detailed  analytical  equations  used  for  each  of  the  flow  regime 
transitions.  Additional  design  maps  in  dimensional  coordinates  can  be  generated  if  desired 
using  the  equations  presented  in  this  section. 

Slug-to-Stratified  Regime  Transition.  This  flow  regime  transition  involves  the 
simultaneous  solution  of  two  equations.  One  equation  is  a  dimensionless  two— phase 
momentum  equation  for  a  separated  flow  and  the  other  is  a  dimensionless  transition  criterion. 


II 


Figure  1.6.  DIMENSIONLESS  DESIGN  MAP  FOR  SLUG-TO-ANNULAR  FLOW  REGIME 
TRANSITION  (TURBULENT  FLOW) 


The  dimensionless  transition  criterion  in  this  model  is  based  upon  a  balance  between  gas 
inertia  and  buoyancy  in  the  gas  and  liquid  phases.  Figure  1.7  illustrates  the  concept.  If  the 
gas  flow  rate  is  large  enough,  the  pressure  drop  created  by  the  flow  of  the  gas  over  a  wave  (the 
Bernoulli  effect)  is  sufficient  to  lift  the  waves  in  a  stratified  flow  to  the  top  of  the  pipe.  A 
critical  gas  flow  rate  is  necessary  for  this  condition  to  occur. 

To  map  this  transition  in  the  design  method,  the  dimensionless  liquid  height  h*  is  varied 
across  the  range  from  0  to  1.  For  a  given  value  of  the  dimensionless  liquid  height,  the 
criterion  for  the  critical  gas  velocity  is  determined  by: 


F„  =  II- 


A/ 


'0.5 


J j - -  i - 

'  U„* 1  idA, 


o  -  - 


;  /dh* 


_ J  \ 


where  Table  1.3  defines  the  other  dimensionless  geometric  parameters  which  are  also  based 
upon  the  liquid  height,  h*. 

The  critical  gas  velocity  at  the  transition  is  then  determined  from  the  definition  of  the 
gas  Froude  number  Fg  by: 

jg  =  Fg[aD(p  |  —  pg)/pg]°  • 5  (1-10) 
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Figure  1.7.  ILLUSTRATING  FORCE  BALANCE  FOR  SLUG-TO-STRATIFIED  ROW  REGIME 
TRANSITION 


Table  1.1 

3. 

DIMENSIONLESS  GEOMETRY  PARAMETERS  FOR  SLUG-TO- 
STRATIFIED  FLOW  REGIME  TRANSITION 

h*  = 

hj  /D,  (Figure 

1.7) 

> 

n 

■w— 

{* 

-  cos~1(2h:* 

-  1)  +  (2h*  - 

1)[1  -  (2h*  -  l)2]i} 

A  *  - 
Ag  ~ 

1 

jcos-1(2h*  - 

-  1)  -  (2h*  -  1) [ 1  -  (2h*  -  l)2]i J 

Si*  = 

[jt  -  cos-1(2h* 

-  D] 

c  * 

=  [cos~1(2h*  - 

1)] 

S,*  =  [1  - 

(2h*  -  l)2]i  = 

(dA^*/dh*) 

i*  =  :'k/4al*'|  =  l/u-cc) 


U2*  =  f^/4Ag*1  =  1/a 


4  \* 


D,*  = 


V  = 


4A0‘ 


[V  *  Si 
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where 


a 

D 

Pi 

Pg 


is  the  gas  Froude  number  from  Equation  1—9, 
is  the  magnitude  of  the  acceleration  (m/s2), 
is  the  pipe  diameter  (m), 
is  the  liquid  phase  density  (kg/m3),  and 
is  the  gas  phase  density  (kg/m3). 


The  critical  gas  velocity  is  then  used  to  calculate  a  Reynolds  number  for  the  gas  phase 
by: 


Regs  =  (PgDjg/M-g)  (1-11) 

where 


pg  is  the  gas  phase  density  (kg/m3) 

L  is  the  superficial  velocity  of  the  gas  phase  (m/s), 
D  is  the  pipe  diameter  (m),  and 
Pg  is  the  gas  phase  viscosity  (kg/m— s). 

Then  a  gas  phase  friction  factor 


fwg  =  CgRegs-«  (1-12) 

is  calculated,  where,  if  Regs  >  1500  then 

Cg  =  0.046  and  (1-1 3a) 

m  =  0.2  (1—1 3b) 

and  if  Regs  <  1500  then 

Cg  =  16  and  (l-14a) 

m  =  1.0.  (1— 14b) 


After  the  friction  factor  is  calculated,  a  dimensionless  parameter  defined  by: 


Y 


a(p  i  -  pg)cos9 


j^j*g 


(1-15) 


is  calculated,  where 

0  is  the  inclination  of  the  acceleration  vector 

relative  to  the  direction  of  the  flow  in  the  pipe  (°). 

The  parameter  Y  is  then  used  in  a  dimensionless  two— phase  momentum  equation  to 
determine  a  dimensionless  parameter  related  to  the  liquid  flow  rate: 
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X2  = 


(ue*Dg.r(ue*)2 


|4 

LAg  LxwgJ  l/M  . 


-  4Y 


[Ui*D]*]  _n(U  i*)2(S  i*/A,*) 


(1-16) 


where  the  dimensionless  geometry  parameters  (starred  parameters)  are  given  in  Table  1.3  as  a 
function  of  the  dimensionless  liquid  level  and 


f;/fWf 


is  a  constant  with  a  recommended  value  of  10, 
has  two  possible  values  (0.2  or  1.0). 


Thus,  two  values  of  X2  are  calculated,  one  with  each  value  of  n  assumed. 

The  liquid  velocity  can  then  be  calculated  from  the  definition  of  the  parameter  X2  by: 


Ji  = 


X2 


f4C , 


PgDj..g 

Pg 


r4£ 


PjDl-n 

L  M-  iJ 


V 


1/(2— n) 


(1-17) 


where: 

Ci  is  0.046  with  n  =  0.2  and  is  16  with  n  =  1.0 

and  two  values  of  the  superficial  liquid  velocity  are  then  calculated,  one  with  each  set  of 
values  for  the  parameters  Ci  and  n.  The  superficial  liquid  velocity  at  the  transition  is  the 
minimum  of  the  two  values. 

Equations  1-9  and  1—16,  with  the  equations  in  Table  1.3,  produce  the  parameters  Fg  and 
X  for  dimensionless  flow  regime  plots  such  as  Figures  1.3a  and  1.3b.  This  assumes  that 
constant  values  of  Y  and  f  i  /fw g  are  selected  for  use  in  Equation  1-16.  The  range  of 
recommended  values  for  f  i  /fw g  is  from  1  to  10.  The  value  of  1  represents  a  smooth  gas-liquid 
interface.  The  value  of  10  represents  a  rough  (wavy)  gas— liquid  interface. 

Equations  1—10  and  1—17  convert  the  dimensionless  parameters  into  superficial 
velocities  (volumetric  fluxes)  for  each  phase.  For  each  value  of  h*  selected  for  the 
calculations,  a  pair  of  values  ( j  i,  jg)  can  be  determined.  Figures  1.8a,  b,  and  c  show  the 
calculated  slug— to— stratified  regime  transition  for  one  set  of  conditions  at  earth  gravity.  The 
different  figures  show  the  results  at  angles  of  180°,  90°,  and  20°.  Compared  with  the  results 
where  the  flow  direction  and  the  acceleration  vector  are  90°  apart  (corresponding  to  a 
horizontal  pipe  on  the  ground),  the  range  of  the  stratified  flow  regime  shrinks  as  the  flow 
direction  and  the  acceleration  vector  oppose  each  other  more  (corresponding  to  upflow  on  the 
ground).  The  stratified  regime  expands  as  the  flow  direction  and  the  acceleration  vector 
become  colinear  (corresponding  to  downflow  on  the  ground). 

Figure  1.9  shows  the  effect  of  the  magnitude  of  the  acceleration  level  on  the  results.  The 
stratified  region  also  diminishes  with  reduced  gravity  for  flow  regime  maps  in  these 
coordinates.  In  fact,  the  stratified  flow  regime  is  of  little  interest  at  acceleration  levels  of 
l(Hg  to  KHg  or  lower. 
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Note:  Example  of  map  for  specific  conditions  -  not  for  general  use 
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Figure  1.8.  EFFECT  OF  PIPE  INCLINATION  ON  FLOW  REGIME  MAP  AT  EARTH  GRAVITY 
(SUPERFICIAL  VELOCITY  COORDINATES) 
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Note:  Example  of  map  for  specific  conditions  -  not  for  general  use 
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Figure  1.9.  EFFECT  OF  GRAVITY  ON  FLOW  REGIME  MAP  (SUPERFICIAL  VELOCITY 
COORDINATES) 
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Slue— to— Bubbly  Flow  Regime  Transition.  This  flow  regime  transition  involves  two 
parts.  The  first  part  considers  a  balance  between  the  inertia  of  the  two  phase  flow  and 
buoyancy  (Figure  1.10).  A  limiting  criterion  in  this  condition  is  the  maximum  packing  of 
bubbles  in  the  pipe  —  the  void  fraction  can  only  be  so  large  before  the  bubbles  touch  and 
coalesce  into  large  bubbles  in  a  slug  flow  pattern. 

For  the  first  transition  criterion,  Bamea  (1988)  suggests  an  equation  which  has  the 
dimensionless  form: 


where 


Wej  = 


'  1.2k  ‘ 
fw!0-4. 


L67Bq0.835 


(1-18) 


fw  1  is  a  friction  factor  with  a  value  of  =  0.005,  and 

k  is  a  constant  in  the  range  of  0.725  to  3.7.  (Taitel  and  Dukler  (1982) 

recommend  a  value  of  k  =  1.14). 


Figure  1.4  plots  this  dimensionless  equation  for  the  flow  regime  transition.  The  range  of 
uncertainty  in  the  value  of  k  from  0.725  to  3.7  is  shown  in  the  design  map. 


Using  the  definition  of  the  two-phase  Weber  number  in  Equation  1-4,  it  can  be 
rewritten  to  solve  for  the  superficial  liquid  velocity 


ji 


We  j  a 


Ip  id 


0.5 


Jg 


(1-19) 


Figure  1.10.  ILLUSTRATING  FORCE  BALANCE  FOR  SLUG-TO-BUBBLY  FLOW  REGIME 
TRANSITION 
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Equation  1—19  can  be  used  to  map  the  transition  as  (jj,  jg)  pairs  by  selecting  a  range  of  values 
of  j»  and  calculating  the  corresponding  j  \  values,  using  the  Weber  number  from  Equation 
1-18. 


The  second  criterion  for  the  slug— to— bubbly  transition  comes  from  a  basic  drift— flux 
equation  (Zuber— Findlay,  1965)  for  bubbly  flow 


= 


C0(jg  +  ji)  +  1-41 


ac(pi  -  pg) 


Pi2 


0.25 


(1-20) 


where 

ac  is  a  critical  value  of  the  void  fraction  for  maximum  packing  of  the 

bubbles.  Dukler  (1988)  recommends  a  value  of  ac  =  0.45,  though  other 
recommended  values  range  from  0.4  to  0.52  (Tai tel— Dukler,  1976). 

C0  is  a  constant  typically  near  1.2. 


Equation  1—20  can  be  written  in  dimensionless  form  as 


Jl' 

r  i  ] 

i 

1.41 

Pg' 

jg. 

i 

C0ac 

Kg 

.Pi. 

(1-21) 


The  design  map  in  Figure  1.5  plots  the  dimensionless  regime  transition  of  Equation  1—21  using 
the  parameters  defined  by  Equations  1-6  to  1—8.  The  design  map  uses  the  recommended 
value  C0  =  1.2  and  shows  the  unceuaimy  due  to  the  variation  in  the  critical  void  fraction  from 
0.4  to  0.52. 


Equation  1—20  can  be  written  to  solve  for  the  superficial  liquid  velocity  as 


ji  = 


Similarly  to  Equation  1—19, 
appropriate  range  of  values. 


-  1 


Cnar 


Jg 


1 .411 


l  C0  J 


aa(Pi-  pg)l 


0.25 


(1-22) 


values  of  (j),  jg)  pairs  can  be  determined  by  varying  jg  across  an 


The  transition  from  slug  to-  bubbly  flow  regimes  is  mapped  by  selecting  the  pair  of 
coordinates  with  the  higher  vaiue  of  jj  for  a  given  value  of  jg.  Said  another  way,  for  the  flow 
regime  map,  the  transition  criterion  with  a  higher  liquid  velocity  at  a  given  gas  velocity 
(Equation  1—19  or  1—22)  takes  precedence  in  the  flow  regime  transition. 

Figures  1.8  and  1.9  also  plot  the  slug— to— bubbly  flow  regime  on  regime  maps  in 
superficial  velocity  coordinates.  Figures  1.8a,  b,  and  c  show  that  pipe  inclination  has  very 
little  effect  on  the  transition.  This  is  confirmed  by  data  at  earth  gravity  conditions  (Bamea, 
1986).  Figure  1.9  shows  that  as  the  magnitude  of  the  acceleration  is  reduced,  the  portion  of 
the  transition  criterion  related  to  turbulence  (Equation  1-19)  disappears  and  the  portion  of  the 
transition  related  to  void  fraction  (Equation  1-22)  remains. 
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Slue— to— Annular  Flow  Regime  Transition.  This  regime  transition,  proposed  by  Bamea 
(1986),  also  consists  of  two  parts.  Both  describe  blockage  of  an  annular  core,  but  by  different 
mechanisms.  As  Bamea  explains,  the  first  mechanism  ("neutral  stability")  is  an  instability  of 
the  annular  flow  configuration,  and  the  second  occurs  when  the  liquid  film  is  large  enough  to 
cause  blockage  as  a  result  of  flow  in  the  liquid  film. 

The  first  part  of  the  transition  involves  the  simultaneous  solution  of  two  equations.  In 
dimensionless  form,  the  two  equations  are: 

r  (fi/fwg)  if1! 

Y  =  - - - X2  (1-23) 

a2-5(l  -  a)j  [(1  -  a)3 


[2 -1.5(1 -a)] 

Y  =  - X2  (1-24) 

(1  -cx)3[l-  1.5(1  -  a)] 


Equation  1—23  is  the  two— phase  void  fraction  equation  for  the  annular  flow  regime  (Wallis, 
1969).  (See  Section  1.2.3.)  Equation  1—24  is  the  transition  criterion  derived  by  Bamea 
(1986). 

If  the  void  fraction  is  varied  between  values  of  unity  and  zero,  the  Equations  1—23  and 
1-24  map  out  coordinates  of  (X,  Y)  as  shown  by  the  upper  portion  of  the  curve  with  positive 
values  of  Y  in  the  design  map  of  Figure  1.6. 

The  second  transition  criterion  involves  using  a  critical  void  fraction  ac  =  0.76  in 
Equation  1-23.  Coordinates  in  (X,Y)  pairs  can  be  mapped  out  when  this  substitution  is  made. 
This  leads  to  the  lower  portion  of  the  transition  criterion  with  negative  values  of  Y  in  the 
design  map  of  Figure  1.6. 

The  parameter  (fj/fw»)  in  Equation  1—23  represents  the  interfacial  friction  between  the 
gas  and  the  liquid  phases.  For  the  purpose  of  these  design  maps,  the  value  proposed  by  Wallis 
(1969)  based  on  tests  at  earth  gravity  is  used: 

(fi/fwg)  =  [1  +75(1  -a)]  (1-25) 

which  is  also  used  by  Dukler  (1988)  for  microgravity  comparisons.  Various  models  for  this 
appear  in  the  literature,  but  the  Wallis  (1969)  model  appears  to  be  a  reasonable  limit  for  most 
data.  The  minimum  value  of  this  parameter  would  be  unity,  representing  a  smooth  gas— liquid 
interface.  Figure  1.6  reflects  this  range  of  uncertainty  in  this  key  parameter. 

Once  values  of  the  coordinates  Y  and  X  are  determined  from  the  equations  for  the 
transition,  they  can  be  readily  converted  to  dimensional  form  using  the  definitions  of  these 
parameters  in  Equations  1—3  and  1—2,  respectively.  Figures  1.8  and  1.9  also  show  the 
slug— to-annular  transition  on  flow  regime  maps  in  (jj,  jg)  coordinates.  The  set  of  plots  in 
Figure  1.8  shows  that  this  transition  is  approximately  linear  in  a  horizontal  pipe  at  earth 
gravity,  but  the  lower  portion  of  the  transition  approaches  a  constant  liquid  velocity  for 
downward  inclinations  and  approaches  a  constant  gas  velocity  for  upward  inclinations.  Bamea 
(1986)  shows  that  this  is  consistent  with  experimental  flow  regime  data  at  earth  gravity. 

Figure  1 .9  shows  that  the  transition  is  approximately  the  same  under  microgravity  conditions 
as  it  is  for  a  horizontal  pipe  at  earth  gravity.  This  is  readily  apparent  because  the  parameter  Y 
approaches  zero  in  the  case  of  either  low  inclinations  (0  =  0)  or  low  acceleration  level  (a  =  0). 
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Alternate  Regime  Map  Coordinates.  Simple  transformations  can  convert  the  flow  regime 
maps  in  (jlt  jg)  coordinates  as  described  above  to  mass  flux  versus  flow  quality  (G,  x) 
coordinates: 


G  =  (Pgig  +  Piii) 


(1-26) 


PgJg 

x  =  - 

(Pg, jg  +  Piii) 


PgJg 

.  G  . 


(1-27) 


where 

G  is  the  mass  flux  (kg/m2— s),  and 
x  is  the  quality  (ratio  of  vapor  mass  flow  rate  to  the 
total  mass  flow  rate). 

Figures  1.11  and  1.12  show  how  the  regime  maps  in  Figures  1.8  and  1.9  appear  when 
converted  to  these  coordinates.  It  is  particularly  interesting  to  note  that  under  microgravity 
conditions  (Figure  1.12c),  the  transitions  from  slug— to— bubbly  and  slug— to— annular  flow 
regimes  are  are  seen  to  be  lines  at  approximately  constant  flow  quality. 


1.1.4  Validation  with  Microgravity  Flow  Regime  Data 

As  Bamea  (1986)  points  out,  the  methods  and  models  described  here  apply  well  at  earth 
gravity  to  the  full  range  of  pipe  inclinations  from  vertical  upward,  to  horizontal,  to  vertical 
downward.  Examples  of  the  flow  regime  comparisons  arc  given  in  Bamea  (1986),  Bamea  et 
al.  (1982)  and  Taitel  (1980). 

The  purpose  of  this  section  is  to  illustrate  that  the  same  models  can  be  extrapolated  to 
microgravity  conditions.  At  this  time,  there  arc  four  key  sets  of  experimental  flow  regime  data 
obtained  under  microgravity  conditions.  These  include: 

•  Heppner  (1975,  1978), 

•  Lee  (1987)  (also  see  Kachnick  et  al.,  1987;, 

•  Dukler  et  al.  (1988),  and 

•  Chen  et  al.  (1988). 

All  represent  experiments  on  aircraft,  where  acceleration  levels  of  1 0- 2g  can  be  achieved  for 
periods  of  approximately  30  seconds  as  the  aircraft  negotiates  a  parabolic  trajectory.  Dukler  et 
al.  have  also  included  some  experimental  data  from  drop  tower  experiments  where  the  low 
acceleration  levels  are  for  periods  of  a  few  seconds.  No  flow  regime  data  in  pipes  are 
available  for  long-term  microgravity  conditions. 

As  discussed  by  Dukler,  the  Heppner  data  provided  the  initial  qualitative  concepts  for 
two-phase  flow  regimes  in  microgravity.  Unfortunately,  the  reconciliation  of  the  reported 
flow  regimes  with  recent  observations  of  the  flow  regimes  depends  upon  information  no  longer 
available  about  those  experiments.  For  that  reason,  the  validation  of  the  models  here  includes 
data  from  the  other  three  experiments  only. 


21 


Note:  Example  of  map  for  specific  conditions  -  not  for  general  use 
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Figure  1.11.  EFFECT  OF  PIPE  INCLINATION  ON  FLOW  REGIME  MAP  AT  EARTH  GRAVITY 
(MASS  FLUX  AND  QUALITY  COORDINATES) 
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Note:  Example  of  map  for  specific  conditions  -  not  for  general  use 


Figure  1.12.  EFFECT  OF  GRAVITY  ON  FLOW  REGIME  MAP  (.MASS  FLUX  AND  QUALITY 
'  COORDINATES) 


Table  1.4  summarizes  the  conditions  for  each  of  the  experiments  compared  below. 
Figures  1.13,  1.14,  and  1.15  compare  the  measured  flow  regime  data  from  these  experiments 
on  flow  regime  maps  in  velocity  coordinates.  (Since  the  dimensionless  groups  are  different  for 
each  flow  regime  transition,  it  is  not  possible  to  plot  data  from  different  regimes  on  a  single 
dimensionless  plot.) 

The  comparisons  with  the  Dukler  et  al.  flow  regime  data  (Figure  1.13)  show  quite  good 
agreement  with  the  mechanistic  analysis  methods.  It  should  be  noted  that  reducing  the 
acceleration  level  in  the  mechanistic  models  proposed  by  Bamea  (1986)  leads  to  effectively 
the  same  result  as  the  microgravity  modelling  methods  proposed  by  Dukler  et  al.  (1988). 
Although  the  starting  point  in  each  approach  is  different,  the  end  result  under  low  gravity 
conditions  is  not  significantly  different. 

The  comparisons  with  the  Lee  et  al.  flow  regime  data  (Figure  1.14)  also  agree  well  with 
the  mechanistic  analysis.  This  is  not  unexpected  since  the  test  conditions  are  quite  close  to 
those  of  the  Dukler  et  al.  experiments  (Table  1.4). 

The  limited  data  from  the  experiments  by  Chen  (at  one  mass  flux  with  varying  quality) 
show  the  location  of  the  transition  from  slug  to  annular  flow  regimes  (Figure  1.15). 

Comparison  with  the  locations  of  the  predicted  and  measured  transitions  for  the  previous  data 
(e.g.  Figure  1.13)  shows  that  the  transition  occurs  at  a  liquid— to— gas  velocity  ratio  about  5 
times  greater  in  the  Chen  experiment.  The  reason  is  that  the  density  ratio  (liquid  to  gas)  is 
about  28  times  larger  in  the  Chen  experiments  (with  Freon  R— 1 14  at  elevated  pressure)  than 
for  the  other  experiments  (with  air  and  water  at  atmospheric  pressure).  By  analysis  (Equations 
1—23  and  1-2)  it  can  be  shown  that  the  transition  depends  upon  the  square  root  of  the  density 
ratio.  Thus,  the  limited  experimental  data  confirm  the  scaling  of  the  low  gravity  flow  regime 
results  with  gas  density. 

With  these  low  gravity  comparisons  and  the  previous  comparisons  of  the  mechanistic 
model  with  data  at  earth  gravity,  it  is  seen  that  the  proposed  design  method  is  valid  for  a  wide 
range  of  inclinations  and  gravity  (acceleration)  levels. 


1.2  Pressure  Drop 

Once  the  flow  regime  is  selected  from  the  four  choices  according  to  the  methods 
described  in  Section  1.1,  or  if  the  flow  regime  to  be  analyzed  is  known  in  advance,  the 
pressure  drop  can  then  be  calculated.  This  section  describes  the  methods  used  to  calculate  the 
pressure  drop  in  the  four  regimes  identified: 

•  Annular, 

•  Stratified. 

•  Slug,  and 

•  Bubbly. 

For  space  systems  operating  under  microgravity  conditions,  the  annular  flow  regime  is  likely 
to  be  the  most  common.  This  is  because  the  annular  flow  regime  occurs  at  flow  qualities 
above  x  =  0.05  to  0.1.  Flow  qualities  approaching  1.0  take  maximum  advantage  of  the  latent 
heat  of  vaporization  of  the  liquid,  which  is  of  course  the  purpose  in  using  two-phase  systems 
in  the  first  place.  The  slug  and  bubbly  flow  regimes  will  be  less  common  during  steady-state 
operation  because  they  occur  at  low  flow  flow  qualities.  Slug  and  bubbly  flows  may  be 
encountered  at  turndown  conditions,  during  startup  and  shutdown,  or  in  short  regions  of 
high— shear  condensers  and  evaporators. 
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Table  1.4.  CONDITIONS  FOR  MICROGRAVITY  FLOW  REGIME  EXPERIMENTS 


EXPERIMENTER 

DUELER  ET  AL.  (1988) 

LEE  (1987) 

i  CHEN  ET  .AL  (  1988) 

1 

MODE 

LEAR  JET 

DROP  TOWER 

KC-135 

i  KC-135 

ACCELERATION  LEVEL 

=  10-2g 

s  l(Hg 

=  l<Hg 

=  NHg 
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AIR 

WATER 
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WATER 
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Figure  1.13.  FLOW  REGIME  COMPARISONS  FOR  DUKLER  ET  AL.  (1988) 
MICROGRAVITY  DATA 
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Section  1.2.1  provides  a  brief  introduction  to  the  recommended  models  and  the 
approaches  for  pressure  drop  calculations  in  pipes  of  constant,  circular  cross-section. 

Section  1.2.2  summarizes  the  information  which  the  designer  must  have  in  hand  in  order 
to  calculate  the  pressure  drop.  It  also  presents  dimensionless  parameters  and  graphs  which  can 
be  used  to  quickly  calculate  void  fraction  and  pressure  gradient  for  a  specific  regime.  This 
section  is  set  up  so  that  the  user  does  not  need  to  know  the  details  of  the  theory  or  the 
equations.  The  designer  approaches  the  calculations  with  knowledge  of  the  flow  regime,  the 
phase  flow  rates,  physical  geometry,  and  fluid  conditions.  Simple  dimensionless  parameters 
are  used  to  determine  a  void  fraction  or  a  two— phase  multiplier  on  plots.  The  pressure 
gradient  is  then  calculated  by  using  this  information  in  explicit  equations. 

Section  1.2.3  presents  the  basic  equations  for  the  pressure  gradient  calculations  in  each 
of  the  four  flow  regimes.  The  equations  are  presented  in  sufficient  detail  that  the  designer 
could  recreate  the  dimensionless  design  plots  presented  here  or  incorporate  the  calculations 
into  a  computer  routine. 

Finally,  Section  1.2.4  compares  the  pressure  gradient  calculations  for  the  recommended 
design  methods  against  limited  data  available  from  microgravity  experiments. 


1.2.1  Introduction 

Summary  of  Recommended  Models.  Table  1.5  summarizes  the  main  features  of  the 
recommended  Design  Manual  models  for  pressure  gradient.  The  model  used  for  each  of  the 
flow  regimes  is  mechanistic.  In  contrast  to  empirical  correlations,  which  abound  in  the 
two^thase  flow  literature,  the  mechanistic  approach  considers  the  governing  physics  in  the 
momentum  balances,  therefore  the  mechanistic  models  should  scale  well  over  a  range  of 
acceleration  magnitudes  (g-levels)  and  vectors. 

These  are  not  the  only  mechanistic  modelling  approaches  in  the  literature.  Generally, 
alternate  mechanistic  models  can  be  derived  at  various  levels  of  detail  for  each  of  the  four 
flow  regimes.  However,  the  models  recommended  here  are  sufficient  for  design  purposes. 

The  recommended  models  have  been  extensively  compared  against  experimental  data  at 
earth  gravity  in  various  flow  orientations  ranging  from  vertical  upward  to  vertical  downward. 
In  the  absence  of  gravity,  the  hydrostatic  component  of  the  pressure  gradient  becomes 
negligible  and  the  frictional  component  dominates,  so  the  models  reduce  to  the  appropriate 
limits  under  microgravitv  conditions. 

Bnef  Review  of  Modelling  Approaches  for  Two-Phase  Flow.  Table  1.5  shows  that  the 
mechanistic  models  for  the  slug  and  the  bubbly  flow  regimes  are  based  upon  a  "drift— flux" 
approach  and  that  the  models  for  the  annular  and  stratified  flow  regimes  are  based  upon  a 
"separated"  approach.  For  the  designer  unfamiliar  with  the  modelling  of  two-phase  flows,  the 
following  paragraphs  briefly  review  the  basic  modelling  approaches  which  include: 

•  Homogeneous, 

•  Drift-flux,  and 

•  Separated  (or  two— fluid). 
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Table  1.5.  RECOMMENDED  MODELS  FOR  PRESSURE  DROP  CALCULATIONS 


FLOW 

REGIME 

RECOMMENDED 

MODEL 

APPROACH 

KEY  PARAMETER 
FOR 

UNCERTAINTY 

BASELINE 
VALUE  OF 

KEY  PARAMETER 

RANGE  OF 
UNCERTAINTY 

FOR  KEY 
PARAMETER 

ANNULAR 

HEWITT  (1982) 
WALLIS  (1969) 

SEPARATED 

FLUIDS 

INTERFACIAL 
FRICTION 
(fi  /f-g) 

1  +  75(l-a) 

1  to 

rl  +  75(l-a)] 

|  STRATIFIED 

1 

i 

TAITEL-DUKLER 
( 1976) 

SEPARATED 

FLUIDS 

INTERFACIAL 
FRICTION 
(fi  /f„g) 

10 

1  to  10 

SLUG 

WALLIS  (1969) 

DRIFT-FLUX 

DISTRIBUTION 

PARAMETER 

(C0) 

,.3 

1.0  to  1.3 

BUBBLY 

WALLIS  (1969) 

DRIFT-FLUX 

DISTRIBUTION 

PARAMETER 

(C0) 

1-2 

1.0  to  1.2 

For  more  detailed  information,  the  reference  texts  on  two— phase  flow  modelling  by  Wallis 
(1969)  and  Collier  (1986)  are  recommended  reading. 

In  the  homogeneous  modelling  approach,  both  the  liquid  and  vapor  phases  are  assumed 
to  have  equal  velocities.  The  two-phase  flow  is  treated  as  if  it  is  a  single-phase  flow,  with 
mass,  momentum  and  energy  conservation  equations  using  effective  (or  average)  properties. 
The  techniques  of  solution  draw  upon  standard  procedures,  and  the  formulations  resemble 
familiar  ones  from  gas  dynamics  for  example.  The  key  is  to  define  suitable  average  properties 
for  the  two— phase  mixture.  Applying  this  method  is  most  appropriate  when  the  flow  is 
relatively  uniform  and  well-mixed.  Examples  include  emulsions  or  a  dispersion  of  small 
bubbles.  This  approach  is  also  useful  for  providing  quick  calculations  for  flow  conditions 
which  remain  to  be  analyzed  in  detail.  This  approach  is  the  basis  for  several  large  computer 
programs  for  steady-state  and  transient  thermal/hydraulics  (e.g.  RELAJP4/MOD5  and 
SINDA/FLUIND. 

The  gas  and  liquid  phases  need  not  necessarily  move  at  the  same  velocity.  The 
drift— flux  model  is  a  way  of  modifying  the  homogeneous  approach  to  account  for  relative 
motion  between  the  phases.  It  is  less  complicated  than  the  separated  flow  model  i.see  below > 
and  can  often  be  deduced  from  the  separated  flow  model  by  deriving  the  relative  motion  from 
a  force  balance  or  similar  simplifying  assumption.  It  has  found  most  practical  use  where  the 
relative  motion  is  due  to  a  body  force  field,  such  as  gravity.  That  is  why  it  is  used  for  the  slug 
and  bubbly  flow  regimes  here.  Equations  describing  the  relative  motions  between  the  phases 
in  slug  and  bubbly  flow  regimes  have  been  weil  studied.  Thermal/hvdrauiic  computer  codes 
such  as  RELAP4/MOD7,  RETRAN,  TRAC-PD2  and  COBRA/TRAC  use  this  approach  in 
one-dimensional  components. 

In  the  separated  (or  two— fluid)  approach  the  phases  are  treated  as  interwoven  continua, 
flowing  side  by  side  and  interacting  with  each  other.  Separate  conservation  laws  (mass, 
momentum,  and  energy)  are  written  for  each  phase,  thus  doubling  the  number  of  equations 
compared  with  the  homogeneous  approach.  The  variables  describing  the  morion  of  each 
phase,  such  as  velocity  and  pressure,  are  averaged  over  some  volume  or  surface  larger  than  the 
size  of  the  individual  entities  (e.g.  bubbles)  but  smaller  than  the  size  of  the  physical  system. 
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The  two-fluid  model  has  the  virtue  of  providing  a  self-consistent  framework  for  conservation 
laws  that,  because  of  greater  flexibility,  can  represent  many  more  phenomena  than 
homogeneous  or  drift— flux  models.  However,  knowledge  of  the  interaction  terms  in  the 
equations  (which  are  often  empirical),  closure  laws  (if  they  exist)  and  conditions  for  numerical 
convergence  are  still  at  an  elementary  level.  Fortunately,  empirical  interaction  terms  have 
been  derived  which  work  well  for  the  annular  and  stratified  flow  regimes  for  flow  in  pipes. 
Reasonable  limits  to  the  interaction  terms  can  be  defined  in  order  to  define  the  range  of  the 
uncertainty  in  calculating  the  pressure  gradient  and  void  fraction.  The  thermal/hydraulic  codes 
RELAP5,  ATHENA  and  TRAC— PF1  are  examples  of  computer  programs  which  use  this 
modelling  approach. 


1.2.2  Dimensionless  Design  Calculations 

This  section  illustrates  how  quick  estimates  of  the  void  fraction  and  the  pressure  gradient 
can  be  performed.  In  addition  to  the  flow  regime,  the  information  required  in  order  to 
calculate  these  quantities  includes: 


Fluid  Properties 

p  i  Liquid  phase  density 

pg  Gas  phase  density 

Pl  Liquid  phase  viscosity 

pg  Gas  phase  viscosity 

a  Surface  tension 


Geometry 

D  Pipe  diameter 


Acceleration 

a  Magnitude  of  acceleration 

9  Direction  (angle  between  acceleration  vector  and  the 
flow  direction) 


Flow  Rates 

ji  Liquid  phase  volumetric  flux  (Q/A) 
jg  Gas  phase  volumetric  flux  (Qg/A) 
j  The  total  volumetric  flux  (jg  +  ji),  and 

G  The  total  mass  flux  (pgjg  +  pji) 

See  the  section  below  for  the  annular,  stratified,  slug,  or  bubbly  flow  regime  as  appropriate. 

Annular  Flow  Regime.  The  procedure  to  calculate  the  pressure  gradient  in  the  annular 
flow  regime  is: 
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•  Calculate  the  dimensionless  parameters  (in  Equations  1-28,  1-29,  and 
1-30), 

•  Obtain  the  void  fraction  (see  Figure  1.16  or  1.18), 

•  Obtain  the  two-phase  multiplier  for  the  frictional  portion  of  the  pressure 
gradient  (see  Figure  1.17  or  1.19),  and 

•  Calculate  the  total  pressure  gradient  (Equation  1-32). 

The  dimensionless  parameters  for  the  annular  flow  regime  are: 


X  = 


Y  = 


m 

fpigjjl 

hi  J 

-  p4*]- 

f4Cg] 

_pr 

PgQig 

l  pg  J 

a(pi  -  pg)cos0 


’4C„ 

TT 


Pgtyj 


P*  =  (P/Pg) 


(1-28) 


(1-29) 


(1-30) 


where 


X 

is 

the 

Y 

is 

the 

P* 

is 

the 

a 

is 

the 

Pi 

is 

the 

Pg 

is 

the 

Ji 

je 

is 

the 

is 

the 

D 

is 

the 

0 

is 

the 

acceleration  vector  (°)  (e.g.  0°  is  colinear,  90°  is 
perpendicular,  and  180°  is  opposing,  corresponding 
to  downflow,  horizontal,  and  upflow  on  the  ground). 


Table  1.2  indicates  how  to  determine  the  values  of  the  friction  factor  coefficients  Cj,  Cg,  n  and 
m  to  use  in  Equations  1—28  and  1-29. 


The  dimensionless  parameters  in  Equations  1-28  and  1-29  are  the  same  as  those 
calculated  for  the  slug-to— annular  flow  regime,  Equations  1—2  and  1-3  (in  Section  1.1.2). 
Therefore  if  these  parameters  have  been  previously  calculated  as  part  of  the  flow  regime 
determination,  then  they  do  not  need  to  be  recalculated  here. 
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Figure  1.16  shows  the  liquid  fraction  (1-a)  for  the  annular  flow  regime.  The  void 
fraction  (one  minus  the  liquid  fraction)  can  be  estimated  using  only  the  dimensionless 
parameters  X  and  Y.  Note  that  the  line  represented  by  Y  =  0  would  apply  to  microgravity 
conditions.  (Generally,  the  absolute  value  of  Y  will  be  within  the  range  ±10  for  a  <  10“2g.) 


A  two— phase  multiplier  can  be  determined  for  the  annular  flow  regime  by  using  Figure 
1.17.  This  two— phase  multiplier  is  defined  as  the  ratio  of  the  frictional  pressure  gradient  in  the 
annular  flow  regime  to  the  equivalent  pressure  gradient  as  if  single-phase  (vapor)  flowed  in 
the  pipe.  That  is. 


‘  (dP/dz)f 
f  f  ELsi-s^l 


(1-31) 


where:  (dP/dz)f 

f 

lwg 


is  the  frictional  pressure  gradient  in  two— phase  flow, 
is  the  single-phase  friction  factor  for  the  gas 
phase  alone, 


and  the  other  parameters  are  defined  above. 


Once  the  frictional  two— phase  multiplier  is  determined,  the  total  pressure  gradient  can  be 
calculated  by: 


(dP/dz)  t  = 


f4f  1 

H  1  wg 

Pg  Jg2 

D 

2 

J 

+  a[(l  -  a)pi  +  apg]cos0 


(1-32) 


This  is  the  basic  pressure  gradient  equation  for  a  separated  (two— fluid)  flow.  In  Equation 
1—32,  the  first  term  on  the  right  hand  side  is  the  frictional  component,  and  the  second  term  is 
the  hydrostatic  (body  force)  component.  The  friction  factor  in  this  equation  can  be  evaluated 
from  a  smooth— pipe  model  such  as  the  Blasius  equation: 


f  16/R e„  for  laminar  flow,  Re2  <  1500 

fw  =  (1-33) 

0.046/Re„0-2  for  turbulent  flow.  Re.,  >  1500 

o  o 

or  any  other  appropriate  model.  For  example,  the  Colebrook  friction  factor  equation  found  in 
standard  fluid  flow  textbooks  could  be  used  if  the  value  of  the  pipe  roughness  would  have  a 
significant  effect  on  the  friction  factor.  Note  that  if  the  g— level  is  near  zero  (microgravitv )  or 
the  vector  angle  is  90°  (horizontal  pipe),  then  the  hydrostatic  term  in  Equation  1-32  becomes 
negligible  and  the  pressure  gradient  is  due  only  to  the  frictional  term. 


The  equations  and  figures  above  give  a  best  estimate  for  the  void  fraction  and  pressure 
gradient  in  the  annular  flow  regime.  Using  Figures  1—18  and  1-19  in  place  of  Figures  1-10 
and  1—17  allows  the  uncertainty  to  be  assessed.  (The  uncertainty  specifically  lies  in  the 
interfacial  friction  between  the  gas  and  liquid  phases,  as  described  in  Section  1.2.3.)  The 
variation  in  the  void  fraction  obtained  using  Figure  1—16  versus  1-18  and  in  pressure  gradient 
values  from  Figure  1—17  versus  1—19  represents  the  range  of  the  uncertainty. 
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Figure  1.16.  BEST-ESTIMATE  VOID  FRACTION  FOR  THE  ANNULAR  FLOW  REGIME 

(WALLIS  MODEL,  f./f  =  LI  +  75(l-<x)J) 
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Figure  1.17.  BEST-ESTIMATE  TWO-PHASE  MULTIPLIER  FOR  THE  ANNULAR  FLOW  REGIME 
(WALLIS  MODEL  f./f  =  11  +  75(l-a)J) 
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Figure  1.18.  LIMITING  VOID  FRACTION  FOR  THE  ANNULAR  FLOW  REGIME 

(SMOOTH  INTERFACE;  f./f  =  1) 
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Figure  1.19.  LIMITING  TWO-PHASE  MULTIPLIER  FOR  THE  ANNULAR  FLOW  REGIME 
(SMOOTH  INTERFACE;  f;/f  =  1) 
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Stratified  Flow  Regime.  The  procedure  for  calculating  the  void  fraction  and  the  pressure 
gradient  in  the  stratified  flow  regime  is  the  same  as  for  the  annular  flow  regime,  only  the 
figures  used  to  estimate  the  void  fraction  and  the  two— phase  multiplier  are  different.  The 
procedure  to  calculate  the  pressure  gradient  in  the  stratified  flow  regime  is: 

•  Calculate  the  dimensionless  parameters  (in  Equations  1—28,  1—29,  and 


•  Obtain  the  void  fraction  (see  Figure  1.20  or  1.22), 

•  Obtain  the  frictional  two-phase  multiplier  (see  Figure  1.21  or  1.23),  and 

•  Calculate  the  total  pressure  gradient  (Equation  1-32). 

Figure  1.20  shows  that  the  void  fraction  a  in  the  stratified  flow  regime  can  be  estimated 
using  the  dimensionless  parameters  X  and  Y,  just  as  it  was  for  the  annular  flow  regime.  Note 
that  the  line  represented  by  Y  =  0  is  appropriate  for  microgravity  conditions. 

Figure  1.21  plots  the  two— phase  multiplier  for  the  frictional  portion  of  the  pressure 
gradient  in  the  stratified  flow  regime.  Once  the  two— phase  multiplier  is  determined,  Equation 
1—32  can  again  be  used  to  determine  the  ac,uul  pressure  gradient  in  the  stratified  flow  regime. 

Using  Figures  1.22  and  1.23  instead  of  Figures  1.20  and  1.21  allows  the  uncertainty  in 
modelling  the  void  fraction  and  the  pressure  gradient  to  be  assessed.  (The  uncertainty 
specifically  lies  in  the  interfacial  friction  between  the  gas  and  liquid  phases,  as  described  in 
Section  1.2.3.)  The  uncertainty  in  the  void  fraction  is  given  by  the  difference  in  values 
obtained  with  Figures  1.20  and  1.22.  The  uncertainty  in  the  pressure  gradient  is  given  by 
using  the  values  obtained  from  Figures  1.21  and  1.23. 

Slug  Flow  Regime.  The  procedure  to  calculate  the  pressure  gradient  in  the  slug  flow 
regime  is: 

•  Calculate  the  dimensionless  velocity  ratio  (Equation  1—34), 

•  Obtain  the  void  fraction  (see  Figure  1-24  or  1-25),  and 

•  Calculate  the  pressure  gradient  (Equation  1—35). 

The  dimensionless  velocity  parameter  required  for  the  slug  flow  regime  is  the  ratio: 

M  =  (j  j/jg;  0-34) 


Figure  1.24  shows  that  the  void  fraction  can  be  estimated  using  this  parameter.  The 
drift— flux  relationship  shown  here  is  simplified  —  as  discussed  in  Section  1.2.3  —  but  should 
provide  adequate  estimates,  especially  for  microgravity  conditions. 

Once  the  void  fraction  is  determined,  the  total  pressure  gradient  can  be  readily  evaluated 
from  the  simple  equation: 


(dP/dz)  (  =- 


2fw  lg  il2 


(1 -a)  +  a[(l  -  a)pi  +  apg]cos0 


(1-35) 
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Figure  1.20.  BEST-ESTIMATE  VOID  FRACTION  FOR  THE  STRATIFIED  FLOW  REGIME 
(f./f  =  10) 
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Figure  1.21.  BEST-ESTIMATE  TWO-PHASE  MULTIPLIER  FOR  THE  STRATIFIED  FLOW 

REGIME  (f./f  =  10) 
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?ure  1.22.  LIMITING  VOID  FRACTION  FOR  THE  STRATIFIED  FLOW  REGIME 
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Figure  1.23.  LIMITING  TWO-PHASE  MULTIPLIER  FOR  THE  STRATIFIED  FLOW  REGIME 
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Figure  1.24.  BEST-ESTIMATE  VOID  FRACTION  FOR  THE  SLUG  FLOW  REGIME  (C  =  1.3) 


The  first  term  on  the  right  hand  side  of  this  equation  is  the  frictional  component  of  the 
pressure  gradient,  and  the  second  term  is  the  hydrostatic  component.  Note  that  if  the 
magnitude  of  the  acceleration  becomes  small,  as  it  would  under  microgravity  conditions,  then 
only  the  frictional  component  remains. 

Using  Figure  1.25  instead  of  Figure  1.24  allows  the  uncertainty  in  modelling  the  void 
fraction  and  the  pressure  gradient  to  be  assessed.  (The  uncertainty  specifically  lies  in  the 
velocity  difference  between  the  gas  and  liquid  phases,  as  described  in  Section  1.2.3  )  Figure 
1.25  represents  a  iimiting  case  of  homogeneous  flow  (see  Section  1. 2.1).  The  variation  in  the 
void  fraction  and  pressure  gradient  values  obtained  from  using  Figure  1.25  in  place  of  Figure 
1.24  illustrates  the  range  of  the  uncertainty.  Void  fraction  values  will  be  higher  and  hence  the 
pressure  gradient  ;0'ver  using  the  homogeneous  approach.  The  magnitude  of  the  ertect  is 
about  30%  on  ihe  void  fraction  and  50%  on  both  pressure  gradient  components. 

Bubblv  Flow  Regime.  The  procedure  to  calculate  the  pressure  gradient  in  the  bubbly 
flow  regime  is: 

•  Calculate  the  dimensionless  velocity  ratio  (Equation  1—34), 

•  Obtain  the  void  fraction  (see  Figure  1.26  or  1.25),  and 

•  Calculate  the  pressure  gradient  (Equation  1-36), 

which  is  very  similar  to  the  procedure  in  the  slug  flow  regime. 
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Figure  1.25.  LIMITING  (HOMOGENEOUS)  VOID  FRACTION  FOR  THE  SLUG  OR  BUBBLY  FLOW 
REGIME  (C  =  1.0) 

Figure  1.26  shows  that  the  void  fraction  can  be  estimated  using  the  dimensionless 
velocity  ratio  M.  The  dimensionless  velocity  ratio  required  for  the  bubbly  flow  regime  is  the 
same  as  for  the  slug  flow  regime  (Equation  1—34).  The  drift— flux  relationship  shown  in  this 
figure  is  simplified  -  as  discussed  in  Section  1.2.3  -  but  should  provide  adequate  estimates, 
especially  for  microgravity  conditions. 

Once  the  void  fraction  is  determined  the  total  pressure  gradient  can  be  calculated  from 

(dP/dz)t  =  -  |  +  a[(l  -  a)pi  +  apg]cos0  (1-36) 

As  before,  the  first  term  on  the  right  hand  side  of  this  equation  is  the  frictional  component  or 
the  pressure  gradient,  and  the  second  term  is  the  hydrostatic  component. 

Using  Figure  1.25  instead  of  Figure  1.26  allows  the  uncertainty  in  modelling  the  void 
fraction  and  the  pressure  gradient  to  be  assessed.  vTne  uncertainty  specifically  lies  in  the 
velocitv  difference  between  the  gas  and  liquid  phases,  as  described  in  Section  1.2.3.)  Figure 
1.25  repiesents  the  limiting  case  of  homogeneous  flow  (see  Section  1.2.1).  The  variation  in 
the  void  fraction  and  pressure  gradient  values  obtained  from  using  Figure  1.25  in  place  of 
Figure  1.26  illustrates  the  range  of  the  uncertainty.  The  magnitude  of  the  effect  is  about  20% 
in  the  void  fraction.  This  does  net  affect  the  frictional  component  of  the  pressure  gradient  at 
all,  but  does  affect  the  hydrostatic  component  by  20%. 
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Figure  1.26.  BEST-ESTIMATE  VOID  FRACTION  FOR  THE  BUBBLY  ROW  REGIME 

(Co=  1-2) 


1.2.3  Equations  for  Pressure  Drop  Calculations 

This  section  presents  the  exact  analytical  equations  used  to  calculate  the  void  fraction 
and  the  pressure  gradient  in  each  of  the  four  flow  regimes.  Additional  design  maps  with 
alternate  values  of  the  key  modelling  parameters  can  be  generated  if  desired  using  the 
equations  presented  in  this  section. 


Annular  Flow  Regime.  Derivation  of  the  basic  modei  for  the  annuiar  How  regime  is 
'discussed  m  the  references  by  Wallis  i  1969''  and  Hewitt  ( 19821.  in  the  separated  fluid 
approach  used  in  this  regime,  individual  momentum  equations  are  written  for  balances  on  me 
gas  and  liquid  phases.  These  two  momentum  equations  are  then  subtracted  (to  eliminate  the 
pressure  gradient)  and.  using  the  definitions  of  the  dimensionless  parameters  X  ana  V.  the 
following  equation 

f  (f/fwg)(l  -  a)2] 

X2  = - (1— a)3  Y  (1-371 

l  a2-5 

is  obtained  to  relate  the  void  fraction  a  and  the  dimensionless  parameters  X  and  Y. 


The  parameter  (fi/fwg)  in  Equation  1—37  represents  the  interfacial  friction  between  the 
gas  and  the  liquid  phases.  For  the  purposes  of  this  Design  Manual,  the  empirical  relationship 
proposed  by  Wallis  (1969)* 

(fi/fwg)  =  [1 +  75(1 -a)]  (1-38) 

is  recommended.  This  relationship  is  also  used  by  Dukler  (1988)  in  his  analysis  of  flow 
regimes  in  microgravity.  Figure  1.16,  the  void  fraction  map  for  the  annular  flow  regime,  has 
been  generated  using  Equation  1-38  in  Equation  1-37.  Various  models  for  this  relationship 
appear  in  the  literature  for  the  annular  flow  regime,  but  the  Wallis  (1969)  model  appears  to  be 
a  reasonable  limit  for  most  data. 

The  interfacial  friction  is  the  key  modelling  parameter  in  the  annular  flow  regime.  The 
major  uncertainties  in  the  modelling  are  represented  by  it.  The  minimum  value  of  this 
parameter  would  be  unity  (fj/fwg  =  1),  representing  a  smooth  interface  between  the  gas  and  the 
liquid  phases.  Figure  1.18  is  based  on  this  limiting  value  <Y  the  key  parameter.  Minimum 
values  of  the  void  fraction  are  obtained  using  this  limit. 

The  "separated  model"  equation  for  the  total  pressure  gradient  in  the  annular  flow  regime 
is: 


(dP/dz;  t  =  - 


2f  w  lPij  i2 
D(  1  -  a)2 


+  a[(l -a)Pi  -  apg]cos0 


(1-39) 


which  is  obtained  by  adding  the  momentum  equations  for  each  phase  in  order  to  eliminate  the 
interfacial  friction  factor  ratio.  The  first  term  on  the  right  hand  side  of  this  equation  represents 
the  frictional  portion  of  the  pressure  gradient.  This  equation  could  be  evaluated  directly  from 
known  information  once  the  void  fraction  is  calculated. 


For  convenience  in  this  Manual,  a  two-phase  multiplier  for  the  frictional  portion  of  the 
pressure  gradient  is  defined  as 


4>go2  -  (i-_  a-)2  (1_40) 

Since  the  void  fraction  is  related  to  the  parameters  X  and  Y  by  Equation  1—37,  the  two— phase 
multiplier  can  be  plotted  as  a  function  of  the  same  parameters.  Figures  1.17  and  1.19  have 
been  prepared  using  the  Wal'is  model  1 1  +  75(1  -  a)]  and  the  limiting  value  of  1.0  for  the 
interracial  friction,  respectively.  The  pressure  gradient  is  then  evaluated  using  Equation  1-52. 

More  complicated  models  could  be  written  for  the  annular  flow  regime,  involving  liquid 
droplet  entrainment  and  deposition.  For  most  cases  of  interest  to  thermal  management 
systems,  liquid  dropiet  entrainment  should  not  be  significant.  The  effect  of  liquid  droplet 
entrainment  on  the  pressure  gradient  is  relatively  small  in  any  event. 


Stratified  Flow  Regime.  The  approach  to  the  model  for  the  stratified  flow  regime  is 
similar  to  that  used  for  the  annular  regime.  Taitel  and  Dukler  (1976)  first  developed  this 
approach  for  the  stratified  flow  reg'rne,  although  the  form  used  was  first  proposed  by  Lockhart 
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and  Martinelli  in  the  1950s.  Separate  momentum  equations  are  written  for  the  gas  and  liquid 
phases  as  before.  The  main  difference  compared  with  the  annular  flow  model  is  in  the 
geometric  parameters  which  describe  the  interfaces  between  the  wall  and  the  fluids  and 
between  the  liquid  and  vapor  phases.  The  void  fraction  relationship  for  the  stratified  flow 
regime  is: 


X2  = 


[Ui*Dj*]  ~"(U  i*)2(S  j  */A !*) 


4Y 


(1-41) 


Where  the  dimensionless  parameters  X  and  Y  have  the  usual  definitions.  The  "starred" 
parameters  in  Equation  1—41  are  the  geometric  parameters  for  the  stratified  flow  regime. 

Table  1.3  summarizes  these  parameters  as  a  function  of  the  dimensionless  liquid  level  and 
shows  how  they  are  related  to  the  void  fraction. 

The  values  of  the  parameters  n  and  m  are  determined  from  the  Reynolds  numbers  for  the 
liquid  and  gas  phases.  Table  1.2  provides  the  equations  necessary  to  determine  these  values. 

Just  as  in  the  annular  flow  regime,  the  parameter  (f;/fwg)  represents  the  interfacial 
friction,  the  key  modelling  parameter.  Various  empirical  correlations  have  been  proposed  for 
this  relationship  in  the  stratified  flow  regime.  In  this  Design  Manual  a  constant  value  of 
(f;/fwg)  =  10  is  recommended.  This  selection  is  based  upon  extensive  work  by  the  authors 
reviewing  alternative  modelling  approaches  and  specific  models  representing  the  interfacial 
shear  in  stratified  flows  at  high  gas  densities. 

The  void  fraction  relationship  in  Figure  1.20  represents  Equation  1-41  using  the 
recommended  value  for  the  interfacial  friction.  Values  of  n  =  m  =  0.2  have  also  been  used 
(assuming  turbulent  flow  in  both  phases).  Other  void  fraction  relationships  for  the  stratified 
flow  regime  could  be  derived  using  alternative  models  for  the  interfacial  friction.  Figure  1.22 
is  for  the  limiting  case  of  a  smooth  gas— liquid  interface.  This  limit  gives  the  minimum  void 
fraction.  It  is  also  the  limit  used  by  Taitel  and  Dukler  in  the  original  derivation  of  this  model; 
the  model  has  been  modified  here  to  accommodate  alternative  interfacial  friction  models. 

The  dimensional  equation  for  the  total  pressure  gradient  in  the  stratified  How  regime  is: 


(dP/dz)  t 


2fwiP  1  j  r  '  7  “^WKpgjg"!  7 

-  (U,*)“(S,*)  -  i — - - (U„*)  (S.,*)  U-42) 

l  4D  j  (_  ttD  j 


+  a[(l  -  a)p]  +  apgjcosG 


which  is  obtained  by  adding  the  momentum  equations  for  each  phase  in  order  to  eliminate  the 
interfacial  shear  ratio.  The  first  two  terms  on  the  right  hand  side  of  this  equation  represent  the 
frictional  portion  of  the  pressure  gradient.  This  equation  could  be  evaluated  directly  from 
known  information  once  the  void  fraction  is  calculated. 
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For  convenience,  a  two— phase  multiplier  for  the  frictional  portion  of  the  pressure 
gradient  can  be  defined.  For  the  stratified  flow  regime,  this  multiplier  has  the  form: 

<t>go2  =  [i]  [X^(U  [*)2(S  [*)  +  (Ug*)2(Sg*)]  (1-43) 

Since  the  dimensionless  geometric  parameters  Ups  Ug*,  Si*,  and  S„*  are  related  to  the  void 
fraction,  which  in  turn  is  a  function  of  the  parameters  X  and  Y  by  Equation  1—41,  this  allows 
the  two-phase  multiplier  to  be  plotted  as  a  function  of  the  X  and  Y  parameters,  just  as  for  the 
annular  flow  regime.  Figures  1.21  and  1.23  have  been  prepared  using  the  recommended 
interfacial  shear  model  (fi/fwg)  =  10  and  the  limiting  value  of  unity  for  the  interfacial  friction, 
respectively. 

Defining  the  two-phase  multipliers  for  the  annular  and  the  stratified  flow  regimes  in  the 
manner  shown  allows  a  comparison  of  the  .void  fractions  and  the  pressure  drops  in  the  two 
regimes.  Comparison  of  Figures  1.16  and  1.20  shows  that  the  void  fraction  tends  to  be 
somewhat  greater  in  the  stratified  flow  regime.  Comparison  of  Figures  1.17  and  1.21  shows 
that  the  two— phase  multiplier,  and  hence  the  frictional  pressure  gradient,  tends  to  be 
significantly  greater  in  the  annular  flow  regime.  (This  is  due  to  the  increased  interfacial  area 
in  the  annular  flow  regime.) 

More  complicated  models  could  be  written  for  the  stratified  flow  regime,  involving 
liquid  droplet  entrainment  and  deposition.  For  most  cases  of  interest  to  thermal  management 
systems,  liquid  droplet  entrainment  would  not  be  significant.  The  effect  of  liquid  droplet 
entrainment  on  the  pressure  gradient  is  relatively  small  in  any  event. 


Slug  Flow  Regime.  The  model  for  void  fraction  in  the  slug  flow  regime  is  based  upon  a 
drift— flux  approach,  which  requires  an  empirical  relationship  between  the  void  fraction  and  the 
flow  rates  of  each  phase.  The  basic  modelling  approach  for  the  pressure  gradient  is  a 
homogeneous  one,  as  suggested  by  Wallis  (1969). 

The  drift- -flux  equation  recommended  for  the  slug  flow  regime  has  the  form: 
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where 


C0  is  a  constant  typically  near  1.3,  and 
k  is  a  constant  whose  value  is  given  in  Table  1.6,  and 

the  other  parameters  have  been  defined  previously.  The  values  in  Table  1.6  have  been 
selected  based  upon  our  review  of  almost  two  dozen  experiments. 

Equation  1 — 44  can  be  rearranged  in  the  form 
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Using  the  definitions  of  the  dimensionless  parameters  M,  p*,  and  Fg  (Equations  1—6,  1—7,  and 
1—1),  Equation  1-45  can  be  written  in  dimensionless  form  as 

1  i  k  r  no.5  . 

M  = - 1 -  (1-46) 

C0ct  J  Fg  Lp*_ 

or  rearranging  the  equation  to  solve  for  the  void  fraction 

1 

a  = - (1-47) 

C0(l  +  M)  +  (k/Fgp*o.5) 

Figure  1.24  plots  the  void  fraction  relationship  in  Equation  1-47  using  the  recommended 
value  C0  =  1.3  with  the  approximation  that  k  =  0  under  all  conditions.  The  contribution  of  the 
term  with  the  constant  k  in  Equation  1—47  is  generally  small,  so  the  approximation  with  k  =  0 
is  valid.  The  approximation  should  be  very  good  for  microgravity  conditions  because  this  term 
also  becomes  negligible  when  the  magnitude  of  the  acceleration  approaches  zero. 

The  key  uncertainty  in  the  drift— flux  relationship  of  Equation  1—47  is  the  value  of  the 
parameter  C0.  The  limiting  case  of  C0  =  1.0  with  k  =  0  represents  a  homogeneous  flow,  i.e. 
no  drift-flux.  This  is  the  relationship  plotted  in  Figure  1.25.  This  limit  maximizes  the  void 
fraction  in  the  slug  flow  regime. 

The  basic  homogeneous  form  of  the  equation  for  pressure  gradient  is: 

'2fw  l  pjr 

(dP/dz)t=-  -  +  a[(l  -  a)pi  +  apg]cos0  (1-48) 


where  pm  is  the  homogeneous  mixture  density. 


Table  1.6.  DRIFT-FLUX  PARAMETERS  FOR  THE  SLUG 
FLOW  REGIME 
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Since  the  gas  density  is  generally  much  smaller  than  the  liquid  density,  the  simple 
approximation 


pm  =  (1  -<x)p, 

can  be  used  to  obtain  a  slightly  simpler  form  of  the  equation,  which  is 

f2fwipg2] 


(dP/dz)t  =  - 


(1  -  a)  +  a[(l  -  a)pi  +  apg]cos0 


(1-49) 


(1-50) 


This  is  the  pressure  gradient  equation  solved  in  the  slug  flow  regime. 


More  complicated  models  could  be  written  for  the  slug  flow  regime  also.  These  models 
would  involve  asymmetric  gas  bubbles  or  corrections  for  long  gas  bubbles.  Models  with  as 
many  as  17  parameters  have  been  developed  (Fernandes,  1982).  For  microgravity  flow 
conditions,  these  detailed  models  are  believed  to  be  unnecessary. 


Bubbly  F^w  Regime.  The  drift— flux  approach  is  also  used  to  calculate  the  void  fraction 
in  the  bubbly  flow  regime.  Similarly  to  the  slug  flow  regime,  the  equation  for  the  pressure 
gradient  is  a  homogeneous  one,  as  suggested  by  Wallis  (1969).  The  empirical  drift-flux 
model  used  with  the  pressure  gradient  equation  is  based  on  extensive  research. 


The  drift— flux  equation  recommended  for  the  bubbly  flow  regime  has  the  form 
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(1-51) 


where 

CQ  is  a  constant  typically  near  1.2, 
k  is  a  constant  whose  value  is  given  in  Table  1.7, 

and  the  other  parameters  have  been  defined  previously.  The  values  of  Cc  and  k  are  based 
upon  work  by  Mishima  and  Ishii  (1984),  Harmathy  (1960),  Zuber  and  Hench  (1962),  and 
Martin  (1973)  for  swarms  of  bubbles  flowing  in  a  liquid. 


Equation  1—51  can  be  rearranged  to  obtain 
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Using  the  dimensionless  parameters  M,  p*,  and  Kg  defined  previously  (Equations  1-6,  1—7 
and  1-8),  Equation  1—52  can  be  written  in  the  dimensionless  form 


M  = 
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(1-53) 


and  rearranging  to  solve  for  the  void  fraction  gives 
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(1-54) 


C0(l  +  M)  +  (k/Kgp*0-5) 

Figure  1.26  plots  the  void  fraction  relationship  in  Equation  1-54  using  the  recommended  value 
C0  =  1.2  with  the  approximation  that  k  =  0  under  all  conditions.  For  the  high  phase  velocities 
necessary  to  cause  bubbly  flow  in  pipes,  the  contribution  of  the  term  with  the  constant  k  in 
Equation  1—54  is  small,  so  the  approximation  with  k  =  0  is  valid.  The  approximation  should 
be  very  good  for  microgravity  conditions  because  this  term  also  becomes  negligible  when  the 
magnitude  of  the  acceleration  is  very  small. 

The  key  uncertainty  in  the  drift— flux  relationship  of  Equation  1—54  is  the  value  of  the 
parameter  C0.  The  limiting  case  of  C0  =  1.0  (with  k  =  0)  represents  a  homogeneous  flow,  i.e. 
no  drift— flux.  This  is  the  relationship  plotted  in  Figure  1.25.  This  limit  maximizes  the  void 
fraction  in  the  bubbly  (or  slug)  flow  regime. 

Equation  1-48  is  the  homogeneous  equation  for  the  total  pressure  gradient.  For  the 
bubbly  flow  regime,  this  equation  can  be  modified  using  the  approximation  that 

Pnd  =  (Pgjg  +  Prii)  =  G  (1-55) 

in  order  to  obtain  the  pressure  gradient  equation  in  the  form 

[2fw  i  Gjl 

(dP/dz)t  =  -  -  +  a[(l  -  a)pj  +  ap„]cos0  (1-56) 

D  ] 

which  is  the  equation  solved  in  the  bubbly  flow  regime. 


I  Table  1.7.  DRIFT-FLUX  PARAMETERS  FOR  THE  BUBBLY  j 
FLOW  REGIME 
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More  complicated  models  could  be  written  for  the  bubbly  flow  regime,  involving 
determining  the  sizes  of  vapor  bubbles  and  calculating  the  interfacial  forces  on  gas  bubbles. 
For  bubbly  flow  under  most  conditions,  the  simple  drift— flux  model  is  adequate. 


1.2.4  Validation  with  Microgravity  Pressure  Drop  Data 

The  purpose  of  this  section  is  to  demonstrate  that  the  mechanistic  pressure  drop  models 
can  be  used  for  microgravity  conditions.  Data  obtained  under  microgravity  conditions  are 
presently  available  from  Chen  et  al.  (1989). 

Additional  experiments  are  planned  by  NASA  Lewis  Research  Center  (U.  of  Houston) 
and  the  Weapons  Laboratory  (Creare  and  Foster— Miller);  however,  no  data  have  yet  been 
produced  by  these  experiments.  All  represent  experiments  on  aircraft,  where  acceleration 
levels  of  10~2g  can  be  achieved  for  periods  of  approximately  30  seconds  as  the  aircraft 
negotiates  a  parabolic  trajectory.  No  pressure  drop  data  are  available  for  long-term 
microgravity  conditions. 

Table  1.4  summarizes  the  conditions  for  the  Chen  experiments  with  Freon  R-l  14. 
Experiments  were  performed  at  a  mass  flux  of  approximately  40  kg/m2— s  in  a  tube  15.8  mm  in 
iameter  and  1830  mm  long.  Figures  1.27  and  1.28  compare  the  measured  pressure  drop  data 
from  these  experiments  with  calculations  from  the  models.  The  measured  pressure  drops  are 
in  the  range  of  0.07  to  1.5  kPa  (0.01  to  0.22  psi).  Note  that  the  reported  accuracy  of  the 
transducer  used  to  make  the  pressure  drop  measurements  is  0.33  kPa  (0.05  psi)  so  that  some  of 
the  smaller  pressure  drops  reported  have  a  large  uncertainty. 

In  ground  tests,  the  two-phase  test  section  was  primarily  in  the  stratified  flow  regime. 
Figure  1.27  shows  that  the  calculated  pressure  drops  in  the  stratified  flow  regime  are  in  good 
agreement  with  measured  pressure  drops  at  the  larger,  more  accurate  values.  Using  the  model 
for  the  annular  flow  regime  would  overpredict  the  pressure  drop  by  a  factor  of  three  or  more 
for  the  stratified  data. 

Figure  1.28  shows  that  the  pressure  drop  data  in  the  annular  flow  regime  are  bounded  by 
the  models  suggested  in  this  Design  Manual.  Using  the  Wallis  model  for  interfacial  shear  (see 
Figures  1.16  and  1.17),  the  pressure  drop  is  overpredicted  by  about  50%.  Using  the  limiting 
case  of  a  smooth  interface  (see  Figures  1.17  and  1.18)  the  pressure  drop  is  underpredicted  by 
about  a  factor  of  two.  Chen  et  al.  suggested  a  modified  interfacial  shear  relationship  to  match 
their  experimental  data.  Figure  1.28  also  shows  that  this  model  succeeds  in  matching  the  data, 
though  the  difference  between  this  model  and  the  Wallis  model  is  within  the  accuracy  of  the 
pressure  measurements. 

With  these  low  gravity  comparisons  and  numerous  previous  comparisons  of  the 
mechanistic  models  with  data  at  earth  gravity,  the  proposed  design  method  appears  to  be  valid 
for  the  annular  flow  regime  in  microgravity.  Additional  microgravity  data  in  the  other  flow 
regimes  are  needed  in  order  to  validate  the  models.  The  experiments  in  progress  by  the 
University  of  Houston  and  the  Weapons  Laboratory  should  supply  the  needed  data. 
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Figure  1.27.  COMPARISON  BETWEEN  GROUND  TEST  DATA  (CHEN  ET  AL.)  AND  THE 
PRESSURE  DROP  MODEL  FOR  TOE  STRATIFIED  FLOW  REGIME 
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Figure  1.28.  COMPARISON  BETWEEN  MICROORAVTTY  TEST  DATA  (CHEN  ET  AL.)  AND  THE 
PRESSURE  DROP  MODEL  FOR  TOE  ANNULAR  FLOW  REGIME 
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2  EVAPORATORS:  Forced  Convection  Dominated  Tubes 


Evaporators  are  the  points  of  contact  between  two— phase  thermal  management  systems 
and  spacecraft  components  (computers,  the  environment,  payloads,  power  sources,  sensors, 
etc.)  which  produce  heat  and  therefore  require  cooling.  Heat  is  absorbed  from  spacecraft 
components  into  the  thermal  management  system  by  vaporizing  the  working  fluid,  hence  the 
name  "evaporator."  There  are  two  major  categories  of  evaporators  for  space  thermal 
management  systems:  those  which  produce  a  two— phase  mixture  of  vapor  and  liquid  and 
those  which  produce  only  single— phase  vapor.  Designs  for  forced  convection  evaporators 
include  straight  and  curved  flow  channels.  In  curved  tubes  ("swirl"  flow),  heat  transfer  can  be 
enhanced  by  centrifugal  forces  at  the  expense  of  larger  pressure  drops  than  straight  tubes. 
Devices  to  produce  single— phase  vapor  generally  use  capillary-controlled  evaporators  or 
evaporate  extremely  thin  liquid  films.  The  goal  of  microgravity  evaporator  design  is  to 
provide  heat  removal  in  a  small,  light-weight  component  which  operates  independently  of 
spacecraft  acceleration  and  which  does  not  require  excessive  pumping  power. 

This  section  provides  basic  relationships  between  coolant  properties,  evaporator 
geometry,  fluid  flow,  temperature,  heat  transfer  and  pressure  losses  for  the  evaporator  concept 
with  a  straight  tube  of  circular  cross-section  where  boiling  (vaporization)  is  dominated  by 
forced-convection.  This  topic  was  selected  for  this  Design  Manual  because  it  is  the  most 
basic  approach  and  is  therefore  the  starting  point  for  understanding  alternate  concepts.  Some 
of  the  more  exotic  two— phase  evaporator  concepts  (such  as  swirl  evaporators)  are  simple 
variations  on  the  straight  tube  concept.  Furthermore,  the  approach  described  here  is  used  in 
microgravity  evaporator  experiments,  and  data  will  soon  be  available  for  comparison  with  the 

d'-sign  mode1*!. 
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For  microgravity  applications  these  evaporators  should  be  designed  so  that  heat  transfer 
is  dominated  by  forced  convection.  Heat  transfer  in  the  evaporator  actually  has  two 
mechanisms:  forced  convection  and  nucleate  boiling.  Nucleate  boiling  heat  transfer  may 
depend  to  some  extent  on  the  spacecraft  acceleration,  however  this  dependence  is  still  poorly 
understood  and  no  satisfactory  test  data  are  yet  available.  Forced  convection,  on  the  other 
hand,  is  entirely  independent  of  acceleration.  Therefore  the  approach  recommended  in  this 
Design  Manual  is  to  design  the  evaporator  so  that  forced  convection  is  the  dominant 
mechanism.  An  evaporator  designed  this  way  operates  independently  of  the  spacecraft 
acceleration  regardless  of  possible  acceleration  effects  on  nucleate  boiling  heat  transfer. 

First,  design  maps  and  calculation  methods  for  heat  transfer  are  presented  (Section  2.1). 
This  section  instructs  the  designer  how  to  size  the  evaporator  tube  to  obtain  the  desired  outlet 
quality  and  to  operate  independently  of  spacecraft  acceleration.  Section  2.2  concerns  pressure 
losses  in  a  straight  tube  evaporator.  The  designer  uses  the  maps  and  methods  in  this  section  to 
calculate  the  pressure  drop  through  the  evaporator  in  order  to  estimate  system  pumping  power 
requirements. 

Some  of  the  design  maps  in  this  Section  are  less  general  than  in  the  previous  Section  on 
flow  regimes  and  pressure  losses.  This  is  because  heat  transfer  is  still  largely  an  empirical 
science.  The  present  state  of  the  art  does  not  provide  generalized  mechanistic  design  maps  for 
many  of  the  important  phenomena.  In  such  cases  the  designer  must  carry  out  more  detailed 
calculations  for  specific  design  cases. 
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2.1  Heat  Transfer 

2.1.1  Introduction 

This  section  describes  the  main  phenomena  for  two— phase  heat  transfer  in  an  evaporator 
tube.  In  general,  an  evaporator  has  several  possible  thermal  boundary  conditions:  If  the  inlet 
conditions  are  constant,  then  either  the  temperature  of  the  heated  surface  and/or  the  heat  flux 
on  the  tube  wall  must  be  related  or  specified.  Here  we  assume  a  constant  heat  flux  to  the  tube 
(and  a  constant  inlet  temperature)  so  that  the  temperature  of  the  heated  surface  is  a  dependent 
variable.  Figure  2.1  illustrates  this  situation.  The  alternate  situadon  of  a  fixed  wall 
temperature  with  heat  flux  as  the  dependent  variable  is  not  of  great  interest  to  spacecraft 
thermal  management  systems. 

Heat  Transfer  Regimes.  Heat  transfer  occurs  by  several  mechanisms-  in  the  evaporator 
tube  because  the  fluid  quality  increases  along  the  tube  length.  The  flow  rates  of  each  phase 
along  the  tube  depend  on  the  heat  transfer  which,  in  tum,  depends  on  the  flow  regime.  For  the 
simple  case  of  a  constant  heat  flux  input  along  the  length  of  the  tube,  the  quality  can  be 
calculated  by  a  simple  energy  balance.  The  major  unknown  variable  in  this  case  is  the  wall 
temperature. 

Figure  2.2  illustrates  all  of  the  two-phase  heat  transfer  regimes.  In  a  forced-convection 
evaporator  the  regimes  will  appear  in  the  order  indicated  by  the  figure: 

•  single-phase  liquid, 

•  subcooled  nucleate  boiling, 

•  saturated  nucleate  boiling,  and 

•  two— phase  forced  convection  heat  transfer. 

Liquid  deficient  and  single— phase  vapor  regimes  have  poor  heat  transfer  and  will  not  appear  in 
a  well-designed  evaporator,  though  they  are  shown  in  the  figure  for  completeness. 


Figure  2.1.  EVAPORATOR  TUBE  DESIGN  INPUTS 
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Single-Phase  Liquid.  Single-phase  heat  transfer  occurs  when  subcooled  liquid  enters  the 
evaporator  and  the  wall  temperature  is  too  low  to  nucleate  bubbles.  This  regime  is  treated 
using  standard  correlations  which  relate  the  heat  transfer  coefficient  to  fluid  properties  and  the 
single-phase  Reynolds  number. 

Nucleate  Boiling.  The  two  nucleate  boiling  regimes  occur  between  the  single-phase 
liquid  and  the  two-phase  forced  convection  regimes.  The  wall  temperature  must  be  high 
enough  to  nucleate  bubbles,  and  the  void  fraction  must  be  small  enough  so  that  an  annular 
flow  pattern  is  not  established.  Heat  transfer  coefficients  are  very  high,  however  the  total 
amount  of  heat  transfer  is  usually  small  since  this  regime  exists  over  only  a  short  length  of  the 
tube.  As  described  in  the  detailed  design  section  below,  there  are  actually  two  types  of 
subcooled  nucleate  boiling,  fully  developed  and  partially  developed,  which  differ  in  the  extent 
of  nucleation  on  the  tube  wall. 

Two-phase  Forced  Convection.  For  most  cases  of  practical  interest  this  is  the  only 
important  heat  transfer  regime.  The  major  design  parameters  for  an  evaporator  tube  can  be 
estimated  by  assuming  that  the  entire  tube  is  characterized  by  this  regime.  The  design  maps  in 
Section  2.1.2  use  this  approach. 

In  this  regime  the  fluid  flows  in  an  annular  pattern  with  a  vapor  core  and  a  liquid  film  on 
the  tube  wall.  Vaporization  occurs  by  two  mechanisms:  nucleate  boiling  in  the  liquid  at  the 
wall  and  evaporation  at  the  vapor/iiquid  interface.  Two— phase  forced  convection  is  the 
dominant  mechanism  of  heat  transfer  for  practical  coolants  because  the  void  f.„iion  is  quite 
high  even  at  low  qualities,  implying  an  annular  flow  regime  over  most  of  the  tube  length. 

In  tnis  Design  Manual  the  "saturated  nucleate  boiling"  regime  (Figure  2.2)  is  considered 
a  part  of  the  two-phase  forced  convection  regime  because  they  are  treated  identically  in  most 
cases.  This  regime  of  non— annular  (bubbly  or  slug)  flow  is  typically  very  small  and  has  little 
effect  on  the  overall  design. 

A  useful  concept  for  describing  heat  transfer  in  the  two-phase  forced  convection  regime 
is  the  superposition  principle.  Heat  transfer  is  treated  as  the  sum  of  two  components.  One  is 
due  to  nucleate  boiling  at  the  wall,  and  the  other  is  due  to  forced  convection  heat  transfer  to 
the  iiquid  film.  Quantitatively,  the  total  heat  flux,  q",  is  expressed  as  a  sum: 

4  =  (hfc  +  hnb,)  ATsal  i) 

where  hfc  is  the  "forced  convection"  heat  transfer  coefficient.  hnb  is  the  "nucleate  boiling"  heat 
transfer  coefficient,  and  ATsal  is  the  temperature  difference  between  the  tube  wall  and  the 
fluid. 


The  superposition  principle  is  the  basis  for  an  important  evaporator  design  criterion. 
Since  forced  convection  heat  transfer  is  independent  of  gravity  while  nucleate  boiling  heat 
transfer  has  an  uncertain  gravity  dependence,  a  dr  in  criterion  for  evaporators  in  microgravitv 
is: 


hfc  »  tint.  (2-2) 

In  this  way  operation  of  the  evaporator  is  assured  regardless  of  the  spacecraft  acceleration. 
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Typical  Behavior.  Figure  2.3  shows  predicted  fluid  and  wall  temperature  profiles  in  one 
particular  forced  convection  evaporator  tube.  This  temperature  plot  applies  to  a  6  mm 
evaporator  tube  with  a  uniform  wall  heat  flux  of  2.7  W/cm2.  The  coolant  is  Freon  R— 1 14 
which  enters  with  a  subcooling  of  3.3°C  and  a  flow  rate  of  5.66  g/s.  In  the  single-phase  liquid 
regime,  the  coolant  temperature  increases  steadily  due  to  the  constant  heat  input.  The  wall 
temperature  remains  a  fixed  amount  above  the  fluid  temperature  because  the  single-phase  heat 
transfer  coefficient  is  constant.  In  the  nucleate  boiling  regime,  the  coolant  temperature 
continues  to  rise  but  the  wall  temperature  drops  because  heat  transfer  is  very  effective.  After 
the  fluid  reaches  its  saturation  temperature,  the  two-phase  forced  convection  regime  begins. 
The  wall  temperature  increases  at  first,  as  bubble  nucleadon  is  suppressed,  but  falls  along  the 
length  of  the  tube  as  interaction  between  the  liquid  and  high  velocity  vapor  act  to  increase  the 
heat  transfer  coefficient. 

Figure  2.4  shows  the  predicted  quality  and  void  fraction  in  the  same  6  mm  evaporator 
tube.  The  quality  increases  linearly  over  the  entire  tube  length  because  the  heat  flux  is 
constant.  The  void  fraction  increases  much  more  quickly  because  the  vapor  density  is  much 
less  than  the  liquid  density.  Thus  the  volume  fraction  of  the  vapor  is  larger  than  its  mass 
fraction. 
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Figure  2.3.  TEMPERATURE  PROFILES  IN  A  FORCEl>-CONVECriON  EVAPORATOR  TUBE 
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Figure  2.4.  QUALITY  AND  VOID  FRACTION  IN  A  FORCED  CONVECTION  EVAPORA'DR 


Depending  on  the  heat  flux  level,  some  of  the  heat  transfer  regimes  may  not  occur  in  an 
evaporator  tube.  At  heat  fluxes  which  are  high  relative  to  the  liquid  mass  flux,  various  other 
regimes  may  occur,  as  described  by  Collier  (1986).  However,  these  cases  represent  poor 
evaporator  designs  in  which  the  tube  walls  dry  out  and  reach  very  high  temperatures.  This 
discussion  is  limited  to  well  designed  evaporators  in  which  two— phase  forced  convection  is  the 
major  heat  transfer  mechanism. 

As  will  be  shown  in  Section  2.2.  pressure  drops  in  evaporator  tubes  are  usually  quite 
small  compared  to  the  absolute  pressure  or'  the  cooiant.  As  a  result  the  cooiant  properties  .viil 
not  change  appreciably  along  the  tube. 

Readers  who  wish  funner  background  .n  two— phase  neat  transfer  should  consult  the  text 
by  Collier  ( 1986;  or  the  report  by  Eastman  (,1985). 

2.1.2  Design  Maps  and  Scoping  Calculations 

This  section  presents  methods  for  scoping  calculations  of  the  evaporator  size  and  heat 
transfer.  It  is  assumed  that  the  heat  flux  and  the  fluid  conditions  at  the  inlet  are  known  and 
that  the  primary  design  goals  are  to  achieve  a  desired  fluid  exit  quality  and  to  ensure  operation 
independent  of  spacecraft  acceleration.  The  basic  procedure  to  meet  these  objectives  is: 
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•  Estimate  the  evaporator  tube  length, 

•  Estimate  the  forced  convection  fraction,  and 

•  Check  the  design  for  maximum  heat  flux  and  tube  wall  temperature. 

Two— phase  forced  convection  with  annular  flow  is  the  dominant  regime  in  an  evaporator  tube 
if  the  inlet  subcooling  of  the  fluid  is  only  a  few  degrees.  The  following  discussion  assumes 
that  this  regime  dominates  the  evaporator  behavior.  For  detailed  design  all  of  the  heat  transfer 
regimes  must  be  treated  explicitly,  and  the  detailed  analysis  presented  in  Section  2.1.3  is 
recommended. 

The  information  needed  in  order  to  perform  the  scoping  calculations  includes: 


Fluid  Properties 

Pi  density  of  liquid  phase  (kg/m3) 

(i  i  viscosity  of  liquid  phase  (kg/m-s) 
cpI  specific  heat  of  liquid  phase  (J/kg — K) 
a  surface  tension  (N/m) 

h  fg  latent  heat  of  evaporation  (J/kg) 

pg  viscosity  of  vapor  phase  (kg/m— s) 
pg  density  of  vapor  phase  (kg/m3) 

Operating  Conditions 

Tsat  saturation  temperature  (K) 

ATsub  inlet  subcooling  (K) 

G  mass  flux  (kg/m2— s) 

P  evaporator  absolute  pressure  (Pa) 

x0  evaporator  outlet  quality  (-) 


Geometry 

H  tube  heated  perimeter  (m) 

A  cross-sectional  flow  area  (m2) 


Evaporator  Tube  length.  Figure  2.5  is  a  universal  design  map  (independent  of  heat 
transfer  regime)  which  shows  the  required  evaporator  tube  length  in  terms  of  the  desired 
change  in  fluid  quality  and  an  evaporator  thermodynamic  ratio,  Rev,  which  is  proportional  to 
the  heat  flux  divided  by  the  mass  flux.  A  compact  evaporator  has  a  large  thermodynamic 
ratio,  unfortunately  this  competes  with  design  requirements  for  forced— convection  dominance 
and  low  pressure  drops.  The  optimum  Rev  comes  from  a  tradeoff  between  these  objectives. 


5  EVAPORATOR  THERMODYNAMIC  RATIO,  R0V 


Figure  2.5.  EVAPORATOR  TUBE  LENGTH 


The  procedure  to  use  this  map  includes  three  steps  as  described  below. 

1)  First,  calculate  the  thermodynamic  ratio,  Rev,  for  the  evaporator: 

R  =  4q."  r 

Kcv  “  u“ 

where: 

q '  =  the  evaporator  tube's  (uniform)  heat  dux  (W/mA,  and 
G  and  h  are  defined  above. 

2)  Next,  determine  the  change  in  fluid  quality.  Ax,  whicii  the  evaporator  is  to  produce. 
The  inlet  subcooling  and  exit  quality,  and  thus  the  change  in  quality  within  the  evaporator,  are 
determined  as  pan  of  the  overall  system  design.  If  fluid  enters  the  evaporator  with  a 
subcooling,  ATSUb,  and  is  to  exit  with  a  quality  x0,  then  the  change  in  enthalpy  is: 


Ah  —  cpiATsub  +  x0hfg 


(2-1) 


The  change  in  quality  in  the  evaporator  is  calculated  by: 


Ax  =  Ah/hfg  (2-5) 

3)  Select  the  hr. a  which  corresponds  to  Ax  in  Figure  2.5  and  locate  the  point  with  the 
abscissa  Rev.  The  value  of  the  ordinate  at  this  point  is  the  dimensionless  tube  length,  L*, 
which  will  produce  the  desired  exit  quality.  The  required  tube  length  is: 

L  =  L*Dh  (2-6) 


where: 

L  =  tube  length  (m),  and 

Dh  =  4A/H  (the  hydraulic  diameter)  (m). 


Forced  Convection  Fraction.  The  evaporator  design  should  ensure  that  forced  convection 
is  the  dominant  heat  transfer  mechanism.  The  fraction  of  heat  transfer  due  to  forced 
convection  in  the  two— phase  forced  convection  regime  depends  on  the  coolant  properties,  mass 
flux,  tube  hydraulic  diameter,  heat  flux  and  the  quality.  Figures  2.6  through  2.9  are 
dimensionless  design  maps  which  allow  the  designer  to  estimate  the  forced  convection  fraction 
for  any  situation.  Figures  2.10  and  2.11  show  specific  design  maps  for  Freon  11  and  Freon 
1 14  in  a  6  mm  diameter  tube  with  a  heat  flux  of  2  W/cm2.  Though  the  first  two  maps  are 
more  general,  they  are  considerably  more  difficult  to  use  than  the  specific  maps.  Additional 
specific  maps  can  be  generated  using  the  procedures  detailed  in  Section  2.1.3  which  follows. 

The  forced  convection  fraction  for  any  fluid  in  any  straight  tube  evaporator  can  be 
calculated  using  Figures  2.6  through  2.9.  These  figures  permit  the  designer  to  estimate  the 
forced  convection  and  nucleate  boiling  heat  transfer  coefficients  (h  fc  and  hnb)  as  a  function  of 
fluid  properties  and  condenser  operating  conditions.  The  forced  convection  fraction,  Efc,  is 
simply  a  ratio  of  these  heat  transfer  coefficients. 

There  are  10  steps  to  calculate  the  forced  convection  fraction: 

1.  Select  a  representative  fluid  quality,  x. 

2.  Calculate  the  Martinelli  parameter.  X,t  (Equation  2—7). 

i.  Determine  the  Reynolds  numoer  factor,  F  i Figure  2.6). 

4.  Calculate  the  two— phase  Reynolds  number,  Re^  (Equadon  2—8). 

5.  Determine  the  forced— convection  heat  transfer  coefficient,  hrc  (Figure  2.7  and 
Equation  2-9). 

o.  Determine  the  suppression  factor.  S  (Figure  2.8). 

7.  Estimate  the  wall  superheat,  ATsat  (Equations  2—10  and  2-1 1). 

8.  Determine  the  nucleate  boiling  heat  transfer  coefficient,  hnt)  (Equation  2-12  or 
Figure  2.7). 

9.  Check  the  heat  flux,  and  iterate  (return  to  Step  8)  if  necessary  (Equation  2-171. 

10.  Calculate  the  fraction  of  the  heat  transfer  due  to  forced  convection,  Efc  (Equation 

2-19). 
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Figure  2.9.  DESIGN  MAP  FOR  NUCLEATE  BOILING  HEAT  TRANSFER  IN  TWO-PHASE 
FORCED  CONVECTION 
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These  steps  are  described  in  more  detail  below. 


1)  Select  a  representative  fluid  quality,  x:  The  average  flow  quality  in  the  channel 
will  do.  This  is  simply  one-half  the  exit  quality. 


2)  Calculate  the  Martinelli  parameter,  X  u: 


Xu  = 


ri-x1 

0.9 

0.5 

Ei 

rt£i 

X 

Lp  ij 

LM-gJ 

3)  Determine  the  Reynolds  number  factor:  Th 
function  of  1/Xtl  by  referring  to  Figure  2.6. 

4)  Calculate  the  two— phase  Reynolds  number: 


0.1 

(2-7) 

'  factor  may  be  read  directly  as  a 


Re,p  = 


G(l-x)Dh  1 -25 

- F 

Hi 


(2-8) 


5)  Determine  the  forced— convection  heat  transfer  coefficient,  hfc.  Use  Figure  2.7 
and  read  the  value  of  the  Nusselt  number  for  forced  convection  (Nufc)  which  corresponds  to 
the  two-phase  Reynolds  number  and  the  liquid  Prandtl  number  (Pr[  =  [ijCpi  /k|).  Then 
calculate  the  heat  transfer  coefficient: 


h  f,. 


Nufcki 


Dv 


(2-9) 


6)  Determine  the  suppression  factor  S:  S  is  a  function  of  Re,p  and  can  be  read 
directly  from  Figure  2.8. 


7)  Estimate  the  wall  superheat:  The  wall  superheat  can  then  be  estimated  by 


"0.5 


•X  r sat  fc  /2hnb )-  t  tq  /hnb )  |  (h fc  /-hnb ) 


2—10) 


where 


hnb'  -  0.00122 


k  1:) •  79C p  jO.  J5p  ,0.4 9hr  0.51 


(2-1 1) 


Tsat0.7500.5p,0.29pg0.24Vfg0.75 

(which  is  approximately  hnb/ATsal). 

8)  Determine  the  nucleate  boiling  heat  transfer  coefficient:  The  easiest  way  to  do 
this  is  simply  to  calculate 


hnb  —  ^rib  ^Tsat 


(2-12) 


Alternatively,  Figure  2.9  can  be  used  to  read  the  Nusselt  number  for  nucleate  boiling, 
Nunb.  The  design  map  in  Figure  2.9  is  relatively  impractical  to  use  because  the  heat  transfer 
coefficient  is  a  function  of  several  dimensionless  parameters,  and  Equation  2—12  is 
recommended  instead.  Nevertheless,  Figure  2.9  is  useful  because  it  shows  that  this  procedure 
scales  to  any  fluid.  If  this  figure  is  used,  the  Nunb  is  located  which  corresponds  to  the  liquid 
Prandtl  number  and  the  bubble  Reynolds  number,  Reb: 

1  Cpl  Pg^Tsa 

Reb  = - 

Pi  L  Pihfg 

The  heat  transfer  coefficient  for  nucleate  boiling  is  then: 

Nunbki 

hnb  = -  (2-14) 

ib 

where  rb  is  a  characteristic  bubble  radius  given  by: 

ATsat  [2itk  ip |Cp|Cl £  f  PiH  o . 25 
rb  = - —  -  S  (2-15) 

h  fgPg  .  ^^sat  .  A^satj 

and  APsa[  is  the  difference  in  saturation  pressures  which  corresponds  to  the  wall  superheat: 

APsat  =  ATsat - — -  (2-16) 

Tsat(Vg-Vi) 

where  vg  and  vj  are  the  specific  volumes  of  the  gas  and  liquid  phases,  respectively. 

9)  Check  the  heat  flux  and  iterate  if  necessary:  Verify  that  the  estimated  wall 
superheat  is  accurate  by  comparing  the  known  heat  flux  with  the  heat  flux  calculated  from 
ATsat,  hfc  and  hnb.  Calculate  this  ratio: 

(h  fc  +  hnb;ATsat 

P  = -  12-17) 

q" 

There  are  now  two  possibilities: 

a)  P  is  acceptably  close  to  1.0  -  in  this  case,  go  on  to  step  10  to  calculate 
Elc.  "Acceptably  close"  depends  upon  the  level  of  precision  desired  in 
the  design.  Since  the  Chen  correlation  is  only  accurate  to  within  \0%  of 
the  actual  heat  flux,  it  is  usually  acceptable  to  go  on  if  P  lies  between 
0.9  and  1.10. 


b)  If  p  is  outside  the  acceptable  range,  modify  the  estimate  for  ATsat  and  go 
back  to  step  8.  A  good  next  guess  for  the  superheat  is: 


AT 

^  1  sat  - 


f 


i_hfc  +  ^nbj 
62 


(2—18) 


10)  Calculate  the  forced  convection  fraction:  The  fraction  of  the  heat  transfer  due  to 
forced  convection,  Efc,  is: 


Efc  — 


h  fc 

h  fc  +  hnb 


This  ratio  should  be  near  unity  (E  fc  >  0.9)  for  forced  convection  to  dominate. 


(2-19) 


Figure  2.10  shows  the  forced  convection  fraction  as  a  function  of  flow  quality  and  mass 
flux  for  the  specific  design  case  of  an  R— 1 1  evaporator  operating  at  a  pressure  of  227  kPa,  a 
heat  flux  of  2  W/cm2  and  with  a  diameter  of  6  mm.  Figure  2. 1 1  shows  the  same  information 
for  R— 1 14  at  434  kPa.  These  design  charts  are  simpler  to  use  than  the  dimensionless  maps; 
however,  to  generate  such  plots,  the  designer  must  perform  the  calculations  outlined  in  Section 
2.1.3  which  follows. 


Maximum  Heat  Flux.  There  are  two  limits  which  should  be  checked  for  evaporator 
design.  One  is  the  maximum  heat  flux  and  the  other  is  the  wall  temperature.  To  first  order,  if 
the  fraction  £fc  is  greater  than  0.9,  operation  should  be  possible. 

High  heat  fluxes  tend  to  increase  the  importance  of  nucleate  boiling  within  the 
two— phase  forced  convection  regime.  For  a  given  fluid,  mass  flux  and  tube  diameter  there  is  a 
maximum  value  of  the  heat  flux  above  which  forced  convection  is  no  longer  the  dominant 
mechanism  of  heat  transfer.  For  an  evaporator  tube  which  is  to  operate  in  microgravity  this 
value  is  a  limit  to  the  heat  flux  which  should  not  be  exceeded.  There  are  also  heat  flux  limits 
due  to  critical  heat  flux  considerations  which  are  not  included  in  this  Design  Manual  but  which 
may,  in  some  cases,  be  more  restrictive  than  these. 

Figure  2.12  provides  the  designer  with  quick  checks  on  potential  evaporator  designs  with 
Freon.  To  operate  in  microgravity,  the  tube  heat  flux  should  be  less  than  the  maximum  values 
indicated  by  Figure  2.12.  The  design  map  shown  in  Figure  2.12  presents  limiting  heat  flux 
values  for  R— 1 1  in  tubes  of  various  diameters.  The  mass  flux  is  expressed  as  a  Reynolds 
number  for  the  flow  rate  of  liquid  at  the  inlet: 


Ref  =  GDty'fJ-i  ;  2—20) 

The  heat  flux  indicated  by  Figure  2.12  is  the  value  above  which  nucleate  boiling  accounts  for 
over  10%  of  the  heat  transfer  along  the  tube  length.  Note  that  there  is  roughly  a 
second— power  relationship  between  the  maximum  heat  flux  and  the  mass  flux: 


Figure  2.12  also  shows  the  same  calculation  for  R-l  14.  Section  2.1.3  details  how  to  produce 
similar  design  maps  for  other  coolants. 


Wall  Temperature.  The  evaporator  designer  may  be  faced  with  temperature  limits  for 
the  tube  material.  Figure  2.13  shows  the  wall  temperatures  which  correspond  to  the  maximum 
heat  fluxes  illustrated  in  the  previous  design  maps.  The  wall  temperature  remains  fairly 
constant  until  a  critical  value  of  the  heat  flux  is  reached.  Beyond  this  point  the  wall 
temperature  increases  rapidly.  Design  maps  for  additional  coolants  can  be  prepared  using 
procedures  in  Section  2.1.3. 
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Figure  2.11.  FORCED  CONVECTION  FRACTION  FOR  R-I14  IN  A  6  mra  EVAPORATOR  TUBE 


MAX 


100.000 


cr 


I.COO 


'0.000 


ICC.C! 


INLET  REYNOLDS  NUMBER,  Ref 


Figure  2.12.  MAXIMUM  HEAT  FLUX  FOR  OPERATION  IN  MICROGRAVITY 


Uncertainty.  The  main  uncertainty  in  the  heat  transfer  calculation  is  the  effect  of 
microgravity  on  nucleate— boiling  heat  transfer.  The  only  models  for  forced-convection  heat 
transfer  have  been  developed  based  on  experiments  conducted  under  normal  Earth  gravity. 
The  recommended  approach  is  to  design  the  evaporator  so  that  the  controlling  phenomenon  is 
forced-convection  heat  transfer.  The  designer  should  address  this  uncertainty  by  requiring  a 
high  fraction  of  the  evaporator  heat  transfer  (e.g.  90%)  to  be  due  to  forced  convection. 
Ongoing  microgravity  experiments  will  provide  data  so  that  more  specific  guidelines  can  be 
used  in  the  future. 


2.1.3  Equations  for  Heat  Transfer 

To  first  order,  an  evaporator  tube  can  be  designed  assuming  only  one  heat  transfer 
regime:  two— phase  forced  convection.  However,  for  final  design  calculations  more  precision 
is  often  needed,  and  the  detailed  behavior  of  the  coolant  must  be  considered.  These 
calculadons  should  consider  all  the  possible  heat  transfer  regimes  and  calculate  the  transition 
from  one  regime  to  the  next. 

This  section  presents  detailed  design  relations.  Procedures  are  given  to  generate  design 
maps  for  additional  cases  of  interest  or  to  program  computerized  calculations  for  forced 
convection  evaporators.  The  usual  way  to  use  these  design  relations  is  to  first  divide  the 
evaporator  tube  into  a  number  of  segments  within  which  the  fluid  properties  are  approximately 
constant.  In  each  segment  there  are  three  basic  calculations: 

•  Determine  the  fluid  state  at  the  segment  exit, 

•  Calculate  the  wall  temperature,  and 

•  Check  for  heat  transfer  regime  transition. 

Table  2.1  summarizes  the  equations  for  these  calculations  in  each  of  the  heat  transfer 
regimes.  Figure  2.14  shows  schematically  the  sequence  for  performing  these  calculations. 

The  sections  which  follow  provide  more  detail. 

One  hundred  or  so  segments  are  usually  sufficient  for  the  evaporator  calculation. 
Conditions  in  the  first  segment  are  known  from  the  evaporator  inlet  conditions  and  will  usually 
be  single -phase  liquid  at  a  known  pressure  and  subcooling.  Exit  conditions  calculated  for  one 
segment  become  the  iniet  conditions  for  the  next.  The  following  sections  show  how  to 
calculate  the  fluid  state,  wall  temperature  and  regime  transitions  in  each  of  the  heat  transfer 
regimes  experienced  by  the  fluid  as  it  travels  through  the  evaporator:  single-phase  liquid, 
partial  and  fully  developeu  nucleate  boiling,  and  two— phase  forced  convection. 


Single  Phase  Liquid  Recime.  In  this  regime,  vapor  bubbles  are  not  present  at  all. 


Fluid  State.  In  the  single -phase  regime  all  of  the  heat  is  absorbed  by  the  liquid  as 
sensible  heat.  Therefore  the  temperature  rise,  ATf,  in  a  segment  of  length  A z  is  calculated 
from  a  simple  energy  balance: 


ATf 


4  q"  Az 


(2-22) 
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Figure  2.14.  FLOWCHART  FOR  DFTA1LED  EVAPORATOR  CALCULATIONS 
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Wall  Temperature.  The  evaporator  wall  temperature  in  the  single-phase  regime  depends 
on  the  heat  flux  and  the  single-phase  heat  transfer  coefficient,  hfc: 


Tw  =  Tf  +  q"/hf 


(2-23) 


For  turbulent  flow,  the  Dittus— Boelter  equation  (Collier,  1986)  yields  the  heat  transfer 
coefficient: 


hfc  =  0.023  (k,  /Dh)  (GDh/^0'8  (PrO0'4 


(2-24) 


where 


kj  =  thermal  conductivity  of  the  liquid  (W/m— K),  and 
Pri  =  Prandtl  number  of  the  liquid  (Pt[  =  p.jcp i  /k 


Transition  to  nucleate  boiling.  The  single-phase  heat  transfer  regime  ends  when  the 
wall  temperature  reaches  a  critical  value,  Twonb,  at  which  vapor  bubble  nucleation  and 
subcooled  nucleate  boiling  begins.  The  correlation  of  Davis  and  Anderson  (1966): 


r  —  t  4-  8  o  q  T sat 

w°nb  “  +  [  h  fgk ! p  g  J 


(2-25) 


is  recommended,  or  that  of  Frost  and  Dzakowic  (1967),  which  is  identical  except  that  the 

1/3 

second  term  is  multipl’cd  by  P^ 

This  transition  criterion  can  be  expressed  in  dimensionless  form  as  indicated  in  Table 
2.1.  The  dimensionless  criterion  is: 


k  i(Tw  Tsal) 


9  ^  crit 


-v  fgT'sat^'s 

=  !  1  +  - 

L  ^  fg1"  c  r  it'9 


(2-26) 


where  rrri,  is  a  critical  bubble  radius: 


-v.,T.j3  ,ak;  » 


1 2-271 


and  vf  =  ( l/p„j  -  ( 1/p  ij. 


Partially  Developed  Nucleate  Boiling  Regime.  The  partially  developed  nucleate  boning 
regime  is  the  first  to  occur  following  the  single -phase  liquid  regime.  It  is  characterized  by 
incomplete  activation  of  nucleation  sites  on  the  evaporator  so  that  some  fraction  of  the  heat 
transfer  occurs  by  normal  single— phase  convection  between  patches  of  bubbles.  Vapor 
bubbles  appear  at  the  wall  surface  but  do  not  detach  and  enter  the  flow. 


Fluid  State.  In  this  regime  there  is  no  net  vapor  generation,  so  that  fluid  heating  is  the 
same  as  single-phase,  Equation  2-22. 
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Wall  Temperature.  The  recommended  analytical  approach  is  to  treat  the  heat  flux  in  this 
regime  as  the  superposition  of  a  subcooled  boiling  heat  flux  and  a  single-phase  forced 
convection  heat  flux.  Bowring  (1962)  suggests  a  simple  additive  approach: 

q"  =  q'sub  +  q”fc  (2-28) 

=  q"  sub  +  h  fc(T  sa  f-T  f)  (2-29) 

The  subcooled  boiling  heat  flux  is  calculated  using  the  correlation  of  Thom  (1965): 

q  "sub  =  1.95x10-3  (Tw-Tsat)2  exp(0.023P)  (2-30) 

(where  T  and  P  must  be  expressed  in  SI  units  —  temperature  in  K  and  pressure  in  bars).  The 
single— phase  heat  flux  is  calculated  using  the  Dittus— Boelter  relation  (Equation  2—24). 
Combining  Equations  2-24,  2—29,  and  2-30  and  solving  for  the  wall  temperature  gives: 

Tw  =  Tsat  +  [q”— hfc(Tsa£-T^  [22.65exp(— 0.014P)]  (2-31) 

With  Tw  known,  the  subcooled  boiling  and  forced  convection  heat  fluxes  are  calculated 
using  Equations  2—30  and  2-28,  respectivefy. 

Transition  from  partially  to  fully  developed  subcooled  boiling.  Saha  and  Zuber  (1974) 
recommend  a  model  for  the  transition  from  partially  to  fully  developed  subcooled  boiling 
(FDB).  This  model  expresses  the  fluid  subcooling  at  the  transition  as  a  function  of  the  heat 
flux: 

(Tsar-T^fdb  =  (0.0022)  q"  DAi  (2-32a) 

(when  Re  f  Pr ]  <  70,000),  or 

q" 

(TSat-Tf)fdb  -  053.8) -  (2— 32b) 

Gcpl 

(when  Re,-  Pr i  >  70.000) 

The  transition  criterion  can  also  be  written  in  dimensionless  form: 

455 

St  = -  Re,  Pr,  <  70.000  :2-33a) 

Re  fPr , 

St  =  6.5  x  10-3  Ref  Pr,  >  70,000  (2— 33b) 

where  St  is  the  Stanton  number: 

Nu 

St  = -  (2-34) 

Re  fPr, 


Fully  Developed  Nucleate  Boiling  Regime.  In  this  regime  there  is  net  vapor  generation. 
Bubbles  detach  from  the  wall  and  flow  with  subcooled  liquid. 
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Fluid  State.  In  this  regime  there  is  no  division  of  the  heat  flux  between  different  heat 
transfer  mechanisms  since  the  entire  surface  of  the  tube  is  covered  with  bubbles.  Nevertheless 
the  liquid  may  still  be  subcooled  and,  if  so,  a  fraction  of  the  heat  addition  will  increase  the 
liquid  temperature.  This  occurs  by  condensation  of  vapor  bubbles  in  the  subcooled  liquid 
away  from  the  wall. 

Ivey  and  Morris  (1962)  suggest  that  the  fraction  of  the  heat  flux,  Ef,  which  goes  to 
heating  the  fluid  is: 


1  +  0.1(pi  /pg)°-75cpi(Tsar-T£)/hfg 


-1 


(2-35) 


In  this  case  the  fluid  temperature  rise  in  a  segment  is: 

AT'=  <2-36' 

Note  that  the  bulk  fluid  temperature  will  not  exceed  the  saturation  temperature,  Tsal. 
Once  the  liquid  reaches  Tsat ,  the  heat  transfer  regime  changes  to  two— phase  forced 
convection,  as  discussed  below. 

Wall  Temperature.  Thom  (1965)  recommends  this  correlation  for  the  wall  superheat  in 
fully  developed  boiling: 


(Tw— Tsat)  =  22.65  (q")°-5  exp(— 0.014P)  (2-37) 

(where  T  and  P  must  be  expressed  in  SI  units  —  temperature  in  K  and  pressure  in  bar). 

Transition  from  subcooled  to  saturated  boiling.  The  nucleate  boiling  regime  ends  when 
the  two— phase  forced  convection  regime  begins.  To  be  precise,  one  could  calculate  a 
transition  to  annular  flow  based  on  the  local  void  fraction,  and  this  would  represent  the 
transition  to  two— phase  forced  convection  heat  transfer.  However  the  transition  void  fraction 
is  exceeded  almost  immediately  as  soon  as  the  two— phase  fluid  achieves  a  positive  quality.  So 
it  is  acceptably  accurate  and  simpler  to  assume  that  the  nucleate  boiling  regime  ends  when  the 
liquid  reaches  its  saturation  temperature: 


Tf  =  Tsal  (2-38) 

Two— Phase  Forced  Convection  Regime.  Typically  most  of  the  evaporator  will  operate 
in  this  regime.  This  section  will  first  give  procedures  for  calculating  the  fluid  state  and  wail 
temperature.  The  generation  of  design  maps  for  forced  convection  fraction,  maximum  heat 
flux  and  maximum  wall  temperature  (such  as  Figures  2.6  through  2.1 1)  will  also  be  described 


Fluid  State.  In  this  regime,  ail  of  the  heat  flux  is  used  to  generate  vapor  and  increase  the 
quality  of  the  two— phase  mixture.  The  change  in  quality,  Ax,  in  a  segment  of  length  of  Az  is: 


To  calculate  the  pressure  drop  in  the  evaporator  tube  (Section  2.2)  it  will  be  useful  to 
have  an  expression  for  dz/dx.  Rearranging  Equation  2—39  and  taking  the  limit  as  Az 
approaches  zero  gives: 


dz  Gh  fgDh 
dx  4  q" 


(2-40) 


Wall  Temperature.  The  wall  temperature  is  calculated  using  the  correlation  developed 
by  Chen  (1963).  There  are  eight  steps  described  below. 

1.  Calculate  the  Martinelli  parameter,  Xtt: 


Xtt  = 


1— x 


0.9 


0.5  r..  -i 0.1 


ite' 


(2-41) 


2.  Calculate  the  Reynolds  number  factor,  F.  F  can  be  determined  either  by  examination 
of  Figure  2.6,  which  shows  F  as  a  function  of  1/Xtt,  or  by  using  the  following  approximate 
equations: 


=  1.0  (for  1/Xtt<  0.1) 

=  2.35  (0.213  +  1/Xtt)0.736  (f0r  i/xtt>0.1) 

3.  Calculate  the  effective  two— phase  Reynolds  number,  Re^,. 

Re  =  Fi .  25 

v  4 1 

4.  Calculate  the  forced  convection  heat  transfer  coefficient,  hfc: 

hfc  =  (0.023)(ky'Dh;(Retp/>-SfPr1)0-33  F 


(2-42) 


(2-43) 


(2-44) 


5.  Calculate  the  suppression  factor.  S.  S  can  be  read  directly  from  Figure  2.3,  where  it 
is  shown  as  a  function  of  the  two— phase  Reynolds  number.  Or  S  may  be  calculated  from  these 
approximate  equations: 

/ 

0.90  Re^  <  1.5x1 04  (2-45) 


S  =  <  3.09  -0.234  ln(Re,p) 


1.5x10*  <  Re^  <  3.5x105 


0.1 


3.5x10s  <  Req, 
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6.  Calculate  the  factor  hnb',  which  is  equal  to  the  nucleate  boiling  heat  transfer 
coefficient  divided  by  ATsat°-": 

kj0-79c_j0.45pj0.49hf  0.51 

hnb'  =  0.00122  - S  (2-46) 

Tsa  ,0  ■  75(j0 . 5p  ,0 .  29pg0 . 24v  fg0 . 75 

(Note  that  this  correlation  by  Chen  suggests  that  the  nucleate  boiling  component  of  heat 
transfer  is  independent  of  gravity.  However,  other  workers  such  as  Shah  (1976)  maintain  that 
gravity  effects  are  significant.) 


7.  Solve  for  ATsat: 

ATsati-99hnb'  +  ATsathfc-q"  =  0  (2-47) 

To  be  precise,  Equation  2-47  needs  to  be  solved  numerically;  however,  it  is  usually 
accurate  enough  to  approximate  this  equation  as  a  quadratic  in  ATsal,  in  which  case  the 
solution  is: 

ATsat  =  [(hfc/2hnb')2  +  (q"/hnb')]°-5-(hfc/2hrib’)  (2-48) 

This  approximation  is  accurate  for  the  range  of  heat  fluxes  normally  encountered  in 
evaporators.  If  the  value  of  ATsat  given  by  Equation  2-48  does  not  solve  Equation  2-47  to  an 
acceptable  degree  of  accuracy,  then  a  few  iterations  will  be  enough  to  converge  on  an  accurate 
estimate. 

8.  Calculate  the  wall  temperature,  Tw: 

Tw  =  Tsat  +  ATsat  (2—49) 


Limits  to  Evaporator  Heat  Transfer.  In  a  well-designed  evaporator,  two— phase  forced 
convection  is  the  final  heat  transfer  regime.  Following  two— phase  forced  convection  comes  a 
liquid  deficient  region  with  droplet  flow  and  then  a  single-phase  vapor  regime.  Both  of  these 
regimes  are  charactenzed  by  poor  heat  transfer  coefficients  and  correspondingly  high  wall 
temperatures.  A  substantial  body  of  literature  is  available  on  predicting  the  end  of  the 
two— phase  forced  convection  regime,  which  is  often  terr.  ed  "dryout."  However,  this  analysis 
is  beyond  the  scope  of  this  design  manual  report.  Interested  readers  should  refer  to  Collier 
(1986)  for  a  complete  discussion  of  dryout.  Below,  simple  limits  to  check  are  recommended. 


Forced  Convection  Fraction  A  design  criterion  for  microgravitv  evaporators  is  that  the 
fraction  of  heat  transferred  by  forced  convection  should  be  close  to  1.0.  This  fraction  can  be 
calculated  using  the  above  correlation  for  heat  transfer  in  the  two— phase  forced  convection 
regime. 

The  dimensional  design  maps  in  Figures  2.10  and  2.11  apply  for  specific,  constant 
fi"'d  properties,  tube  diameter  and  heat  flux.  If  these  parameters  are  known,  then  the  above 
procedure  (Equations  2—41  through  2—46)  shows  how  to  calculate  hfc  and  hnb'  as  a  function  of 
the  fluid  quality,  x.  It  will  always  be  acceptable  to  use  the  qr  .dratic  approximation  of 
Equation  2—48  since  high  forced  convection  fractions  correspond  to  low  values  of  ATsat.  The 
fraction  of  heat  transfer  due  to  forced  convection,  Efc,  is  then: 
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Efc(x)  =  = 


1/2  + 


(l/4)+(q”hnb7hfc2) 


0.5 


i  -1 


(2-50) 


Maximum  Heat  Flux.  The  maximum  heat  flux  design  maps  (Figure  2.12)  give  the 
maximum  heat  flux  for  microgravity  operation.  Higher  heat  fluxes  than  indicated  by  these 
maps  require  too  high  a  fraction  of  heat  transfer  by  nucleate  boiling.  The  maps  are  generated 
by  finding  the  heat  flux  which  yields  a  forced  convection  fraction  equal  to  the  design  limit 
(90%  for  Figure  2.12). 

The  maximum  heat  flux  is  found  iteratively  at  each  Reynolds  number.  First  guess  a  heat 
flux  assuming  that  all  heat  transfer  is  by  forced  convection  alone.  Then  calculate  Efc(x)  for  a 
range  of  qualiries  representing  different  positions  along  the  evaporator  tube.  The  average 
value  of  Efc  along  the  tube  should  equal  the  design  limit  (e.g.  90%)  at  the  maximum  heat  flux. 
If  the  average  forced  convection  fraction  is  less  than  the  design  limit,  then  modify  the  next 
guess  for  the  heat  flux  by  decreasing  it,  otherwise  increase  the  heat  flux.  Recalculate  the 
average  Efc  until  it  is  acceptably  close  to  the  design  limit  This  yields  the  maximum  heat  flux. 


Wall  Temperature.  The  maximum  wall  temperature  design  maps  (Figure  2.13)  are 
simple  to  generate  once  the  maximum  heat  transfer  calculations  are  complete.  These  maps 
simply  show  the  average  wall  temperature  which  corresponds  to  the  maximum  heat  flux.  The 
average  wall  temperature  is  simply  the  arithmetic  average  of  the  segment  wall  temperatures 
found  in  the  above  calculation. 


2.1.4  Validation  with  Microgravity  Data 

The  only  experiments  which  have  been  performed  using  a  forced-convection  evaporator 
of  this  type  in  microgravity  are  those  of  Lee  (1987).  Unfortunately,  there  are  few  data  from 
these  tests  which  are  useful  for  validating  the  microgravity  evaporator  models.  These  tests 
were  performed  on  a  KC— 135  aircraft  which  flew  parabolic  trajectories  to  provide  short 
periods  (~30s)  of  reduced  (but  unsteady)  gravity.  The  test  loop  could  not  achieve  steady  state 
under  these  conditions.  Furthermore,  the  test  conditions  and  results  are  not  well  documented. 

Recently,  Cuta  and  Krouuk  ;  1988)  have  studied  heat  transfer  data  from  a  few  selected 
tests  in  an  evaporator  flown  as  part  of  this  experiment.  Code  calculations  of  the  evaporator 
temperatures  are  close  to  the  measured  data  for  two  of  three  tests  examined.  Since  the  model 
for  saturated  boiling  in  the  computer  code  used  in  the  analysis  (COBRA/TRAC)  is  the  same  as 
recommended  here,  there  is  preliminary  evidence  that  the  model  is  satisfactory.  Additional 
experiments  and  analysis  are  needed  for  validation,  however. 


2.2  Pressure  Losses 
2.2.1  Introduction 

The  pressure  drop  through  an  evaporator  determines  the  pumping  power  requirement. 
Generally,  good  heat  transfer  performance  requires  higher  pressure  drops,  so  system  level 
tradeoffs  are  necessary  to  balance  component  size  and  pumping  power.  For  example,  the 
pressure  drop  through  an  evaporator  tube  varies  as  the  mass  flux  cubed  while  the  maximum 
heat  flux  varies  as  the  mass  flux  squared.  Thus,  for  an  evaporator  tube  there  will  be  a  certain 
mass  flux  beyond  which  small  increases  in  heat  transfer  are  accompanied  by  much  greater 
pressure  losses. 
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This  section  presen* dimensionless  design  maps  to  calculate  pressure  losses  in 
evaporator  tubes.  The  dimensionless  design  variables  provide  physical  insight  and  the  maps 
are  suitable  for  scoping  level  design.  These  maps  ignore  the  less  important  heat  transfer 
regimes  which  contribute  only  a  small  fraction  of  the  total  pressure  drop.  More  precise 
estimates  of  the  pressure  drop  account  for  all  the  various  heat  transfer  regimes,  and  design 
equations  are  provided  which  are  suitable  for  detailed  computer  calculations. 

Pressure  losses  in  an  evaporator  tube  are  due  to  the  phenomena  illustrated  in  Figure  2.15, 
which  shows  the  forces  acting  on  elemental  control  volume  in  the  evaporator  tube.  The  forces 
are: 

•  Wall  friction, 

•  Body  forces,  and 

•  Fluid  acceleration. 

The  total  pressure  drop  in  the  evaporator  is  the  sum  of  the  pressure  drops  due  to  these  three 
components.  Pressure  drops  due  to  wall  friction  and  body  forces  are  similar  to  the  pressure 
losses  in  a  transport  line  (Section  1.2).  The  difference  is  that  the  fluid  changes  phase  in  an 
evaporator  so  the  effective  fluid  properties  vary  along  the  tube.  The  pressure  drop  due  to  fluid 
acceleration  is  also  caused  by  vaporization  of  liquid.  Since  the  vapor  is  less  dense  than  the 
liquid,  it  must  have  a  higher  velocity.  To  accelerate  the  fluid  into  the  vapor  phase  a  force 
must  be  exerted  on  the  liquid.  This  force  is  the  acceleration  pressure  loss  in  the  evaporator 
tube. 


The  components  of  the  pressure  drop  in  one  particular  R~1 14  evaporator  tube  are 
illustrated  in  Figure  2.16.  Note  that  almost  all  of  the  pressure  drop  occurs  in  the  two— phase 
forced  convection  regime.  The  friction  and  fluid  acceleration  components  of  the  pressure  drop 
are  roughly  equal.  The  pressure  loss  due  to  body  forces  is  negligible  for  this  case,  in  which 
the  gravity  has  magnitude  O.OOlg  (directed  opposite  the  flow  vector).  The  acceleration 
pressure  drop  increases  linearly  along  the  tube  because  the  heat  flux  is  constant,  which  causes 
a  constant  rate  of  increase  in  quality  with  distance.  Frictional  pressure  losses  are  roughly 
proportional  to  the  square  of  the  fluid  velocity.  The  pressure  gradient  due  to  friction  increases 
along  the  evaporator  tube  because  the  fluid  velocity  increases  as  its  density  is  reduced  by 
vaporization. 


2.2.2  Oesign  Maps 

This  section  presents  design  maps  to  estimate  pressure  losses  in  an  evaporator  tube. 
These  maps  assume  that  the  major  pressure  losses  in  the  tube  are  due  to  the  two— phase  forced 
convection  regime.  If  there  is  a  very  large  subcooling  at  the  entrance  to  the  evaporator,  then 
there  may  be  a  significant  pressure  loss  in  the  single-phase  liquid  regime  as  well. 
Single-phase  methods  can  be  used  if  that  is  the  case. 

The  basic  procedure  to  estimate  the  evaporator  pressure  drop  is: 

1.  Estimate  the  pressure  drop  due  to  friction  in  the  two— phase  region  (Figure  2.17 
and  Equations  2-5 1  to  2-54), 

2.  Estimate  the  pressure  drop  due  to  body  forces  in  the  two— phase  region  (Figure 
2.18  and  Equation  2-55), 

3.  Estimate  the  pressure  drop  due  to  fluid  acceleration  in  the  two— phase  region 
(Figure  2.19  and  Equation  2-56), 
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DISTANCE  FROM  EVAPORATOR  INLET,  z(m) 


Figure  2.16.  PRESSURE  LOSSES  IN  A  VERTICAL  FORCED  CONVECTION  EVAPORATOR 
TUBE  (O.OOl  g) 


4.  Calculate  the  total  pressure  drop  in  the  two— phase  region  (Equation  2—57). 

5.  Estimate  the  pressure  drop  in  the  single-phase  region  (Equations  2-58  and 
2-59). 

6.  Estimate  the  total  pressure  drop  in  the  evaporator  (Equation  2-60). 

The  following  sections  describe  the  procedure  to  obtain  a  best— estimate  calculation  for 
the  pressure  drop  in  the  evaporator. 

The  main  uncertainty  in  these  calculations  is  the  interfacial  friction  between  the  liquid 
and  vapor  phases  in  the  two— phase  forced  convection  region  of  the  evaporator.  The  pressure 
drop  can  be  bounded  by  choosing  high  or  low  values  for  the  interfacial  friction. 


Calculating  Pressure  Losses.  Most  of  the  pressure  drop  typically  occurs  in  the 
two-phase  region.  Procedures  to  calculate  the  three  components  of  the  pressure  drop  are 
given  below.  All  of  the  dimensionless  pressure  drops  are  expressed  as  functions  of  the  change 
in  quality  m  the  evaporator. 

1)  Friction  in  the  Two-Phase  Region.  The  procedure  to  calculate  the  frictional  pressure 
drop  is: 

•  Calculate  the  ratio  of  the  liquid  phase  density  to  the  vapor  phase  density.  Pi  /pg. 

•  Locate  the  curve  in  Figure  2.17  which  corresponds  to  this  value  of  the  density 
ratio  and  locate  the  point  on  the  curve  which  corresponds  to  the  evaporator  outlet 
quality.  Call  the  ordinate  on  the  map  at  this  point  APfe*  /4fw  [. 

•  Calculate  Rev,  the  thermodynamic  ratio  of  the  evaporator: 

4  q" 

Rev  =  -  (2-51) 

G  h  fg 

•  Calculate  fw  the  friction  factor  for  liquid— only  flow  through  the  evaporator. 

This  friction  factor  is  a  function  of  Ref,  the  Reynolds  number  for  liquid— (only 
How: 


GDh 

Re,-  =  -  (2-52) 


Fhe  liquid— only  friction  factor  is: 


16 


-0.2 

[  0.046  Ref 


for  Re  f  <  1 500 


for  Re  f  >  1 500 


(2— 53a) 

(2— 53b) 
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QUALITY,  x 


Figure  2.17.  DIMENSIONLESS  FRICTIONAL  PRESSURE  DROP  IN  THE  EVAPORATOR 


•  Calculate  APfe,  the  pressure  drop  in  the  two— phase  section  of  the  evaporator  due 
to  wall  friction: 

’G2/2pi]  TAPfe*' 

APfe  =  4fwl  -  -  (2-54) 

P  ev  _4  f  w  | 

2)  Jociv  Forces  :n  :ne  Two-Phase  legion.  The  procedure  to  eaicuiate  the  pressure  urop 
due  to  body  forces  is: 

•  Locate  the  curve  in  Figure  2.18  which  corresponds  to  the  iiuuiu-to— vapor  density 
ratio,  p;/p...  Find  the  point  on  this  curve  which  corresponds  to  the  evaporator 
outlet  qualify.  Tire  ordinate  of  this  plot  is  the  dimensioniess  pressure  drop 

•  Calculate  the  dimensional  pressure  drop.  APge: 

p  i  Dh  a  cos©' 

APge  =  - - -  [APgc*]  (2-55) 

3)  Acceleration  in  the  Two-Phase  Region.  The  procedure  to  calculate  the  pressure  drop 
due  to  fluid  acceleration  in  the  evaporator  is: 

•  Locate  the  curve  in  Figure  2-19  which  corresponds  to  the  ratio  of  liquid — to— 
vapor  density,  pi  /pg.  Find  the  point  on  this  curve  which  corresponds  to  the 
evaporator  outlet  quality.  The  ordinate  is  the  dimensionless  pressure  drop  APa*. 
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DIMENSIONLESS  PRESSURE  DROP, 


"GRAVITY"  PRESSURE  LOSSES  IN  EVAPORATOR 


Figure  2.18.  DIMENSIONLESS  PRESSURE  DROP  DUE  TO  BODY  FORCES  IN  THE  EVAPORATOR 


ACCELERATION  PRESSURE  LOSS 


Figure  2.19.  D [.MENS IONLESS  PRESSURE  DROP  DUE  TO  FLUID  ACCELERATION  IN  THE 
EVAPORATOR 


Calculate  the  dimensional  pressure  drop,  APa: 


APae  =- 


G2 


L2Pi 


(APa*) 


(2-56) 


4)  Total  Pressure  Drop  in  Two-Phase  Section.  The  total  pressure  drop  in  the  two— phase 
section  of  the  evaporator,  APq,,  is  the  sum  of  the  frictional,  body  force,  and  fluid  acceleration 
components: 


AP  tp  =  APfe  +  APge  +  APae 


(2-57) 


5)  Pressure  Drop  in  the  Single-Phase  Liquid  Section.  This  pressure  drop  will  be 
insignificant  unless  the  subcooling  is  very  high.  There  are  only  two  components  of  the 
pressure  drop  in  the  single-phase  section:  friction  and  body  forces.  The  total  pressure  drop  is 
easily  estimated  using  parameters  already  calculated  in  the  procedure  for  the  two— phase 
section.  The  procedure  is: 


Calculate  Lsc,  the  length  of  the  evaporator  tube  in  which  single— phase  liquid 
flows.  For  these  scoping  calculations  it  is  sufficient  to  estimate  Lsc  using  this 
equation: 


L 


sc 


■*v 


G  Dhcp  )AT sub 
4  q" 


Estimate  APsp,  the  pressure  drop  in  the  single  phase  section: 


(2-58) 


APsp  = 


fL 


4fwl| 


Dh 


G2  ' 

2  Pi. 


—  pi  Lsc  a  cosG 


(2 -5  9) 


(Equation  2—59  can  also  be  used  for  detailed  calculations  in  which  the  evaporator  is 
subdivided  into  many  short  segments.  L.c  is  simply  the  '-eement  lemzth  (  Az)  and  AP is  the 
pressure  drop  in  one  single-phase  segment.; 


6)  Total  Pressure  Drop.  The  total  pressure  drop  in  the  evaporator  tube  is  the  sum  of  the 
pressure  drops  in  the  single- phase  and  two— phase  sections: 

APtfvup  =  AP  ,p  +  AP.,,  1  :-<«)! 

where  APcvap  represents  the  total  pressure  drop. 


Bounding  Calculations.  The  major  components  of  the  total  pressure  drop  are  friction  and 
fluid  acceleration  in  the  two— phase  forced  convection  region.  (Figure  2.16  illustrates  this 
dominance.)  There  is  little  uncertainty  in  the  acceleration  pressure  drop,  which  depends  only 
on  the  difference  between  inlet  and  outlet  conditions.  The  primary  uncertainty  is  the  pressure 
drop  calculation  due  to  interfacial  friction  between  the  liquid  and  vapor  phases  in  the  annular 
heat  transfer  regime  of  two— phase  forced  convection. 
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The  design  map  for  friction  pressure  losses  in  Figure  2.17  uses  the  correlation  of  Wallis 
(1969)  for  the  interfacial  friction  factor.  This  is  probably  a  conservative  estimate  when  applied 
to  space  thermal  management  systems,  since  it  is  based  on  air— water  experiments  in  which  the 
ratio  of  liquid  to  vapor  phase  densities  is  very  large. 

A  lower  bound  estimate  of  the  frictional  pressure  losses  in  the  two-phase  forced 
convection  regime  can  be  made  using  alternate  models  for  the  interfacial  friction.  A  limiting 
case  is  when  the  interfacial  friction  factor  is  set  equal  to  that  for  vapor  flowing  in  a  smooth 
pipe.  Frictional  pressure  drops  are  roughly  a  factor  of  3  lower  than  those  calculated  using  the 
Wallis  correlation.  Alternate  interfacial  friction  models  are  not  presented  in  this  Design 
Manual. 


2.2.3  Detailed  Equations  for  Pressure  Losses 

This  section  presents  pressure  drop  equations  which  can  be  used  in  detailed  computer 
calculations.  The  primary  model  is  separated  flow  for  the  annular  regime  in  two— phase  forced 
convection.  Homogeneous  flow  is  assumed  for  the  less  important  nucleate  boiling  regimes. 
Further  information  is  available  in  Collier  (1986)  or  Wallis  (1969).  Kakac  (1988)  has  shown 
that  this  method  works  well  for  boiling  systems  on  Earth.  Table  2.2  summarizes  the  pressure 
drop  calculations  for  all  regimes  in  an  evaporator  tube.  Figure  2.14  shows  a  scheme  for 
calculating  the  pressure  losses  during  the  detailed  heat  transfer  computation.  Equations  arc 
provided  for  friction,  body  force  and  fluid  acceleration  pressure  drops  in  the  various  heat 
transfer  regimes  of  an  evaporator  tube.  The  following  sections  provide  more  detail,  and  the 
basis  for  the  dimensionless  design  maps  is  also  explained. 

Friction.  The  basic  equation  for  the  pressure  gradient  due  to  wall  friction  is: 

(dP/dz)  f  =  -f  (4/Dh)  (G2/2p  0  (2-61) 

in  which  f  is  a  friction  factor,  (4/Dh)  is  the  perimeter  to  area  ratio,  and  (G2/2pi)  is  the  kinetic 
energy  per  unit  volume  of  the  fluid.  For  two— phase  flow  this  expression  is  modified  as 
follows: 


(dP/dz)  f  =  -fw  i  02  ( 1-x)2  (4/Dh)  (G2/2p  0  (2-62) 

Now  the  friction  factor  and  kinetic  energy  terms  are  calculated  as  if  the  flow  were 
single— phase  liquid.  Two— phase  effects  are  lumped  into  ri2  which  is  the  two-phase  multiplier. 
The  value  of  O2  depends  on  the  fluid  properties,  the  quality  and  the  flow  regime  as  described 
below. 


Two-Phase  Farced  Convection  Regime.  This  heat  transfer  regime  is  characterized  by 
annular  flow  and  the  two— phase  multiplier  is  (Collier,  1986): 

f2  -  «-«> 

where  a  is  the  local  void  fraction.  Then  the  pressure  drop  due  to  friction  in  a  short  tube 
segment  of  length  Az  can  be  calculated  using  this  expression: 

APfc  =  n%jz(l-x)2(4/Dh)(G2/2p1)  Az  (2-64) 
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where  APfe  is  a  pressure  drop  defined  by: 


zrh 

APfe  =  -J  (dP/dz)  f  dz  (2—65) 

Z  j 


Thus,  a  positive  value  of  AP  represents  a  pressure  loss,  and  a  negative  value  implies  a  pressure 
recovery.  The  void  fraction  a  is  the  only  unknown.  The  void  fraction  depends  on  the  local 
quality,  the  liquid  to  vapor  density  ratio,  and  f;/fwg.  If  the  interfacial  shear  ratio  is  known,  the 
void  fraction  is  related  to  the  quality  by: 


(2—66) 


This  equation  can  be  solved  iteratively  for  the  void  fraction  if  an  expression  for  f;/fwg  is 
known.  Wallis  (1969)  recommends  the  following  expression  for  the  shear  ratio: 


fi/fwg  =  [1  +75(l~a)] 


(2—67) 


Other  Heat  Transfer  Regimes.  In  the  nucleate  boiling  heat  transfer  regimes,  a 
homogeneous  flow  model  is  recommended  to  estimate  the  pressure  drop.  This  model  assumes 
that  the  liquid  and  vapor  have  the  same  velocity  and  that  the  two— phase  fluid  properties  can  be 
estimated  by  a  quality-weighted  average  of  the  liquid  and  vapor  properties. 

McAdams  (1942)  recommends  this  value  for  the  two-phase  multiplier  based  upon  a 
quality— weighted  viscosity  for  the  two-phase  fluid: 


4>2  = 


1  +  X 

gi-il' 

1  +  X 

yj.  -il 

.Pg  J. 

.M-g  J. 

(2-68) 


The  segment  pressure  drop  is  then: 

APfc  =  fw  ,q>2  (1— x)2  (4/Dh)  {GVlpu  Az  (2-o9) 

Dimensionless  Formulation.  The  dimensionless  formulation  for  frictional  pressure  losses 
assumes  that  the  entire  two— phase  section  of  the  evaporator  can  be  treated  as  an  annular  flow. 
The  rate  of  pressure  decrease  with  fluid  quality,  dP/dx,  can  then  be  written  as: 


dP  _  dP  dz 
dx  dz  <3x 

=  -  rrw ( 1  -x)2  (4/Dh)  (G2/2p  0  & 

(2 — 70a) 

(2— 70b) 

=  -[T^o-^n/RcXGWpi) 

(2-70c) 
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Using  the  dimensionless  form  of  the  pressure  drop  from  Equation  2—54,  Equation  2-70c 
is  rearranged  and  integrated: 


APfe* 


4fi 


(l-x)2  dx 

(l-cO* 


(2-71) 


For  the  design  map,  assume  that  the  entire  two— phase  section  of  the  evaporator  is 
characterized  by  annular  flow.  Then  the  integral  can  be  evaluated  numerically  using  Equation 
2—66  to  relate  a  to  x  and  Equation  2—67  for  the  interfacial  shear  ratio.  Figure  2.17  plots  the 
results. 


Body  Forces.  The  basic  differential  relation  for  the  pressure  gradient  drop  due  to  body 
forces  is: 


(dP/dz)g  =  [apgr  (1-a)  pi  ]  a  cos0  (2-72) 

Evaluation  of  this  pressure  drop  in  an  incremental  pipe  segment  is  thus: 

APge  =  -  [  a  pg  +  (1-a)  pi  ]  a  cos0  Az  (2-73) 

(The  pressure  drop  is  positive  when  the  inlet  pressure  is  greater  than  the  outlet  pressure.)  The 
pressure  drop  depends  on  the  void  fraction,  a,  which  in  turn  depends  on  the  heat  transfer 
regime  as  discussed  below. 

Two-Phase  Forced  Convection  Regime.  In  this  regime  the  void  fraction  should  be 
evaluated  using  the  same  procedure  as  described  above  for  the  case  of  the  friction  pressure 
drop.  (See  Equations  2—66  and  2—67.) 

Oth°r  Heat  Transfer  Regimes.  In  nucleate  boiling  heat  transfer  the  flow  is  assumed  to  be 
homogeneous  and  the  void  fraction  is  a  simple  function  of  the  quality  and  phase  densities: 


"(l-xXpg/pi)  +  x 

frt  the  single— phase  liquid  regime  the  void  fraction  is  zero. 

Dimensionless  Formulation.  The  dimensionless  expression  for  the  pressure  drop  due  to 
body  forces  is  obtained  as  follows: 

(E  =  =  [apg  +  tl-aJpiJacosOjj^h  (2-75) 

=  Pl-a  gh.c°se[a(p  /pi>  +  (l-ot)] 

^ev 
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Using  Equation  2—55,  we  can  rearrange  and  integrate  Equation  2—75  to  obtain: 


AP  * 


[a(pg/p  0  +  (1— a)]  dx 


(2-76) 


The  dimensionless  map  for  the  pressure  loss  due  to  body  forces  (Figure  2.18)  is  obtained  by 
again  assuming  that  the  entire  evaporator  can  be  characterized  by  the  dominant  heat  transfer 
regime  (two— phase  forced  convection)  and  using  the  void  fraction/quality  relationship  for 
annular  flow  (Equations  2—66  and  2—67)  to  numerically  evaluate  the  integral  in  Equation 


Acceleration.  The  basic  equation  for  pressure  gradient  due  to  fluid  acceleration  is: 


(dP/dz)a  =  -G2 


d  [  X  2  1  ( 1  — X  )2  1_  ' 

2z  [  a  pg  (1-a)  p, 


Thus  the  pressure  drop  in  a  pipe  segment  due  to  fluid  acceleration  is: 


(2-77) 


APa  =  G2A 


'  x2  1  t  (1— x)2  1 
a  pg  (1-a)  p! 


(2-78) 


where  the  symbol  "A"  is  used  to  represent  the  value  of  the  bracketed  quantity  at  the  tube  exit 
minus  the  value  at  the  inlet.  (Once  again,  the  pressure  drop  is  positive  when  the  pressure 
gradient  is  negative.)  The  pressure  drop  is  independent  of  the  segment  length  and  depends 
only  on  the  fluid  quality  and  void  fraction  at  the  inlet  and  outlet  of  the  segment.  This  pressure 
drop  is  evaluated  using  the  same  void  fraction/quality  relations  described  above  (Equations 
2—66  and  2-67)  for  the  pressure  losses  due  to  friction  and  body  forces.  The  dimensionless 
form  is  obtained  by  dividing  Equation  2—78  by  (G2/2pi),  as  in  Equation  2—56,  to  obtain: 


AP  *  =  A 


X2  e_L  .  (1-x)2 
a  pg  (1-a) 


(2-79) 


Figure  2—19  plots  this  dimensionless  result.  The  curves  on  the  plot  can  be  generated  by 
subtracting  the  values  of  the  expression  inside  the  brackets  of  Equation  2—79  when  it  is 
evaluated  at  the  evaporator  inlet  and  outlet  conditions. 


2.2.4  Validation  with  Microgravity  Data 

Lee  (1988)  measured  pressure  drops  in  a  straight— tube,  forced— convection  evaporator 
under  simulated  microgravity  conditions.  Unfortunately,  the  pressure  drop  data  vary  greatly 
with  the  variable  acceleration  on  the  experimental  package.  Cuta  and  Krotiuk  (1988)  show 
that  code  calculations  (using  COBRA/TOAC)  are  poor  compared  with  the  data.  Analysis  of 
these  data  is  beyond  the  scope  of  this  Design  Manual. 
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3  CONDENSERS:  High  Vapor  Shear  Condensers 

Condensers  transfer  heat  between  two— phase  heat  transport  systems  and  spacecraft  heat 
rejection  systems  such  as  radiators.  The  condenser  allows  for  controlled  condensation  of  the 
two-phase  working  fluid  so  that  heat  can  be  efficiently  and  reliably  removed  from  the 
two— phase  loop.  Overall  design  goals  are  similar  to  those  of  the  evaporator.  The  condenser 
should  be  small  and  lightweight,  have  small  pressure  losses,  and  operate  regardless  of  the 
spacecraft's  acceleration.  Small  component  size  requires  high  heat  fluxes  during  condensation. 
Since  condensate  films  are  typically  the  main  heat  transfer  resistance,  condensers  require  thin 
films  or  large  temperature  differences.  For  spacecraft  thermal  management  systems,  thin  films 
and  small  temperature  differences  are  preferred  in  order  to  minimize  radiator  area. 

The  absence  of  large  body  forces  in  microgravity  makes  phase  control  an  important 
design  issue.  Gravity  is  often  used  on  Earth  to  drain  condensate  from  the  heat  transfer  surface, 
but  in  space  the  condensate  tends  to  accumulate  and  impede  heat  transfer.  A  number  of 
different  schemes  have  been  devised  for  condensate  control  with  two  main  approaches:  high 
vapor  shear  (friction)  forces  and  capillary  forces. 

Vapor  shear  is  a  simple  way  to  control  the  thickness  of  the  condensate  layer  and  is  the 
focus  of  this  Design  Manual.  In  a  high  vapor  shear  condenser,  the  vapor  moves  at  high  speed 
and  exerts  a  high  drag  force  on  the  condensate  film.  The  film  is  thus  propelled  at  high 
velocity  and  its  thickness  and  resistance  to  heat  transfer  are  reduced.  This  section  presents 
recommended  design  methods  for  high  vapor  shear  condensers  where  the  flow  geometry  is  a 
straight  tube  with  uniform  cross-section.  An  alternate  concept,  not  covered  in  this  Design 
Manual,  is  a  tapered  flow  channel  which  maintains  high  vapor  velocity  as  condensation  occurs 
and  the  vapor  mass  flow  rate  decreases.  Thin  films  arc  maintained  over  a  longer  distance  than 
for  channels  with  uniform  cross-section.  The  design  trade-off  for  the  tapered  geometry  is 
increased  heat  transfer  (and  condensation  rate)  at  the  expense  of  increased  pressure  drop. 
(Furthermore,  a  condenser  with  a  tapered  tube  operates  best  at  only  one  flow  condition.)  Other 
proposed  concepts  involve  suction  of  liquid  through  the  condenser  wall  to  maintain  thin  liquid 
films.  The  condenser  models  presented  here  cover  basic  phenomena  which  are  similar  in  these 
alternate  designs. 

The  capillary  approach  uses  the  surface  tension  of  the  liquid  to  separate  the  phases.  No 
models  for  this  type  of  condenser  are  addressed  in  this  design  manual. 

Control  of  non-condensible  gases  is  another  important  design  issue.  If  these  gases 
accumulate  within  the  condenser,  they  can  impede  heat  transfer.  Thus  condensers  must 
provide  for  continuous  removal  of  non— condensible  gases  from  heat  transfer  surfaces.  In  high 
vapor  shear  condensers,  the  high  velocity  vapor  continuously  purges  non-condensible  gases 
from  the  condenser  tube. 

This  section  reviews  the  important  phenomena  during  condensation  heat  transfer  and 
presents  two  types  of  design  tools  for  high  vapor  shear  condensers.  Design  maps  have  been 
developed  for  scoping  design,  and  equations  are  presented  for  detailed  computer  calculations 
of  condenser  performance. 

This  section  focuses  on  the  high  vapor  shear  approach  to  condenser  design.  Shear 
condensers  are  used  in  the  Boeing/S undstrand  and  OAO  thermal  management  systems  being 
developed.  This  section  specifically  addresses  condensation  inside  a  circular  tube  with  a 
constant  diameter.  Understanding  this  basic  case  is  also  the  first  step  in  designing  more 
complex  components,  such  as  the  tapered  condensers  with  decreasing  cross-sections  which 
maintain  high  vapor  velocity. 
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The  discussion  is  divided  into  two  major  sections:  heat  transfer  (Section  3.1)  and 
pressure  drop  (Section  3.2).  The  heat  transfer  section  instructs  the  designer  in  selecting  tube 
length,  temperature  differences,  mass  flux,  and  other  thermodynamic  design  variables.  The 
pressure  drop  section  shows  the  designer  how  to  estimate  pressure  losses  so  that  the  pumping 
power  can  be  calculated. 


3. 1  Heat  Transfer 

The  condenser  has  heat  transfer  regimes  which  are  analogous  to  those  in  the 
straight— tube  evaporator.  The  dominant  regimes  are  annular  condensation  and  single— phase 
liquid  cooling  as  described  in  Section  3.1.1.  Section  3.1.2  presents  design  maps  which  enable 
the  designer  to  quickly  estimate  the  required  tube  length  in  both  of  these  regimes.  Section 
3.1.3  presents  detailed  models  for  all  heat  transfer  regimes  which  can  be  used  in  computerized 
calculations  for  final  design  or  optimization.  Finally,  Section  3.1.4  reviews  the  available  data 
from  microgravity  condenser  experiments. 


3.1.1  Introduction 

In  general,  there  are  several  possible  boundary  conditions  for  a  condenser  tube.  This 
section  focuses  on  the  case  of  a  fixed  temperature  on  the  outer  wall  of  the  condenser.  The 
input  variables  for  design  in  this  case  are  illustrated  schematically  in  Figure  3.1.  They  are  the 
mass  flux,  tube  diameter,  inlet  and  outlet  fluid  conditions,  and  the  temperature  difference 
between  the  condensing  fluid  and  the  condenser  secondary  coolant.  This  is  the  most  common 
situation  for  condensers,  approximating  such  physical  situations  as  radiators  and  heat 
exchangers  with  a  high  heat  capacity  in  the  flow  on  the  secondary  (coolant)  side.  This  case  is 
of  particular  interest  for  space  thermal  management  systems  because  heat  pipes  (which  are 
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Figure  3.1.  DESIGN  INPUTS  FOR  A  HIGH  VAPOR  SHEAR  CONDENSER  TUBE 
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constant  temperature  devices)  may  be  used  to  reject  heat.  The  basic  differential  relations 
presented  in  Section  3.1.3  are  generally  applicable  regardless  of  the  actual  boundary  conditions 
and  can  be  used  as  the  building  blocks  of  customized  design  tools  for  any  specific  application. 


Heat  Transfer  Regimes.  As  vapor  becomes  liquid  within  the  condenser  tube,  the  void 
fraction  and  quality  decrease,  leading  to  a  succession  of  flow  and  heat  transfer  regimes.  Figure 
3.2  illustrates  the  general  sequence  of  flow  and  heat  transfer  regimes  in  a  high  vapor  shear 
condenser.  In  general  the  sequence  from  the  condenser  inlet  to  the  outlet  will  be: 

1.  Annular  condensation, 

2.  Slug  condensation, 

3.  Bubbly  condensation,  and 

4.  Single— phase  liquid  cooling. 

Most  of  the  condensation  heat  transfer  occurs  in  the  annular  regime,  so  high  vapor  shear 
condensers  can  be  characterized  as  operating  mainly  in  annular  flow.  The  lengths  of  the  slug 
and  bubbly  regimes  are  typically  very  short  in  relation  to  the  overall  tube  length.  The  length 
of  the  single— phase  liquid  section  depends  on  the  degree  of  outlet  subcooling  which  is  desired. 
High  subcooling  requires  a  long  single-phase  section. 

Annular  Condensation.  This  is  the  first  heat  transfer  regime  in  the  condenser, 
characterized  by  high  vapor  velocity,  high  void  fraction  and  high  flow  quality.  Vapor  flows  in 
a  cylindrical  core  and  liquid  condensate  flows  in  an  annular  film  on  the  tube  wall.  Heat 
transfer  is  limited  by  conduction  across  the  liquid  film,  which  is  typically  in  laminar  flow.  As 
more  vapor  is  condensed,  the  thickness  of  the  liquid  film  increases  which,  in  turn,  reduces  the 
heat  flux.  High  gas  velocity  results  in  high  vapor  shear  which  improves  heat  transfer  by 
reducing  the  film  thickness. 
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Figure  3.2.  HEAT  TRANSFER  REGIMES  IN  A  HIGH  VAPOR  SHEAR  CONDENSER 
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The  annular  condensation  regime  ends  when  there  is  insufficient  vapor  left  in  the  tube  to 
support  annular  flow.  Usually  this  occurs  for  void  fractions  of  about  0.75,  corresponding  to  a 
flow  quality  of  about  0.10.  Condensate  bridges  the  tube,  and  the  slug  condensation  regime 
begins. 

Slug  and  Bubbly  Regimes.  Slug  and  bubbly  regimes  are  low  quality  regimes  intermediate 
between  the  annular  condensation  and  single-phase  liquid  regimes.  Typically  these  flow 
regimes  are  not  present  over  very  much  of  the  tube  length  and  can  be  neglected  in  first— order 
scoping  calculations. 

In  the  slug  condensation  regime,  vapor  flows  in  Taylor  bubbles  in  the  center  of  the  tube, 
separated  by  liquid  slugs.  Heat  transfer  occurs  by  two  mechanisms:  conduction  across  the 
liquid  film  which  surrounds  the  bubbles  and  single-phase  convection  into  the  liquid  slugs. 
The  heat  flux  in  this  regime  can  be  quite  high  because  the  slugs  move  very  quickly  initially, 
leading  to  high  convective  heat  transfer  coefficients.  As  vapor  condenses  and  the  void  fraction 
decreases,  the  slugs  become  bigger  and  the  bubbles  smaller.  Eventually  the  large  Taylor 
bubbles  become  unstable  and  break  up  into  smaller  bubbles,  initiating  the  bubbly  condensation 
regime. 

In  the  bubbly  condensation  regime  the  vapor  is  present  as  small  bubbles  in  a  continuous 
liquid  phase.  The  void  fraction  and  quality  are  both  low.  Heat  is  transferred  to  liquid  at  the 
tube  walls,  then  to  the  bubbles  which  condense. 

One  concern  with  this  type  of  condenser  is  that  the  condensation  process  for  vapor 
bubbles  may  be  unsteady  in  the  slug  and  bubbly  flow  regimes.  The  concern  is  that  this  would 
lead  to  unsteady  behaviors  similar  to,  but  not  as  severe  as,  waterhammer  phenomena. 
Preliminary  analysis  and  experimental  data  indicate  that  this  concern  is  unfounded. 

Single-phase  Liquid.  After  all  the  vapor  has  condensed,  the  liquid  approaches  the 
temperature  of  the  secondary  coolant.  Heat  transfer  in  this  regime  is  simple  forced  convection 
to  single-phase  liquid.  The  degree  of  subcooling  desired  and  the  temperature  of  the  secondary 
coolant  determine  the  length  of  the  "subcooler"  portion  of  the  condenser. 


Typical  Behavior.  The  calculated  flow  quality  and  void  fraction  in  an  R— 1 14  condenser 

are  illustrated  in  Figure  3.3.  This  condenser  tube  is  6  mm  in  diameter,  and  2.8xl()--'  kg/s  of 

condensing  fluid  flows  through  at  a  saturation  pressure  of  435  kPa.  The  temperature 
difference  driving  condensation  is  40  K. 

The  figure  shows  the  dominance  of  annular  conde  .sation.  In  this  regime,  the  void 

fraction  remains  high  (over  80  percent)  because  the  liquid  is  confined  to  a  very'  thin  film  on 

the  tube  wall.  The  quality  decreases  steadily  due  to  condensation.  When  the  void  fraction 
falls  to  a  critical  value  (about  0.73  m  into  the  tube)  liquid  slugs  start  to  form.  The  quality  at 
this  point  is  roughly  15  percent  and  decreases  rapidly  due  to  high  heat  transfer  in  the  slug 
regime.  The  large  Taylor  bubbles  condense,  and  bubbly  flow  is  present  for  a  few  centimeters. 
Single— phase  liquid  fills  the  final  15  cm  of  the  tube.  Only  a  small  amount  of  subcooling 
(indicated  by  the  small,  negative  quality  at  the  tube  exit)  is  achieved  in  this  short  length. 

Figure  3.4  shows  the  predicted  heat  flux  profile  in  this  condenser.  The  total  amount  of 
heat  transfer  is  proportional  to  the  area  under  the  heat  flux  curve.  Approximately  85  percent 
of  the  heat  transfer  occurs  in  the  annular  regime,  while  the  remaining  15  percent  occurs  in  the 
slug,  bubbly,  and  single-phase  regimes  combined.  In  the  annular  regime,  the  heat  flux  is 
inversely  proportional  to  the  thickness  of  the  liquid  film.  Thus  the  heat  flux  is  very  high  at  the 
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FLOW  QUALITY  OR  VOID  FRACTION 


Figure  3.3.  VOID  FRACTION  AND  FLOW  QUALITY  IN  A  HIGH  VAPOR  SHEAR  CONDENSER 


Figure  3.4.  HEAT  FLUX  PROFILE  IN  A  HIGH  VAPOR  SHEAR  CONDENSER 
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tube  entrance,  where  the  film  is  quite  thin,  but  decreases  along  the  length  of  the  tube  as  more 
vapor  condenses  and  the  film  thickens.  The  heat  flux  increases  in  the  slug  regime  because  the 
high  velocity  liquid  slugs  have  high  convective  heat  transfer  coefficients.  The  heat  flux 
decreases  sharply  in  the  short  bubbly  regime  and  reaches  its  lowest  value  in  the  single-phase 
liquid. 


3.1.2  Design  Maps 

The  most  common  condenser  design  problem  is  to  determine  the  tube  length  which  is 
necessary  to  condense  and  cool  the  working  fluid.  This  section  presents  dimensionless  maps 
which  enable  the  designer  to  estimate  the  length  of  the  condenser  given  the  major  design 
parameters. 

The  quick  design  procedure  is  to: 

•  Check  that  the  simple  design  maps  will  apply  (Equation  3—1), 

•  Calculate  the  thermodynamic  ratio  for  the  condenser  (Equation  3— 2), 

•  Estimate  the  length  of  the  two— phase  region  of  the  condenser  (Figure  3.5 
or  3.7  and  Equation  3-4), 

•  Estimate  the  length  of  the  single— phase  region  of  the  condenser  (Figure 
3.6  and  Equation  3-8),  and 

•  Calculate  the  total  length  of  the  condenser  (Equation  3—9). 

To  perform  the  calculations  presented  in  this  section,  the  following 

design  parameters  are  needed: 


Fluid  Properties 

Pi  Liquid  phase  density  (kg/m3) 

pg  Vapor  phase  density  (kg/m3) 

p.  i  Liquid  phase  viscosity  (kg/m-s) 

pg  Vapor  phase  viscosity  (kg/m-s) 

ki  Liquid  phase  thermal  conductivity  (W/m-  K) 
hr„  Latent  heat  of  vaporization  (J/kg) 

a  =  Surface  tension  (N/n.y 


Condenser  Tube  Properties 

tw  Condenser  wall  thickness  (m) 

kw  Wall  thermal  conductivity  (W/m— K) 

Geometry 

D  Condenser  tube  diameter  (m) 


Flow  Rates 

G  Mass  flux  of  condensing  fluid  (kg/m2— s) 
x0  Inlet  quality  of  condensing  fluid  (-) 
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Temperatures 

Tsat  Saturation  temperature  of  condensing  fluid  (K) 

T2  Temperature  of  the  secondary  condenser  coolant  (K) 
Tsub  The  condensing  fluid's  outlet  temperature  (K) 


Check  Tube  Thermal  Resistance.  The  design  maps  presented  here  will  apply  if  the 
thermal  resistance  of  the  condenser  tube  itself  is  less  than  that  of  the  condensate  film: 

-5-  >  200-4^  (3-1) 

l  w  Kw 

If  this  condition  is  not  satisfied  then  the  detailed  thermal  models  presented  in  section  3.1.3 
should  be  used  to  design  the  condenser. 


Condenser  Thermodynamic  Ratio.  In  order  to  use  the  design  maps  presented  below,  the 
value  of  the  condenser  thermodynamic  ratio,  Rc  must  be  calculated.  The  thermodynamic  ratio 
is  defined  as: 


Rc 


_  4  k,  ATsat 

~  "  G  D  H7g  - 


4  Ja 
Ref  Prj 


(3-2) 


where 


Ja 

Ref 

Pr, 


=  Tsat-T2, 

=  the  Jakob  number  (cplATsa/hfg), 

=  the  Reynolds  number  based  only  on 
the  liquid  flow  (GD/pi),  and 
=  the  Prandtl  number  of  the  liquid 

(M-iCpi/kf). 


The  thermodynamic  ratio  is  proportional  to  the  ratio  of  the  condensing  mass  flux  to  the 
inlet  mass  flux.  The  design  maps  below  only  apply  for  values  of  Rc  less  than  a  critical  value: 


Rc  <  1(H 


Basically,  this  is  the  range  for  which  the  design  curves  are  not  themselves  functions  of  Rc.  For 
thermodynamic  ratios  greater  than  1(H,  the  detailed  design  modeis  presented  in  Section  3.1.3 
are  recommended. 


Dimensionless  design  maps  can  be  made  for  thermodynamic  ratios  larger  than  iO4, 
however  in  this  case  the  maps  themselves  become  a  function  of  Rc. 


Length  of  Two-Phase  Condenser  Region.  Figure  3.5  shows  the  dimensionless  condenser 
length  required  to  completely  condense  the  fluid  from  the  quality  at  the  inlet  of  the  condenser 
(i.e.,  from  x  =  x0  to  x  =  0).  The  procedure  to  calculate  the  two-phase  condenser  length  is: 

1)  Calculate  the  ratio  of  the  liquid  to  vapor  density,  pi  /pg. 

2)  Locate  the  curve  in  Figure  3.5  which  corresponds  to  this  density  ratio. 
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CONDENSER  LENGTH 


Figure  3.5.  DESIGN  MAP  FOR  THE  LENGTH  OF  THE  TWO-PHASE  (ANNULAR)  REGION  OF 
THE  HIGH  VAPOR  SHEAR  CONDENSER  (fj/fwg  =  [1  +  75  (1-a)]) 


S'!  Find  the  point  on  that  curve  which  corresponds  to  the  required  reduction  in  fluid 
quality.  Determine  the  value  of  the  point  on  the  y-axis  which  corresponds  to  this  point.  This 
is  the  dimensionless  length.  The  dimensionless  length  is  determined  by  evaluating  an  integral 
''see  Section  3.1.3),  and  is  called  the  dimensionless  integral  I. 

4)  Calculate  L^,,  the  two-phase  length  required  to  condense  the  fluid: 


u 


Using  Equation  3—2  in  Equation  3-4  shows  that  the  two— phase  condenser  length 
depends  on  the  key  design  variables  as  follows: 


L? 


G  D2  hf, 

"MT," 


(3-5) 


For  example,  the  condensing  length  can  be  halved  by  doubling  the  temperature  difference 
between  the  two— phase  and  secondary  coolants  or  by  decreasing  the  tube  diameter  by  30 
percent.  The  length  of  the  condenser  is  doubled  if  the  mass  flux  is  doubled. 


Length  of  Single-Phase  Heat  Transfer  Region.  The  length  of  the  single-phase  portion 
of  the  condenser  depends  upon  the  secondary  coolant  temperature  and  the  degree  of  exit 
subcooling  which  is  required.  ("Subcooling"  is  the  amount  by  which  the  fluid  exit-  temperature 
is  less  than  the  saturation  temperature.)  If  only  a  small  amount  of  subcooling  is  required,  then 
the  single-phase  portion  of  the  condenser  does  not  need  to  be  very  long.  However,  in 
two— phase  loops  using  mechanical  pumps  the  subcooling  may  need  to  be  relatively  high  in 
order  to  provide  the  pump  with  adequate  net  positive  suction  head  (NPSH).  In  these  cases  the 
single— phase  length  can  be  a  significant  fraction  of  the  total  condenser  length  because  the  heat 
transfer  in  this  regime  is  much  less  efficient  than  in  the  condensing  section.  (Actual  thermal 
management  systems  may  have  a  separate,  single-phase  "subcooler"  to  cool  the  condensate 
more  efficiently  than  the  condenser  tube.) 

The  length  of  the  single-phase  section  depends  upon  the  liquid  Reynolds  number,  the 
liquid  Prandtl  number,  the  secondary  coolant  temperature  and  the  liquid  outlet  temperature. 
The  procedure  to  calculate  the  single-phase  length  is: 

1)  Calculate  the  liquid  Prandtl  number,  Prj: 

Prs  =  (3-6) 

2)  Calculate  the  liquid  Reynolds  number,  Re  (. 

Ref  =  GD/n,  (3-7) 

3)  Find  the  curve  in  Figure  3.6  which  corresponds  to  the  liquid  Prandtl  number  and 
locate  the  point  on  the  curve  which  corresponds  to  the  liquid  Reynolds  number.  This 
dimensionless  length  is  called  Isp. 


LIQUID  REYNOLDS  NUMBER, R*f 

Figure  3.6.  DESIGN  MAP  FDR  THE  LENGTH  OF  THE  SINGLE-PHASE  REGION  OF  THE 
CONDENSER 


4)  Calculate  L^,  the  required  length  of  the  single— phase  portion  of  the  condenser: 


Lsp  —  I$p  D  In 


sat' 


Tsub~T2 


(3-8) 


The  dependence  of  the  single-phase  tube  length  on  the  design  variables  is  somewhat 
counter-intuitive  because  condenser  tubes  with  larger  Reynolds  numbers  require  longer  tubes 
to  achieve  the  same  subcooling.  This  is  particularly  true  in  the  laminar  flow  regime.  The 
reason  is  that  increased  flow  through  the  condenser  tube  acts  primarily  to  reduce  the  length  of 
time  that  each  element  of  fluid  can  be  cooled.  There  is  a  compensating  effect  due  to  increased 
heat  transfer  coefficients  for  higher  flow  rates  (in  the  turbulent  flow  regime),  but  this  is  not 
enough  to  reduce  the  required  tube  length. 

Total  Condenser  Length.  The  total  length  of  the  condenser,  L^,,  is  the  sum  of  the 
two-phase  and  single— phase  lengths: 


Ltot  —  Ljp  +  Lgp  (3  9) 

Uncertainty.  The  major  uncertainty  in  this  design  calculation  is  in  the  value  of  the 
"interfacial  friction  factor".  This  parameter  quantifies  the  drag  between  the  vapor  core  and  the 
liquid  film  in  the  annular  condensation  regime.  The  design  map  which  has  been  presented 
(Figure  3.5)  gives  a  best-estimate  of  the  condenser  length  using  Wallis’s  correlation  for 
interfacial  friction  (Wallis,  1969).  This  correlation  represents  data  from  horizontal  and  vertical 
tubes  at  earth  gravity. 

If  the  tube  length  is  a  critical  parameter,  an  upper  bound  can  be  obtained  using  the 
alternate  design  map  shown  in  Figure  3.7.  This  figure  shows  the  same  design  information  as 
the  best-estimate  design  map,  except  that  the  vapor/liquid  interface  is  assumed  to  be  perfectly 
smooth  in  the  alternate  case  (i.e.  the  interfacial  friction  factor  is  the  same  as  if  the  gas  were 
flowing  alone  in  a  smooth  pipe).  The  procedure  to  use  this  alternate  design  map  is  the  same  as 
for  the  best-estimate  map. 

Lengths  calculated  from  the  alternate  map  are  typically  twice  the  lengths  calculated 
using  the  best-estimate  map.  These  conservative  estimates  of  the  condensing  length  are 
greater  than  the  best-estimate  values  because  lower  intert'acial  friction  maximizes  the  liquid 
film  thickness  and  thus  reduces  the  rate  of  condensation.  The  alternate  design  map  uses  the 
minimum  possible  interfacial  friction  and  thus  yields  the  maximum  possible  condensing 
length. 

One  reason  for  considering  this  bounding  case  for  interfacial  friction  and  condensation 
heat  transfer  is  to  bound  the  range  of  uncertainty  due  to  the  effects  of  gravity.  Preliminary 
data  from  microgravity  experiments  (Chen  et  al.,  1988)  indicate  that  intert'acial  friction  may  be 
less  than  at  earth  gravity. 

Another  reason  for  considering  the  bounding  model  is  the  effect  of  gravity  on  the  liquid 
distribution  around  the  perimeter  of  the  tube  in  the  two— phase  region.  For  this  low  gravity 
application,  the  condensate  is  assumed  to  be  distributed  uniformly  around  the  perimeter  of  the 
tube.  On  earth,  or  under  acceleration,  the  liquid  may  be  distributed  asymmetrically.  A 
detailed  analysis  by  Keshock  and  Sadeghipour  (1981)  concludes  that  heat  transfer  may  be 
slightly  worse  under  microgravity  conditions  than  at  earth  gravity  due  to  the  non-uniform 
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Figure  3.7.  DESIGN  LIMIT  FOR  THE  LENGTH  OF  THE  TWO-PHASE  (ANNULAR) 
CONDENSER  REGION  (fj  /fwg  =  1.0) 

distribution  of  the  liquid  film  (because  the  average  film  thickness  will  be  slightly  thicker  under 
microgravity  conditions).  Microgravity  experiments  aboard  an  aircraft,  however,  led  Keshock 
(1974)  to  conclude  that  the  average  heat  transfer  rates  were  about  the  same  under  microgravity 
and  earth  gravity.  So  the  simple  assumption  of  a  uniform  liquid  film  may  be  sufficient,  based 
upon  preliminary  data. 

Laminar  flow  of  liquid  in  the  liquid  film  has  been  assumed  so  that  thermal  conduction 
controls  the  heat  transfer.  If  the  Reynolds  number  in  the  liquid  film  is  large  (Ref  = 
(1  —  x)GD/pi  >  4000),  then  turbulence  in  the  liquid  film  can  increase  the  heat  transfer  in  the 
liquid  film.  The  condition  with  a  turbulent  liquid  film  is  not  treated  in  this  Design  Manual, 
although  the  approach  to  include  the  effect  is  outlined  in  the  Section  3.1.3. 

Another  modelling  uncertainty  is  the  effect  of  liquid  entrainment  on  heat  transfer  in  the 
condenser.  Entrainment  of  liquid  will  reduce  the  thickness  of  the  liquid  film  and  improve  heat 
transfer.  (The  effect  on  the  pressure  drop  wiii  be  small.)  This  effect  is  not  accounted  for  in 
this  Design  Manauai. 


3.1.3  Equations  for  Heat  Transfer 

To  first  order,  a  high— vapor  shear  condenser  can  be  sized  assuming  only  one  or  two  heat 
transfer  regimes  (annular  condensation  and  single-phase  liquid,  as  shown  in  Section  3.1.2). 
For  final  design  or  optimization,  however,  the  detailed  behavior  of  the  condensing  fluid  in  all 
heat  transfer  regimes  should  be  considered. 


This  section  presents  detailed  design  relationships  for  each  heat  transfer  regime. 
Procedures  are  given  to  generate  design  plots  for  specific  cases  or  to  program  computerized 
calculations  for  high  vapor  shear  condensers.  The  usual  way  to  use  these  equations  is  to 
divide  the  evaporator  tube  into  a  number  of  segments  and  to  assume  that  the  fluid  properties 
are  constant  within  each  segment.  In  each  segment  there  are  three  basic  calculations: 

•  Calculate  the  heat  flux, 

•  Determine  the  conditions  at  the  segment  exit,  and 

•  Check  for  regime  transition. 

Table  3.1  summarizes  the  equations  for  these  calculations  in  each  of  the  heat  transfer 
regimes.  Figure  3.8  shows  schematically  the  sequence  for  performing  these  calculations.  The 
following  sections  provide  more  detail. 

About  one  hundred  segments  are  usually  sufficient.  The  quality  and  flow  rate  entering 
the  first  segment  are  known  from  the  condenser  inlet  conditions.  Exit  conditions  calculated  for 
one  segment  become  the  inlet  conditions  for  the  next.  The  following  sections  show  how  to 
calculate  the  fluid  state,  heat  flux  and  regime  transitions  in  each  of  the  heat  transfer  regimes 
experienced  by  the  fluid  as  it  travels  through  the  condenser:  annular,  slug  and  bubbly 
condensation  followed  by  single-phase  liquid  cooling. 


Annular  Condensation  Regime.  The  major  phenomena  which  govern  heat  and  mass 
transfer  in  this  regime  are  conduction  of  heat  across  the  liquid  film  and  interfacial  shear. 
Figure  3.9  shows  the  major  variables  used  to  analyze  this  regime.  A  similar  analysis  has  been 
developed  by  Chow  (1986),  and  Collier  (1986)  also  discusses  film  condensation  with  vapor 
shear. 


The  procedure  requires  three  main  inlet  conditions  for  each  segment  to  be  calculated: 
liquid  velocity  u  {,  vapor  velocity  ug,  and  film  thickness  5.  The  procedure  below  shows  how  to 
calculate  these  three  parameters  at  the  segment  exit  as  a  function  of  the  inlet  values.  The  exit 
values  from  one  segment  can  be  the  inlet  values  for  the  next.  Or  a  more  complex  iterative 
scheme  is  possible  in  which  the  inlet  and  exit  values  of  velocity  and  film  thickness  for  each 
segment  are  used  to  estimate  the  average  values  in  that  segment.  The  procedure  which  follows 
can  be  adapted  for  either  method. 

Inlet  Conditions.  To  begin  calculations  in  the  annular  regime,  the  following  parameters 
must  be  known: 


•  Fluid  properties  —  pi,  pg,  Pi,  p„,  hfg,  kh  a, 

•  Temperatures  —  Tsal,  T-2, 

•  Flow  rates  -  G,  x0  (inlet  quality), 

•  Velocities  -  ug0,  Uf0  (gas  and  vapor  inlet  velocities),  and 

•  Inlet  film  thickness  —  o0. 
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TRANSFER  CAI£U1-ATI0NS  FOR  A  HIGH  VAPOR  SHEAR  CONDENSER 
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Figure  3.8.  FLOWCHART  FOR  DETAILED  CONDENSER  CALCULATIONS 
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Figure  3.9.  ANALYSIS  VARIABLES  FDR  ROW  IN  THE  ANNULAR  CONDENSATION  REGIME 


The  fluid  properties,  temperatures  and  flow  rates  should  all  be  known  independently  of 
the  condenser  design.  The  inlet  velocities  and  film  thickness  must  be  calculated.  The 
procedure  is: 


1) 


2) 


3) 


Calculate  ugo,  the  gas  inlet  velocity,  assuming  homogeneous  flow: 


i  x. 


u„0  =  G  i  — 


I 


iPs  Pi  . 

Calculate  the  gas/wall  friction  factor,  fwg: 

rpgug0Dj-°-“ 


fW2  =  0.046 


Pg 


Solve  for  50,  the  inlet  film  thickness: 


Sg 


1  +300 


50 

D 


GD(l-x0)p, 

7 

f  w  gPlPgugo 


(3—1  O') 


(3-11) 


(3-12) 
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(This  equation  equates  the  shear  stress  at  the  film  surface  with  the  interfacial  shear  using 
Wallis's  (1969)  correlation  for  f  i  /fwg-) 

4)  Calculate  the  inlet  liquid  velocity,  Uf0: 

GD(l-Xo) 

Ufo  = -  (3-13) 

4p]50 


Heat  Flux.  The  heat  flux  in  each  segment  depends  on  the  film  thickness  and  possibly 
the  condenser  tube  thickness: 

(3-14) 


In  many  cases  it  is  reasonable  to  neglect  the  tube  wall  thermal  resistance  (tw/kw)  compared  to 
that  of  the  liquid  film  (S/kj),  especially  if  the  tube  wall  is  metal.  Since  the  flow  is  often 
observed  through  a  glass  tube  in  experiments,  the  wall  resistance  could  be  significant  in  this 
case  and  should  then  be  included. 

For  turbulent  liquid  films,  the  term  8/ki  for  thermal  conduction  in  the  liquid  film  should 
be  replaced  or  modified  by  1/h,  where  h  is  a  convective  heat  transfer  coefficient  for 
condensation  on  liquid  films.  Collier  (1986)  presents  a  correlation  by  Soliman  et  al.  which  has 
been  found  to  agree  well  with  experimental  results. 

A  word  about  notation:  In  the  stepwise  process  of  analyzing  successive  segments  in  the 
condenser,  the  inlet  quantities  (subscript  o)  are  used  to  determine  the  values  at  the  exit  of  the 
first  segment  (subscript  k  for  the  current  segment,  where  k  =  1  for  the  first  segment).  From 
that  point  on,  values  at  the  exit  of  the  preceding  segment  (subscript  k— 1)  become  the  inlet 
values  for  the  current  segment  being  analyzed. 


Fluid  State.  In  each  segment  the  quality  and  film  thickness  lor  void  fraction)  at  the 
segment  exit  must  be  calculated  as  functions  of  the  inlet  parameters  and  the  segment  heat  flux. 
The  procedure  is: 

1)  Calculate  Gc,  the  condensation  mass  flux: 

Gc  =  q"/hfg  (3-15) 

2)  Calculate  xv  the  interfacial  shear,  using  this  expression: 

Pg(Ug-Uf)2 

Xi=  fwg  (fi/fwg) - +  Gc(ug— Uf)  (3-16) 

2 


’pgUgDl-0-2 

where  fW2  =  0.046  -  (3-17) 

L  Pg  . 
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(3-18) 


fi/fwg  =  [1  +  300(8/D)] 

and  8  is  the  inlet  (or  average,  see  below)  film  thickness  in  the  segment. 

3)  Calculate  the  film  thickness  at  the  segment  exit.  Assuming  constant  interfacial 
shear  in  the  segment,  the  film  thickness  at  the  exit  can  be  calculated  using  this  equation: 

r  r  3Azkp.2  Ja  1]  1/3 

5k  =  (6,-03  +  - 1 -  (3-19) 

[  bik  Pi  PriJ. 

where  the  subscript  "k"  refers  to  the  segment  exit  and 

t  j,  =  interfacial  shear  in  the  segment  (N/m2), 

Ja  CpjAT sat  /h  fg,  and 

Pri  =  4icpi/ki- 

4)  Calculate  the  average  velocity  in  the  liquid  film  at  the  segment  exit: 

8k — ]  ^3  cAzk 

Uflc  =  Uflc-j  -  + -  (3-20) 

.  8k  J  p [8k 

where  u^  is  the  film  velocity  at  the  exit  of  the  previous  segment,  and  Azk  is  the  length  of 
segment  "k". 

5)  Calculate  the  vapor  velocity  at  the  exit  of  the  segment: 

4GcAzk 

Ugk  =  ugk_! -  (3-21) 

PgD 

Initially  this  procedure  can  be  used  in  a  "marching  scheme",  in  which  the  interfacial 
shear  and  heat  flux  are  calculated  using  the  inlet  conditions  (film  thickness  and  phase 
velocities)  applied  over  the  entire  length  of  the  segment.  However,  once  the  entire  condenser 
has  been  calculated  in  this  way  the  designer  will  have  estimates  for  the  film  thickness  at  the 
inlet  and  exit  of  each  segment.  Averages  for  these  parameters  in  the  segments  can  oe 
estimated  by  taking  the  mean  of  the  inlet  and  exit  values,  then  the  interfacial  shear  and  the 
heat  flux  can  be  recalculated  based  on  segment  average  values  (of  5,  Uf  and  ug).  This 
procedure  can  be  repeated  until  the  results  converge  to  any  desired  degree.  This  iterative 
procedure  leads  to  more  consistent  results  which  are  less  sensitive  to  the  size  of  the  segments. 

Dimensionless  Design  Map.  All  the  heat  transfer  phenomena  in  the  annular  condensation 
regime  can  be  summarized  in  a  single  design  map,  as  shown  in  Figure  3.5.  The  curves  on  this 
map  are  derived  from  a  dimensionless  energy  balance  in  which  conduction  heat  transfer 
through  the  condensate  film  balances  the  energy  removed  by  the  vapor  during  condensation. 

Figure  3.10  shows  a  control  volume  of  length  dz.  Equating  the  change  in  fluid  enthalpy 
with  the  heat  flux  across  the  tube  walls: 

G  (tcD2/4)  hfgdx  =  -q"(7tD)dz  =  ~kl  (TtD)dz  (3-22) 
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Figure  3.10.  CONTROL  VOLUME  FOR  HEAT  TRANSFER  ANALYSIS 
Rearranging  Equation  3—22  and  using  the  approximation  that  5  =  (l-a)D/4  gives: 

(l-oc)dx  =  ~16gkjj)2AhSfa!  ~dz  =  —  (dRc/D)  dz  =  -4  Rc  dz"  (3-23) 

(where  z*  =  z/D).  An  important  result  which  we  will  use  later  to  calculate  pressure  losses  is: 

dz  (1-cOD 
ax  = 

Integrating  Equation  3—24  from  quality  xo  at  z  =  0  to  quaiitv  x  at  z  =  L  gives: 


4  Rc  L*  =  -  ■  i  1  — ot)  dx  =  '  i  1  — <jc;  vJx  3-25 1 

The  final  quality,  xs,  can  be  taken  as  zero  tor  a  first  approximation  or  set  to  0. i, 
corresponding  to  the  transition  from  annular  to  slug  condensation.  To  evaluate  the  integral  a 
relationship  between  the  flow  quality,  x,  and  the  void  fraction,  a,  is  necessary.  An  appropriate 
model  has  already  been  described  for  the  annular  flow  regime  in  Equation  2^66: 


0  i  (1— 0t)2 

Ps  a-5 


(fi/fwg) 


(3-26) 
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The  integral  in  Equation  3—25  is  evaluated  numerically  using  Equation  3—26  for  various 
constant  values  of  (pi/pg).  The  results  are  presented  in  the  design  map  of  Figure  3.5. 

A  simple,  closed-form  solution  can  be  derived  if  we  assume  that  f,/fWg  is  constant  and 
that  a1-25  «  1.0  over  the  range  of  the  integral  in  Equation  3—25.  The  closed-form  expression 
is: 


I  =  4RCL* 


1 

4(pi  /pg)*  (fj  /fwg^ 


where  xQ 

xs 


=  flow  quality  at  the  condenser  inlet,  and 
=  0.1,  corresponding  to  the  slug  regime  transition. 


(3-27) 


If  we  take  f;  /fwg  to  equal  approximately  10.0,  then  the  results  of  this  simple  equation  are 
within  20%  of  the  design  map. 


Transition  to  t he  slug  regime.  A  convenient  rule  of  thumb  for  microgravity  condensers 
is  to  locate  the  end  of  the  annular  condensation  region  using  a  critical  void  fraction  of  0.80. 
When  the  void  fraction  decreases  below  this  value,  assume  that  the  slug  condensation  regime 
has  begun. 

This  simple  rule  is  an  approximation  to  the  more  precise  method  for  calculating  the 
regime  transition.  The  end  of  the  annular  regime  can  be  found  using  the  same  methods 
described  in  Section  1.1  of  this  Design  Manual.  At  the  exit  of  each  segment,  calculate  the 
Martinelli  parameter,  X.  Figure  1.5  shows  the  value  of  X  which  corresponds  to  the 
annular/slug  transition  as  a  function  of  the  dimensionless  slope  parameter,  Y.  When  the 
transition  value  of  X  has  been  exceeded,  the  condenser  has  entered  the  slug  condensation 
regime. 


Slug  Flow  Regime.  Heat  transfer  in  the  slug  flow  regime  is  calculated  as  a  weighted 
average  of  film  conduction  into  the  vapor  bubbles  and  single— phase  convection  into  the  liquid 
slugs.  Figure  3.11  illustrates  the  analysis  of  this  regime.  The  segment  of  interest  is  denoted 
by  the  subscript  "k.”  Entrance  conditions  are  denoted  by  the  subscript  "k— 1"  and  exit 
conditions  have  the  subscript  k. ' 

The  calculation  procedure  for  the  slug  regime  requires  four  inlet  conditions  for  each 
segment.  These  are  the  quaiitv,  xk_:,  void  fraction.  total  volumetric  flux.  jk_;.  and  a 
parameter  equal  to  the  fraction  of  the  pipe  cross-section  which  is  occupied  by  the  liquid  trim 
which  surrounds  the  Taylor  bubbles.  mk_t.  The  heat  flux  in  a  segment  is  calculated  based  on 
these  inlet  conditions.  The  heat  flux  implies  a  certain  rate  of  condensation  in  the  segment 
which  enables  the  designer  to  calculate  the  segment  exit  conditions. 

Inlet  Conditions.  At  the  end  of  the  annular  regime  the  tlow  quality,  x0,  is  known.  The 
designer  must  calculate  the  other  three  inlet  conditions  for  the  initial  slug  segment  (a0,  j0  and 
m0)  using  the  following  procedure. 

1)  The  total  volumetric  flux,  jOI  is: 


1-Xol 
Pi  j 


1 3-28) 
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Figure  3.11.  ANALYSIS  OF  SLUG  CONDENSATION 


To  calculate  m0,  first  calculate  the  bubble  velocity  ubo  (Wallis,  1969): 


Vj 


O  ! 


0 


Toe  ..rea  ratio  can  then  be  read  from  r-igure  or  calcaiatea  by: 

m.Uko  •.).:->?'•  - 

;r, ,  nun. 0.1 6.  0.56<  !— e.;p  -2.od - 


Calculate  the  initial  void  fraction,  a0  (Wallis,  1969): 


x0(l-m0) 


Heat  Flux.  The  heat  flux  in  segment  "k"  is  calculated  based  on  an  average  heat  transfer 
coefficient  which  accounts  for  conduction  from  Taylor  bubbles  across  the  liquid  film  and 
forced  convection  from  the  liquid  slugs.  The  procedure  to  calculate  the  heat  flux  is  as  follows: 


1)  Calculate  the  bubble  length  ratio,  Lb  /Lt.  Lt  is  the  length  of  a  Taylor  bubble  plus 
its  associated  liquid  slug,  and  Lb  is  the  bubble  length.  The  ratio  is: 


a 


1— m 


(3—32) 


2)  Calculate  the  forced  convection  heat  transfer  coefficient.  hfc: 


r  fk 


I  (0.023)!  — 

ID 


P  ijD] 


o.s 


Pr  [0  •  33 


hfc  =  \ 

I 


L  Pi  J 


(3.66) (k  ;/D; 


P  ijD 

for - >  1500  (tu  r  bulent)  (3-33ai 

Pi 

P  IJD 

for - <  loOO  (.laminar)  :3-33b) 

Pi 


3)  Calculate  hf,  the  heat  transfer  coefficient  representing  conduction  across  the 
liquid  film: 


_ 2k  i 

hf  =  r  0.5'.  (3—34) 

|D[1  -  (l-mk)  I 


Figure  3.12.  VELOCITY  AS  A  FUNCTION  OF  jn,  /a  (Suo  and  Griffith,  1963) 
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4)  Calculate  the  average  heat  transfer  coefficient,  hav: 


f  Lb 

r  Lb  I 

hf  + 

i 

.  Lt  j 

.  Lt  . 

5)  Calculate  the  heat  flux: 

ATsa  t 

q"  =  - : 

1  i  ^  w 

R —  +  IT 

|_u  a  v  Kw. 


(3-35) 


(3-36) 


Fluid  State.  The  heat  flux  from  the  fluid  results  in  decreasing  quality  along  the  tube 
length.  The  rate  of  condensation  controls  the  segment  outlet  conditions.  To  calculate  the 
outlet  condition,  use  the  following  procedure: 


1)  Calculate  the  flow  quality  at  the  outlet  of  segment  "k:" 


*k  =  *k-i 


4  Azk 


2) 


G  D  hfg 

where  Azk  =  the  length  of  segment  "k". 

Calculate  jk,  the  total  volumetric  flux  at  the  exit  of  the  segment: 


(3-37) 


jk  =  G 


xk  l-xk 
—  + - 


Pg  P  i 


(3-38) 


3)  Calculate  mk,  the  film  area  ratio  at  the  segment  exit.  First  calculate  the  bubble 
velocity  at  the  exit  of  the  segment,  ubkt 


r  r  r  fpiiO0-8^ 

Ubk  =  jk  V  +  1-27'  i  —  exp  1  —3.8 i -  *3-39) 

i  l  i  i  o  j  jjj 

Use  this  value  of  the  bubble  velocity  to  calculate  mk: 

|  !  !  M- 1 u  v>k  1 0  ‘  ~ 6 ' 1  ;  I 

mk  =  min  ;0.16,  0.56-f  1  —  exp  i  — 2.64  I - U  (3-40) 

[  ■  i_  (  a  _  j 

xk(l-mk) 

4)  Calculate  ak  =  -  (3-41) 

Xk+(Pg  /Pi)0-*k) 

These  steps  (Equations  3-37  through  3^11)  can  be  used  initially  in  a  marching  method: 
the  film  thickness,  bubble  length  ratio  and  liquid  heat  transfer  coefficient  are  calculated  from 
the  inlet  parameters  only  and  assumed  to  apply  over  the  entire  length  of  the  segment.  A  more 
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sophisticated  scheme  will  use  this  marching  calculation  as  the  first  step  in  an  iterative 
procedure.  Once  the  exit  parameters  for  a  segment  have  been  estimated  from  the  marching 
scheme,  more  precise  estimates  of  the  film  thickness,  bubble  length  ratio  and  liquid  heat 
transfer  coefficient  can  be  made  based  on  the  segment  average  values  of  the  quality  and  void 
fraction.  The  refined  film  thickness,  bubble  length  ratio  and  heat  transfer  coefficient  will,  in 
turn,  result  in  better  estimates  of  the  exit  void  fraction  and  quality.  This  procedure  can  be 
continued  until  the  exit  conditions  converge,  and  the  results  will  be  less  dependent  or.  the 
length  of  the  segments. 


Transition  to  the  bubbly  regime.  As  more  vapor  condenses  the  Taylor  bubbles  will 
become  shorter  and  shorter,  and  the  volumetric  flux  will  decrease.  Slug  condensation  ends 
when  the  bubbles  become  unstable  and  break  up  into  smaller  bubbles. 

This  transition  can  be  calculated  using  the  procedure  described  in  Section  1.1.3.  A 
simple  rule  of  thumb,  accurate  at  acceleration  levels  below  10~2g,  is  to  assume  that  bubbly 
flow  begins  once  the  void  fraction  drops  below  0.45. 


Bubbly  Flow  Regime.  In  this  regime  the  vapor  phase  consists  of  small  bubbles  dispersed 
through  the  liquid.  A  homogeneous  model  is  recommended  to  analyze  this  section  of  the 
condenser.  In  this  case  only  the  flow  quality  at  the  segment  inlet  is  required.  The  following 
procedure  shows  how  to  calculate  the  heat  flux  and  exit  quality. 

Inlet  Conditions.  The  flow  quality  at  the  inlet  to  the  bubbly  condensation  region  is 
simply  the  exit  quality  from  the  last  segment  of  the  slug  region. 

Heat  Flux.  The  heat  flux  is  calculated  using  standard  methods  for  forced  convection. 
The  effect  of  two-phase  flow  is  to  increase  the  effective  Reynolds  number.  The  heat  flux  is 
estimated  using  this  procedure: 

1)  Calculate  jl,  the  average  viscosity  of  the  bubbly  mixture: 

r x  i— 

U  =  1 —  + - ! 

lUg  Hi  j 

2)  Calculate  the  effective  two -phase  Reynolds  number.  Re.p 

G  D 

Retp  =  ~ r~ 

'JL 

3)  Calculate  the  two— phase  heat  transfer  coefficient,  h^.  The  heat  transfer 
coefficient  depends  on  the  two-phase  Reynolds  number  as  follows: 


(3-42) 


i  j — o  i 
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h  tp  — 

4) 


(0.023) 


fk,] 


D 


0.8 


(Re  tp)  (PrO0-33  Re,p> 


(3.66)  (k  i/D) 


Re  tp  < 


Calculate  the  heat  flux  in  the  segment: 


1500  (tu  r  bulent) 

1500  (laminar) 


AT  s  a  t 


'  1  ,  tw 

K -  +  \T 

II  (p  Kw 


(3-44a) 

(3-44b) 


(3-45) 


Fluid  State.  Calculate  the  quality  at  the  exit  of  the  segment  Uaing  this  equation: 


4  q  ’’  AZfc 

xk  =  xk_; -  (3-46) 

G  D  h  fg 

As  in  the  previous  flow  regimes,  this  method  can  first  be  applied  in  a  marching  method 
in  which  the  entire  segment  is  characterized  by  the  inlet  quality.  Subsequent  iterations  use  the 
average  qualify  in  the  segment  to  calculate  the  heat  flux. 

Transition  to  single-phase  regime.  The  bubbly  condensation  regime  ends  when  the  void 
fraction  and  flow  quality  are  zero.  From  this  point  to  the  condenser  exit  there  is  only 
single-phase  liquid. 


Single— Phase  Liquid  Regime.  Single-phase  flow  is  considerably  simpler  than 
two— phase  flow  because  the  fluid  properties  are  constant,  and  heat  transfer  can  be  analyzed 
using  standard  methods.  There  is  no  need  to  divide  the  single-phase  region  into  segments 
because  the  entire  single-phase  section  can  be  modeled  analytically.  The  single-phase  regime 
is  the  final  heat  transfer  regime  in  the  condenser. 

Inlet  Conditions.  The  temperature  at  the  inlet  to  the  single-phase  liquid  region  is  the 
saturation  temperature.  Tsa[. 

Heat  Flux.  The  heat  flux  in  the  single— phase  region  depends  on  the  fluid  properties,  the 
flow  rate,  and  the  secondary  cooiant  temperature.  The  heat  flux  can  be  calculated  as  a 
function  of  position  with  this  procedure: 

1)  Caicuiate  Ref,  which  is  the  Reynolds  number  in  the  single-phase  region: 

Ref  =  GD/4,  (3-47) 

2)  Calculate  the  fluid  heat  transfer  coefficient,  hfc.  The  equation  for  the  heat  transfer 
coefficient  depends  on  whether  the  flow  is  laminar  or  turbulent,  as  indicated  by  the  Reynolds 
number. 


Ill 


il  fc  — 


k  1  0.8 

(0.023)  —  Rcf  Pr^-33 
D 

k  i 

(3.66)  — 

D 


for  Ref  >  1500  (turbulent) 

for  Ref  <  1500  (laminar) 


(3-48a) 

(3^48b) 


3)  Calculate  U,  the  overall  heat  transfer  coefficient  which  accounts  for  the  thermal 
resistance  of  the  fluid  boundary  layer  and  the  tube  wall: 


u  =  ffrr  +  ir 

L“  fc 


,-l 


If  the  wall  thermal  resistance  is  negligible  then  U  =  hfc. 


(3-49) 


4)  The  heat  flux  as  a  function  of  position  in  the  single— phase  section  of  the  condenser 
tube  can  be  calculated  using  this  equation: 


q"(zsP) 


U  ATsatexp 


r-4  u  z 
G  Dc 


pi. 


(3-50) 


where  zsp  is  the  position  in  the  tube  measured  from  the  start  of  the  single— phase  liquid  region. 


The  total  heat  transferred  in  the  single-phase  region,  qsp,  can  be  calculated  by: 


flsp  =  Gcpi 


1  -  exp 


[-4  U  L 
G~T)~c 


(3-51) 


where  Lsp  is  the  total  length  of  the  single-phase  region. 

Fluid  State.  The  fluid  temperature  as  a  function  of  position  in  the  single-phase  region  is 
given  by  this  equation: 


T  i"  \zsp,l 


|  ,  r— 4  U  ' 

i  1  ~exp[  CJ  D  cpu  j 


•3-52) 


The  fluid  temperature  at  the  condenser  exit  is  simply  TTLsp). 


Dimensionless  Man.  A  simple  dimensionless  design  map  can  be  constructed  for  the 
common  situation  in  thermal  management  applications,  in  which  the  thermal  resistance  of  the 
tube  wall  is  negligible  in  comparison  with  the  fluid  boundary  layer.  In  this  case,  the  overall 
heat  transfer  coefficient,  U,  is  approximately  the  fluid  heat  transfer  coefficient,  h  fc.  The  fluid 
temperature  in  the  single-phase  region  is  then  given  by: 


Tf(Zsp)  =  Tsat  -  ATS 


4- 


exp 


4  zsp 
G  D  ell 


=  "^sat  sat  •  1  exp 


rlNufc  is* 
Ref  Pr i  D 


(3-53) 

(3-54) 


where  Nufc  is  the  Nusselt  number  for  forced  convection  =  hfcD/k]. 
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This  equation  can  be  rearranged  into  a  form  suitable  to  produce  the  universal  design 

map: 


I 


sp 


1 


Re  f  Pr  i 

TTferr 


(3-55) 


where  is  used  to  represent  the  length  of  the  single-phase  section  with  an  outlet 
temperature  of  T^. 


The  Nusselt  number  is  actually  a  function  of  the  Reynolds  and  Prandtl  numbers, 
depending  on  whether  the  flow  is  laminar  or  turbulent: 


f  (0.023)Re/>-8Pr10-33  (Re f>  2300)  (turbulent) 


No*  = 


[3.66 


(Ref  <  2300)  (laminar) 


(3-56a) 

(3-56b) 


Equations  3.55  and  3.56  are  used  to  generate  the  design  map  (Figure  3.6)  for  the  length  of  the 
single-phase  section  of  the  condenser. 


3.1.4  Validation  with  Microgravity  Data 

No  data  are  presently  available  to  validate  these  models  for  condenser  performance  in 
microgravity.  The  key  issue  to  be  resolved  is  the  possible  effect  of  gravity  on  interfacial 
shear.  Previous  experiments  by  Keshock  (1974)  investigated  a  high  vapor  shear  condenser 
under  microgravity  (aircraft)  and  earth  gravity  conditions.  Insufficient  data  are  available  from 
these  tests  for  quantitative  comparison  with  the  design  models.  Qualitatively,  Keshock 
observed  that  the  heat  transfer  behavior  was  not  significandy  different  at  both  gravity  levels. 
We  have  reviewed  the  data  of  Lee  (1987)  who  operated  a  high— vapor  shear  condenser  in 
simulated  low  gravity  aboard  a  KC-135  aircraft.  The  condenser  did  no;  reach  steady  state  in 
these  tests,  and  consequently  the  data  do  not  appiv  to  the  models  developed  here.  Transient 
models  are  needed  to  predict  these  data.  It  is  expected  that  new  data  to  be  obtaineu  in 
experiments  sponsored  by  the  Weapons  Laboratory  will  provide  data  needed  for  quantitative 
analysis. 


3.2  Pressure  Losses 
3.2.1  Introduction 

Pressure  losses  in  the  condenser  occur  primarily  in  the  annular  condensation  regime. 
Just  as  in  the  case  of  the  forced  convection  evaporator,  there  are  three  major  components  to 
the  pressure  drop: 

•  Friction, 

•  Body  forces,  and 

•  Fluid  acceleration. 
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The  total  pressure  drop  is  due  to  the  sum  of  these  component  pressure  drops.  Because  high 
vapor  shear  leads  to  improved  heat  transfer,  the  frictional  component  should  be  dominant. 
Note  that  for  the  case  of  a  condenser  the  pressure  drop  due  to  fluid  acceleration  is  negative. 
Deceleration  of  the  fluid  from  the  vapor  to  liquid  phase  acts  to  reduce  the  overall  loss  in 
pressure.  If  the  frictional  losses  do  not  dominate  the  pressure  recovery  due  to  acceleration 
effects,  then  the  condenser  may  cause  instabilities  in  the  two-phase  system.  So  a  design 
criterion  for  this  type  of  condenser  is  that  frictional  pressure  losses  should  dominate. 

Pressure  losses  in  an  R— 114  condenser  tube  are  illustrated  in  Figure  3.13.  The  frictional 
and  fluid  acceleration  pressure  drops  are  shown  as  a  function  of  tube  length  in  all  heat  transfer 
regimes.  Note  that  nearly  all  of  the  pressure  losses  occur  in  the  annular  regime  where  the 
velocities  are  high.  The  totai  pressure  drop  is  only  about  half  that  due  to  friction  alone 
because  of  the  pressure  recovery  due  to  fluid  deceleration.  (The  acceleration  level  for  this 
calculation  is  lCh3g  and  the  pressure  losses  due  to  body  forces  are  negligible). 


3.2.2  Design  Maps 

Dimensionless  design  maps  are  provided  for  each  component  of  the  pressure  drop  in  the 
two— phase  annular  condensation  regime.  The  slug  and  bubbly  condensation  regimes 
contribute  only  small  amounts  to  the  total  pressure  drop  and  are  addressed  in  detail  only  in 
Section  3.2.3.  Pressure  losses  in  the  single— phase  region  are  calculated  by  standard  methods. 
The  relevant  equations  are  included  here. 

The  overall  procedure  to  estimate  the  condenser  pressure  losses  is: 

1.  Estimate  the  frictional  pressure  loss  in  the  two— phase  region  of  the  condenser 
(Figure  3.14  and  Equation  3—60), 

2.  Estimate  the  pressure  loss  due  to  body  forces  in  the  two— phase  region  of  the 
condenser  (Figure  3.15  and  Equation  3-61), 

3.  Estimate  the  pressure  loss  due  to  fluid  acceleration  (Figure  3.16  and  Equation 
3-62), 

4.  Add  these  components  together  to  obtain  AP.p,  the  total  pressure  drop  in  :he 
two-phase  section  of  the  condenser  (Equation  3-03), 

5.  Calculate  APsp,  the  pressure  drop  in  the  single— phase  section  (Equation  3—64), 
and 

A.  Add  the  pressure  drops  in  the  two— phase  and  single— phase  sections  to  estimate 
the  total  pressure  drop  (Equation  3—65). 

The  following  sections  detail  these  calculations. 

These  design  maps  presented  below  are  valid  only  for  a  certain  set  of  conditions  related 
to  the  thermal  resistance  of  the  condenser  wall  and  the  condenser  thermodynamic  ratio.  Refer 
to  Section  3.1.2  for  procedures  to  check  the  validity  of  the  dimensionless  design  maps  for  a 
specific  application. 
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Figure  3.13.  PRESSURE  LOSSES  IN  A  HIGH  VAPOR  SHEAR  CONDENSER 


The  design  parameters  needed  to  estimate  the  pressure  drop  in  the  condenser  include 

Fluid  Properties 

p  i  liquid  density  (kg/m3) 

Ui  liquid  viscosity  (kg/m-s) 

k,  liquid  '.hermul  conductivity  tW/m-K) 

hfg  latent  heat  of  evaporation  (  J/kgj 

p3  vapor  density  <kg/m3) 

Geometry 


D  condenser  tube  inner  diameter  mi) 


Operating  Conditions- 

G  condenser  mass  flux  (kg/m2  -s) 

x0  condenser  inlet  quality  (-) 

ATsa[  Temperature  difference  between  condensing  fluid 
and  the  secondary  coolant  (K) 
a  magnitude  of  acceleration  (m/s2) 

9  angle  between  acceleration  and  flow  vectors  (°). 
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Calculating  Pressure  Losses.  The  total  pressure  drop  in  the  annular  condensation  region 
is  the  sum  of  three  component  pressure  drops  due  to  friction,  body  forces  and  acceleration. 
Dimensionless  design  maps  for  these  three  pressure  drop  components  are  shown  in  Figures 

3.14,  3.15  and  3.16,  respectively.  All  of  the  dimensionless  pressure  drops  are  displayed  as 
functions  of  quality,  x.  In  most  cases  it  is  sufficiently  accurate  to  approximate  this  quality  by 
the  condenser  inlet  quality.  The  following  sections  explain  the  use  of  the  design  maps. 

1 )  Friction  in  the  Two-Phase  Region.  Estimate  the  frictional  pressure  drop  using  Figure 

3.14.  The  procedure  is: 

•  Calculate  Rc,  the  condenser  thermodynamic  ratio: 


Rc  = 


4k,  AT, 


(3-57) 


Calculate  the  ratio  of  liquid  to  vapor  densities,  p/Pg- 
Calculate  the  Reynolds  number  due  to  liquid  flow,  Re (. 

Re  f  =  GD/jj.  ] 

Calculate  fw ,,  the  liquid-only  friction  factor: 

fwi  =  0.046  Ref0-2 


(3-58) 


(3-59) 


Locate  the  curve  in  Figure  3.14  which  corresponds  to  the  liquid-to-vapor  density 
ratio. 

Locate  the  point  on  this  curve  which  corresponds  to  the  condenser  inlet  quality, 
x0.  The  ordinate  of  this  plot  is  the  dimensionless  pressure  drop  (APfc*/fw  i)- 

Calculate  APfc,  the  frictional  component  of  the  pressure  drop  in  the  annular  part 
of  the  condenser: 


_  lf  rG2/2p,|  rAPfc*l 

ar  lc  -  4  iw  i  jj  — r  i 

L  ‘'c  J  L  *wl  j 


(3-60) 


2)  Body  Forres  in .  (he  Two-Phase  Region.  Estimate  the  pressure  drop  due  to  body 
forces  using  Figure  3.15.  The  procedure  is: 

•  Locate  the  curve  in  Figure  3.15  which  corresponds  to  the  liquid — to — vapor  density 
ratio. 

•  Locate  the  point  on  this  curve  which  corresponds  to  the  condenser  inlet  quality, 
x0.  The  abscissa  of  the  plot  gives  the  dimensionless  pressure  drop  APgc*. 

•  Calculate  APgc,  the  pressure  drop  due  to  body  forces  by: 


APgc  - 


D  a  cosQ 


AP  * 

urgc 


(3-61) 
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Figure  3.16.  DESIGN  MAP  FOR  CONDENSER  ACCELERATION  PRESSURE  DROP 


3)  Acceleration  in  the  Two-Phase  Region.  Estimate  the  pressure  drop  due  to  fluid 
acceleration  using  Figure  3.16.  The  procedure  is: 

•  Locate  the  curve  in  Figure  3.16  which  corresponds  to  the  liquid— to— vapor  density 
ratio. 

•  Locate  the  point  on  this  curve  which  corresponds  to  the  condenser  inlet  quality,  x. 
The  abscissa  of  the  plot  gives  the  dimensionless  pressure  drop  APaN 

•  Calculate  AP1C,  the  pressure  drop  due  to  fluid  acceleration: 


A1TC  =  - 


{j-'Z.'j 


ap; 


v»2) 


4)  Trial  Pressure  Drop  in  the  Two- Pease  Section.  Now  calculate  A?  the  total 
pressure  drop  in  the  two— phase  section  of  the  condenser: 


AP  £  =  APfc  +  APgC  +  APac  (3—63) 

5)  Pressure  Drop  in  the  Sinsle-Phase  Liquid  Region.  Usually  the  pressure  drop  in  the 
single-phase  region  will  be  very  small  due  to  low  fluid  velocity.  It  is  easily  estimated  using 
parameters  which  have  already  been  calculated  to  find  the  pressure  drop  in  the  annular 
condensation  region.  The  pressure  drop  in  the  single-phase  liquid  region,  AP^,  can  be 
estimated  by: 
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(3-64) 


APsp  =  4  fw ,  (Lsp/D)  (G2/2p  ])  -  piL^  a  cosG 

where  fwi  =  the  liquid— only  friction  factor 

(Equation  3—59),  and 

Lsp  =  the  length  of  the  single-phase  region 
(from  Figure  3.6  and  Equation  3-8). 

6)  Total  Pressure  Drop.  The  total  pressure  drop  in  the  condenser  is: 

APcond  =  AP,p  +  APsp  (3—65) 

Uncertainty.  The  major  uncertainty  in  the  condenser  pressure  drop  calculation  is  due  to 
the  interfacial  friction  between  the  vapor  and  liquid  phases  in  the  annular  condensation  regime. 
The  design  maps  in  Figures  3.14  through  3.16  assume  that  the  interfacial  friction  factor  is 
greater  than  the  smooth— wall  friction  factor.  Specifically,  these  curves  use  the  friction  factor 
correlation  of  Waifs  (1969)  which  is  presented  in  Section  1.2  in  this  manual.  This  is  a 
conservative,  high  e-  timate  for  the  interfacial  friction. 

A  lower  bound  for  the  pressure  drop  can  be  estimated  by  assuming  that  the  interfacial 
friction  factor  is  the  same  as  if  the  vapor  were  flowing  in  a  smooth  pipe.  The  friction  pressure 
losses  are  roughly  a  factor  of  two  less  than  would  be  calculated  using  the  design  map  in  Figure 
3.14  with  the  higher  interfacial  friction.  Design  maps  using  interfacial  friction  relationships 
for  microgravity  conditions  can  be  produced  similarly  when  validated  models  become 
available. 


3.2.3  Detailed  Equations  for  Pressure  Losses 

This  section  presents  pressure  drop  equations  which  can  be  used  in  detailed  computer 
calculations  and  explains  the  basis  for  the  dimensionless  design  maps.  Table  3.2  summarizes 
the  equations  for  friction,  body  force  and  fluid  acceleration  pressure  drops  in  the  four  regimes 
of  heat  transfer  within  a  condenser.  The  following  sections  present  more  detail. 


Friction.  The  basic  equation  for  the  two-phase  pressure  gradient  due  to  wall  friction  is 
the  same  as  for  the  case  of  the  evaporator  tube: 

(dP/dz)f  =  — fw !  O2  ( 1-x)2  (4/D)  (G-/2pi)  »3-f>6i 

All  of  these  parameters  are  design  inputs  except  for  the  single— phase  friction  factor.  f.A ,.  and 
the  two— phase  multiplier,  02.  The  friction  factor  is  calculated  using  Equation  3-59.  The 
value  of  Q2  depends  on  the  fluid  properties,  the  quality  and  the  flow  regime  as  described 
below. 

Annular  Condensation  Regime.  The  two— phase  multiplier  for  annular  flow  is  (Collier. 
1986): 

O2  -  rrW  <3H>7) 
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wheie  a  is  the  local  void  fraction.  Then  the  pressure  drop  due  to  friction  in  a  short  tube 
segment  of  length  Az  can  be  calculated  using  this  expression: 

APf  =  ^2  (1— x)2  (4/D)  (G2/2p  [)  Az  (3-68) 

This  is  identical  to  the  expression  for  an  evaporator  tube  (Equation  2—64).  See  Section  2.2.3 
for  a  procedure  to  evaluate  the  void  fraction. 

Note  that  the  symbol  "APf"  is  defined  as  a  pressure  drop: 


APf  = 


(dP/dz)f  dz 


A  positive  value  therefore  means  a  pressure  loss  and  corresponds  to  a  negative  pressure 
gradient.  (A  negative  pressure  drop,  which  occurs  with  momentum  changes  due  to  phase 
change,  indicates  a  pressure  recovery.) 

Other  Condensation  Regimes.  In  the  slug  and  bubbly  condensation  regimes  we 
recommend  a  homogeneous  flow  model  to  evaluate  the  pressure  drop.  The  two-phase 
mul'iplier  and  segment  pressure  drop  are  calculated  using  the  same  expressions  as  for  the 
nucleate  boiling  sections  of  an  evaporator  tube.  Refer  to  Section  2.1.3  for  this  procedure. 

Dimensionless  Formulation.  The  dimensionless  formulation  assumes  that  the  entire 
two-phase  section  of  the  condenser  can  be  treated  as  an  annular  flow.  The  rate  of  pressure 
decrease  with  tube  length  is  convened  to  a  quality  dependence  using  Equation  3-24  which 
relates  the  change  in  quality  to  distance  along  tne  condenser  tube.  The  rate  of  pressure 
decrease  with  fluid  quality,  dP/dx,  can  be  written  as: 

-i=  U-*)>  <4/D)(G=/2pi)  {g\  <3-69a) 

=  -  (l-x)2  (4/D)  (G2/2p  ])  (3— 69b) 

=  ( 1-x):  !'1/rc)  (G2/-P  l>  . - 


When  Equation  (3— 69c )  is  rearranged  and  integrated,  it  becomes 


APf 


:  l-x):  Jx 
1  (l-oo 


where  xs  is  the  quality  at  the  point  of  transition  to  slug  flow,  and  AP*  is  the  dimensionless 

fc 

pressure  drop  due  to  friction: 


AP  fc  f"  R  c 

AP*  =  - - 

fc  4  [G2/2pi 


(3-71) 
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The  dimensionless  map  in  Figure  3.14  has  been  created  by  using  Equations  3—70  and  3-71, 
solving  the  integral  with  the  void  fraction/quality  relationships  shown  in  Table  3.2.  To 
precisely  calculate  the  pressure  drop  using  Figure  3.14,  the  procedure  is  to  take  the  difference 
in  ordinates  which  correspond  to  qualities  x0  and  xs.  However,  the  ordinate  corresponding  to 
xs  is  so  close  to  zero  that  this  correction  makes  no  practical  difference.  It  is  accurate  to 
approximate  (APfc*/fwi)  as  the  ordinate  in  Figure  3.14  corresponding  to  the  inlet  quality,  x0. 


Body  Forces.  The  basic  differential  relation  for  the  pressure  drop  due  to  body  forces  is: 

(dP/dz)g  =  [  a  pg  +  (1-a)  pi  ]  a  cos0  (3-72) 

Evaluation  of  the  pressure  drop  in  an  incremental  pipe  segment  is  thus: 

APg  =  —  [  a  pg  +  (1— a)  pf  ]  a  cos0  Az  (3-73) 

The  void  fraction,  a,  depends  on  the  heat  transfer  regime  as  follows. 

Two-Phase  Forced  Convection  Regime.  In  this  regime  the  void  fraction  should  be 
evaluated  using  the  same  procedure  as  described  above  for  the  case  of  the  friction  pressure 
drop. 


Other  Heat  Transfer  Regimes.  In  the  slug  and  bubbly  regimes  the  flow  is  assumed  to  be 
homogeneous  and  the  void  fraction  is: 


a  = 


Tl-x)(pg/pi)  +  x 


(3-74) 


Dimensionless  Formulation.  The  dimensionless  expression  for  the  pressure  drop  due  to 
body  forces  is  obtained  as  follows: 


dP  dP  dz  r  ,  „  (l—a)D 

_  dx  ~  -  faPg  +  C1-01)  Pi  1  a  cos9  -4TJ-  — 


dz  dx 


"VT, 


=  ■  a c°-~  iodpg/pi)  -  (l-aiRl-a) 


3-75) 


Rearrange  and  integrate  to  obtain: 


\P  *  =  - 
argc 


[a(pg/pi)  +  (1-a)]  (l-a)dx 


(3-76) 


where  APgc*  is  the  dimensionless  pressure  drop  due  to  body  forces  in  the  condenser: 

Rc 

APgc*  -  APgc 


PiDacos© 


(3-77) 


L 
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The  pressure  drop  due  to  body  forces  is  negative  because  when  the  body  force  is  in  line  with 
the  flow  (cosQ  >  0)  there  is  a  pressure  recovery  in  the  condenser.  The  dimensionless  map  for 
the  pressure  loss  due  to  body  forces  (Figure  3.15)  is  obtained  by  assuming  that  the  entire 
condenser  can  be  characterized  by  annular  condensation.  The  design  map  is  generated  by 
using  the  void  fraction/quality  relationship  for  annular  flow  as  given  in  Table  3.2. 


Acceleration.  The  basic  equation  for  pressure  loss  due  to  fluid  acceleration  is: 


(dP/dz)a  =  -G2 


d  fx2  1  (  ( 1  — x )2  1 

3z  a  pg  (1-a)  pi 


Thus  the  pressure  loss  in  a  pipe  segment  due  to  fluid  acceleration  is: 


(3-78) 


APa 


L  a  pg  (1-a)  pi 


(3-79) 


Note  that  this  expression  is  independent  of  the  segment  length  and  depends  only  on  the  fluid 
quality  and  void  fraction  at  the  iniet  and  outlet  of  the  segment.  This  pressure  drop  is  evaluated 
using  the  same  void  fraction/quality  relationships  described  above  for  the  pressure  losses  due 
to  friction  and  body  forces.  The  dimensionless  form  is  obtained  by  dividing  Equation  (3—79) 
by  (G2/2p !)  to  obtain: 


AP 


* 

ac 


=  A 


'5.2£l  +  (}z*)2' 
a  pg  (1-a) 


(3-80) 


where  APac*  is  the  dimensionless  pressure  drop  due  to  fluid  acceleration  in  the  condenser: 


AP  *  =  AP, 


1 


G2/2p, 


=  -APa* 


(3-81) 


3.2.4  Validation  with  Microgravity  Data 

Only  a  few  microgravity  experiments  for  condensation  behavior  have  been  performed. 
These  are  discussed  by  Keshock  (1974),  Lee  (1987),  and  Cuta  and  Krotiuk  (1988).  The  early 
experiments  by  Kesnock  (197-1)  consider  primarily  neat  transfer  data  and  do  not  have  pressure 
drop  data  for  comparisons.  In  the  experiments  performed  by  Lee,  which  are  the  same 
experiments  discussed  bv  Cuta  and  Krotiuk,  tne  data  show  that  equilibrium  conditions  were 
<  not  achieved  in  the  condenser.  Some  unsteady  behavior  was  observed,  therefore  the  data  are 

unsuitable  to  validate  the  design  methods.  Ongoing  experiments  by  the  Weapons  Laboratory 
shouid  provide  the  data  necessary  to  validate  the  recommended  models. 
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GLOSSARY 


Annular.  A  two— phase  flow  regime  in  which  liquid  flows  as  a  continuous  film  on  the  wall  of 
a  pipe  (in  an  annular  pattern)  and  the  vapor  flows  as  a  continuous  gas  core  in  the 
center. 

Bubbly.  A  two— phase  flow  regime  in  which  vapor  flows  as  discrete  bubbles  in  a  liquid. 

Capillary  Condenser.  A  condenser  in  which  the  heat  transfer  rate  is  governed  by  the  thickness 
of  a  liquid  film,  which  is  in  turn  determined  by  capillary  forces. 

Critical  Heat  Flux  (CHF).  The  point  at  which  the  heat  transfer  coefficient  decreases  and  the 
wall  temperature  rises  dramatically  in  boiling  heat  transfer.  This  is  usually  associated 
with  dryout,  or  formation  of  a  vapor  film  at  the  heated  surface. 

Drift— Flux.  An  approach  to  modelling  two— phase  flows  in  which  the  momentum  equation  is 
written  similarly  to  the  homogeneous  modelling  approach,  but  the  velocities  of  each 
phase  are  related  by  a  separate  empirical  relationship  called  the  "drift— flux" 
relationship. 

Flow  Regime.  A  repeatable  pattern  by  which  liquid  and  vapor  flow  in  a  two— phase  mixture. 

Flow  Regime  Transition.  A  boundary  defining  the  flow  conditions  for  a  change  from  one  flow 
regime  to  another. 

Fully  Developed  Subcooled  Boiling.  Boiling  regime  characterized  by  detachment  of  vapor 
bubbles  from  the  wall  surface.  Vapor  bubbles  and  subcooled  liquid  exist 
simultaneously  in  the  heated  flow  channel. 

High— Shear  Condenser.  A  condenser  approach  in  which  the  heat  transfer  is  governed  by  the 
thickness  of  a  liquid  film  whose  dimension  is  set  by  friction  between  the  vapor  and 
liquid  phases. 

Homogeneous.  An  approach  to  modelling  two— phase  Hows  in  which  the  momentum  equation 
is  written  in  a  manner  analogous  to  single-phase  flow,  except  that  the  density,  velocity, 
and  friction  factor  are  quantities  based  upon  average  fluid  properties.  This  model 
assumes  that  the  vapor  and  liquid  phases  travel  at  the  same  actual  velocity. 

Interfacial  Shear  or  Interfacial  Friction.  The  shear  at  a  gas-liquid  interface,  given  by 
empirical  relationships  in  this  Design  Manual. 

Mass  Flux.  The  total  mass  flux  (liquid  plus  vapor)  per  unit  area  of  the  flow  cross-section. 

Onset  of  Nucleate  Boiling.  The  wall  superheat  (temperature  difference  between  wall  and  fluid 
saturation  temperature)  at  which  vapor  bubbles  begin  to  nucleate  in  cavities  of  the  wall 
surface.  The  phenomenon  defines  the  transition  from  single-phase  to  subcooled 
boiling. 

Partially  Developed  Subcooled  Boiling.  Boiling  regime  in  which  vapor  bubbles  form  and 
condense  at  the  heated  wall  surface,  but  remain  attached  to  cavities  in  the  wall  surface. 
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Quality.  The  fraction  of  the  total  mass  flow  rate  which  is  vapor  in  a  vapor— liquid  two— phase 
flow.  Quality  is  obtained  by  taking  the  mass  flow  rate  of  the  vapor  and  dividing  by  the 
total  mass  flow  rate  of  the  vapor  plus  the  liquid. 

Separated.  An  approach  to  modelling  two^jhase  flows  in  which  the  conservation  equations 
are  written  separately  for  each  phase  ai.d  interaction  terms  are  defined  to  relate  the 
shear  (friction)  between  the  phases.  (See  Interfacial  Shear.) 

Slug.  A  two— phase  flow  regime  in  which  vapor  flows  as  bullet— nosed  bubbles  surrounded  by 
a  thin  liquid  film.  The  bubbles  are  typically  several  pipe— diameters  in  length  and 
separated  by  slugs  of  liquid. 

Stratified.  A  two— phase  flow  regime  in  which  vapor  flows  in  the  upper  portion  of  a  pipe  and 
liquid  flows  in  the  lower  portion  of  the  pipe.  Both  phases  are  continuous  in  the  flow 
direction. 

Superficial  Velocity.  A  superficial  velocity  is  calculated  for  each  of  the  liquid  and  vapor 

phases.  The  superficial  velocity  is  the  velocity  that  each  phase  would  have  if  it  flowed 
alone  in  the  pipe,  e.g.  the  superficial  velocity  for  the  liquid  phase  is  j  i  =  Q/A  and  the 
superficial  velocity  for  the  vapor  phase  is  jg  =  Qg/A.  Because  these  superficial 
velocities  are  based  upon  known  input  quantities,  the  phase  flow  rates  and  the 
cross-sectional  flow  area,  two-phase  flow  models  generally  use  these  parameters  and 
relate  them  to  the  actual  flow  velocities  through  the  void  fraction.  For  example,  the 
actual  velocity  for  the  liquid  phase  Uj  =  j/(l  -  a),  and  the  actual  velocity  for  the  vapor 
phase  is  Ug  =  jg/a.  Determining  the  void  fraction,  or  the  relationship  between  the 
velocities  of  each  phase,  is  most  of  the  problem  in  a  two— phase  flow. 

Void  Fraction.  The  volume  of  a  vessel  or  the  cross-sectional  area  of  a  pipe  occupied  by  the 
vapor  in  a  two-phase  flow  system. 

Volumetric  Flux.  See  Superficial  Velocity.  The  superficial  velocity  is  also  know  as  the 
volumetric  flux  because  it  is  the  volumetric  flow  rate  per  unit  area  of  the  flow 
cross-section. 
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